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SUMMARY 


This  report  contains  a  statement  of  work  accomplished  during 
the  past  year  in  the  area  of  chemical  and  molecular  lasers.  In 
formulating  this  program  five  distinct  problems  were  identified  as 
being  basic  to  progress  in  the  further  development  of  existing 
molecular  and  chemical  lasers  and  to  the  discovery  of  new  laser 
systems.  The  problems  identified,  and  the  methods  employed  in  each 
area  will  now  be  summarized: 

1.  Vibrational  Deexcitation  and  Energy  Transfer 

In  every  molecular  and  chemical  laser  the  lasing  molecules  are 
subjected  to  collisions  with  each  molecular  constituent  of  the  gas 
thereby  providing  a  source  of  loss  of  vibrational  energy.  In 
addition,  in  certain  lasers,  e.g.,  DF-CO^  collisions  between  molecules 
and  the  transfer  thereby  of  vibrational  energy  from  one  molecule  to  another 
is  an  efficient  mechanism  for  achieving  specific  molecular  excitation. 

In  most  known  and  important  molecular  lasers  many  important  rates  of 
deexcitation  and  energy  transfer  were  unknown.  As  a  result  of  our 
effort  we  have  experimentally  determined  many  of  these  rates. 

a.  The  method  of  laser  induced  fluorescence  was  applied  to 
study  these  rates  in  the  KF,  DF,  HF-C02  and  DF-C02  chemical 
lasers.  As  a  result  of  this  work  the  rate  of  vibrational 
deexcitation  of  DF;  the  rates  of  energy  transfer  from  HF 
and  DF  to  C0_,  and  the  rates  of  deexcitation  of  C0o  in  the 

4  tm 

upper  laser  level  by  HF  and  DF  have  been  measured. 

b.  Shock  tube  studies  in  C02  were  devoted  to  the  measurement  of 
vibrational  relaxation  in  C02  induced  by  a  selected  sequence 
of  molecular  collision  partners.  These  partners  were  selected 
to  elucidate  the  effects  of  mass,  moment  of  inertia,  dipole 
moment,  and  energy  mismatch  between  nearly  resonant  energy 
levels  upon  the  deexcitation  efficiency.  Collision  partners 
with  C02  were  H20,  D20,  CH^,  CD^  and  CH^F.  The  rates  were 
measured  from  300°K  to  1200°K  and  were  correllated  with  prior 


work  where  available  and  with  theory, 
c.  Work  was  initiated  to  determine  the  collisional  deactivation 
of  CO2  in  the  upper  laser  level  by  itself  and  other  molecules 
over  a  pressure  range  1  -  100  atmospheres  and  a  temperature 
range  300°  -  600°K.  This  work  is  motivated  by  the  recent 
advent  of  CO2  lasers  operating  at  as  high  a  pressure  as 
1  atmosphere  and  by  the  desirability  of  developing  lasers  that 
operate  at  even  higher  pressures  to  maximize  volumentric 
energy  storage.  Relevent  rates  at  and  above  one  atmosphere 
are  simply  not  known.  The  method  employed  is  again  laser 
induced  fluorescence  for  which  we  have  developed  a  chemical 
HBr  laser  with  relatively  high  power  output  (100  -  1000  watts) 
which  matches  several  strong  absorptions  in  CO2  and  ^0. 

This  work  is  continuing. 

2.  Atom  Production  for  Chemical  Laser  Initiation 

Most  known  chemical  lasers  require  the  presence  of  specific  free 

atoms  to  act  as  Initiators  of  and  participants  in  the  chemical  reactions 

from  which  laser  excitation  is  derived.  It  is  therefore  desirable  to 

establish  methods  of  copious  and  efficient  atom  production.  Atoms  are 

frequently  produced  by  electrical  discharges  of  various  types  in  a 

molecular  gas  whereby  molecular  dissociation  upon  electron  Impact 

provides  one  or  more  atoms.  A  few  cases  are  known  in  which  purely 

chemical  means  may  be  employed  for  atom  production. 

a.  We  have  studied  atom  production  by  various  gas  discharge 

me.>ns  and  have  found  that  a  microwave  discharge  can  be 

operated  with  good  results.  Using  a  1500  watt  cw  magnetron 

operating  at  2.45  GHz  we  have  obtained  up  to  50%  dissociation  of  O2 

at  pressures  of  5-6  Torr  with  the  addition  of  an  optimum  amount 

of  N».  At  4  liters  per  sec  pumping  speed  and  2.3  Torr  the  0  atom 
‘  20 

production  rate  was  1.4  x  10  atoms/sec.  Thus  production 
rate  was  basically  limited  by  pumping  capacity  and  could 
easily  be  increased  by  more  than  an  order  of  magnitude.  At 
these  flow  rates  of  atomic  oxygen  it  becomes  possible  to  design 
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chemical  lasers  such  as  the  CS2  +  0  system  that  produces 
lasing  in  CO  at  the  Kilowatt  level,  cw.  We  are  considering 
the  parametric  study  of  such  a  laser  which  could  be  readily 
undertaken  with  the  atom  generation  means  we  have  developed. 

b.  One  of  the  better  known  purely  chemical  means  of  atom 

production  involves  the  NO  +  F2  NOF  +  F  reaction  by  which 
atomic  fluorine  is  produced.  This  reaction  has  been  applied 
to  the  purely  chemical  DF  -  C02  vibrational  transfer  chemical 
laser.  The  rate  at  which  the  reaction  proceeds  is  not  known 
accurately  and  we  have  begun  a  study  of  it  utilizing  the 
method  of  paramagnetic  resonance  absorption  to  directly 
measure  density  of  F  and  NO  as  the  reaction  proceeds. 

This  work  is  continuing. 

3.  Detailed  Study  of  Chemical  Laser  Systems 

Two  chemical  laser  systems  have  been  studied  in  detail  utilizing 
parametric  variation  of  reactants  and  diluents  as  well  as  spectroscopic 
study  of  emission  wavelengths,  populations  and  time  sequence  of  particular 
emissions. 

a.  The  highly  energetic  reaction  between  aceteline  (C^)  and 
atomic  oxygen  to  yield  vibrationally  excited  CO  in  a  laser 
configuration.  The  results  obtained  indicate  that  a  high 
vibrational  temperature  characterizes  the  CO  produced  and 
that  under  normal  operating  conditions  a  relatively  high 
rotational  temperature  in  excess  of  600°K  also  occurs. 

Addition  of  diluent  decreases  the  rotational  temperature  and 
it  would  appear  that  effective  gas  cooling  should  result  in 
laser  action.  Computer  fitting  of  the  spectral  intensity 
data  obtained  is  in  progress  and  will  result  in  knowledge 

of  vibrational  populations  in  the  chemically  produced  CO. 

b.  A  study  of  the  CS2  +  0  reaction  which  also  leads  to  lasing 
in  CO  is  in  progress.  The  object  in  this  study  is  to 
elucidate  the  mechanisms  that  operate  in  this  reaction.  A 
chemical  CO  laser  utilizing  this  reaction  was  operated -and 
relative  laser  emission  intensity  as  a  function  of  vibrational 
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transition  was  studied  as  were  the  time  delays  after 
initiation  of  the  reaction  for  the  different  vibrational 
transitions.  At  present  non  lasing  chemiluminescent  spectral 
information  is  being  studied.  In  combination  these  data 
will  permit  an  unambiguous  formulation  of  a  mechanism  and 
the  rate  constants  to  be  determined. 

4.  Reaction  Rates  for  Chemical  Lasers 

Many  chemical  reactions  that  are  employed  in  chemical  lasers  have 
unknown  or  poorly  known  rates.  These  rates  must  be  determined  in  order 
to  permit  a  computer  model  for  the  laser  to  be  developed  and  thereby  to 
obtain  reliable  prediction  and  analysis  of  laser  performance.  We  have 
undertaken  to  measure  these  unknown  and  important  rates  by  using  the 
method  of  microwave  electron  paramagnetic  resonance  (EPR)  to  monitor 
the  presence ,  appearance  and  disappearance  of  the  paramagnetic  free 
atoms  and  molecules.  A  fast  flow  vacuum  system  coupled  to  a  high 
sensitivity  GPR  spectrometer  was  assembled  and  means  for  gas  Injection 
suitable  for  chemical  rate  determination  was  designed.  Rate  information 
on  the  reactions  H  +  F£  ■*  HF  +  F  and  H  +  HC1  +  Cl  has  been 

obtained  and  the  teaction  NO  +  Fg  *►  NOP  +  F  is  being  studied.  These 
rates  are  basic  to  an  understanding  of  the  operation  of  HF,  HC1,  and 
DF  -  CO,j  purely  chemical  lasers. 

5.  Turbulent  Induced  Rapid  Gas  Mixing 

Many  important  molecular  and  chemical  lasers  require  rapid  and 
efficient  gas  mixing  into  flowing  gas  streams.  Introduction  of  controlled 
turbulence  into  the  stream  of  mixing  gases  promises  to  permit  faster  gas 
mixing  than  simple  diffusive  mixing.  Reliable  experimental  and  theoretical 
design  information  is  lacking  at  the  subatmospheric  pressures  characteristic 
of  many  laser  systems.  We  have  undertaken  an  experimental  and  theoretical 
study  of  turbulent  mixing  with  the  goal  of  determining  reliable  design 
information  at  the  pressures  of  interest.  A  mixing  wind  tunnel  was  built 
and  turbulent  mixing  studied  by  the  introduction  of  posts,  grids,  and 
injectors.  The  enhancing  effect  of  turbulence  on  the  mixing  was  observed. 
Further  tests  including  more  extensive  variation  of  operating  parameters 
will  supply  the  required  design  characteristics  for  effective  turbulent 
mixing. 
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VIBRATIONAL  DEEXCITATION  AND  ENERGY  TRANSFER 
IN  THE  HF,  DF,  HF-C02,  and  DF-C02  CHEMICAL  USERS 

T.  A.  Cool 

Much  progress  has  been  achieved  recently  In  the  development  on 
continuous-wave  (cw)  chemical  lasers  despite  a  serious  lack  of  under¬ 
standing  of  the  mechanisms  of  those  vibrational  energy  transfer 
processes  which  play  major  roles  In  chemical  laser  operation. 

Prior  to  the  completion  of  the  measurements  described  here, 

reliable  data  were  available  for  the  self  deactivation  of  HF  as 

1  2 
reported  by  Alrey  and  Fried  ,  and  more  recently  by  Solomon,  et  al  , 

3 

and  Bott  and  Cohen  .  No  data  were  available,  however,  concerning 
the  rate  of  vibrational  deactivation  of  DF  or  the  rates  for  the 
Important  processes  by  which  vibrational  energy  is  transferred  from  HF 
and  DF  to  the  C02(001)  upper  laser  level.  In  addition  the  rates  for  the 
vibrational  deactivation  of  C02(00°1)  molecules  by  collisions  with  HF 
and  DF  molecules  were  unknown. 

He  have  used  the  laser  excited  fluorescence  method  ^  to  determine 
the  rate  constants  for  these  key  processes  at  350°K,  a  typical  temperature 
for  chemical  laser  operation.  These  rates  are  given  below: 

(a)  DF  (v-1)  +  C02(000)  DF  (v-0)  +  C02(001);  k  =  17.5  +  2.5  x  104  sec"1  torr' 

(b)  DF  +  CO, (001)  -►  DF  +  CO, (000);  k  ■  1.9  +  0.4  x  104  sec"1  torr"1 

L  1  4-1-1 

(c)  DF  (v-,1)  +  DF  (v-0)  •+•  2DF  (v-0);  k  -  2.0  +  0.2  x  10  sec  torr 

(d)  HF  (v— 1)  +  C02(000)  -*•  HF  (v-0)  +  C02(001);  k  -  3.7  +0.3  x  104  sec’1  torr" 

(e)  HF  +  CO, (001)  ■+  HF  +  CO, (000);  k-3.6  +  0.3xl04  sec"1  torr"1 

L  L  4-1-1 

(f)  HF  (v— 1)  +  nf  (v-0)  -*•  2HF  (v-0);  k  -  5.25  +  0.30  x  10  sec  torr 

The  experimental  results  for  HF  and  DF  vibrational  deactivation 
exhibit  good  qualitative  agreement  with  the  vibration-rotation  theory  of 
Shin,  when  the  effects  of  strong  attractive?,  interactions  are  considered7. 

The  rates  of  vibrational  energy  transfer  from  HF  and  DF  to  C02  are 
too  large  to  be  accounted  for  with  present  approximate  theories.  Interpreta¬ 
tion  of  the  above  results^  suggests  that  the  magnitude  of  these  rates  may 
be  explained  in  terms  of  a  hydrogen  bonded  attractive  intermolecular  potential 


well  for  HF-CO2  and  DF-COj  collision  pairs  which  leads  to  enhanced 
repulsion  at  close  ranges.  A  successful  theory  must  necessarily 
allow  the  major  portion  of  the  vibrational  energy  defect  between 
initial  and  final  states  to  be  taken  up  as  rotational  motion  of  the 
collision  partners. 

The  approximations  Inherent  in  the  first  order  perturbation  theory 
of  multipole  moment  interactions  are  invalid  for  treatment  of  the 
energy  transfers  to  CC^OOl)  from  HX  (v-1)  and  DX  (v»l)  molecules. 

We  have  been  successful,  however,  in  calculation  of  the  previously 

8  6 
measured  rate  of  energy  transfer  between  D£  (v«l)  and  00^(001) . 

For  this  case  an  Interaction  between  the  quadrupole  moment  of  ortho- 

deuterium  and  the  transition  dipole  moment  of  the  CO^  asymmetric 

stretching  vibration  is  responsible  for  the  observed  rate  of  transfer. 

Our  apparatus  is  currently  being  used  for  measurements  of 
V  -*■  V  coupling  between  HF  and  Hj  and  DF  and  D2.  We  plan  on  detailed 
studies  of  the  above  processes  over  the  temperature  range  200  -  700°K 
during  the  coming  year  to  further  elucidate  the  energy  transfer 
mechanisms . 

A  brief  description  of  the  apparatus  we  have  used  is  given  in 
reference  4  which  is  also  included  in  the  appendix  to  this  report. 
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SHQCK  TUBE  STtJ&IES  , 
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S.  H.  Rauer  ,  , 
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A  substantial  effort  has  been  devoted  to  the  study  of  vibrational 

i  j  i  ii  1 

relaxation  In  COj  Induced  by  a  selected  sequence  of  colliders.  Density 
gradients  at  shock  fronts  were  measured  with  a  laser-schlieren  arrangement. 

,  ■  ,  '  i  '  ■  ;  |  1  :i'  ' 

These  data  provide  vibrational  ■'•translation  relaxation  times  for  the 

1  i 

(010)  state.  The  experimental  set,  up  is  shown  in  Figure  5;  operating 
conditions  are  listed  helper.  , 


SHOCK  TUBE  OPERATING  CONDITIONS 


Shock  speed :  /  measured  to  +1  mm/psec 


raster  accuracy  +  0.5  nsec 


} 


'  I 

Background  pressure:'  7  x  10  torr 


Leak  rate:  4  x  10”^  torr /min 


^unrelax!*5  ±  20'K 


sample  pressures 
(5-150  torr) 


Degassing  minimized  by  flushing  tube  with  sample 

CO.  —  research 'grade  [N.  (7.4  ppm);'  0.  (1.7  ppm);  others  {Os  1  ppm)] 

4  <  *  ,  '  *  i  ,  t1 

Estimated  Impurity  level  (lncl.  degass)  ft*  20  ppm  i 

P2O—  99.82%  level;  other  &ases  of  high  quality 

Driver  gases f  He  and  N.  (100  -1  200  psi)  '  , 

■  “  ,  1 

Rise  time  of  optical-electronic  system:  j  is  200  n  sec 

I  ,  1 

Correction  for. finite  rise  tim^  inserted  in  program  1 

t  I  »  f  I 

The  most' difficult  problem  encountered  in  this  investigation,  when  H^O 

and  DjQ  served  as  catalyst  gases  Was  the  determination  of  the  actual 

water  content  in  the  shocked  gas.  After  testing  several  analytical 

devices ;Mr.  M.  Buchwald  finally 'obtained  consistent  results  with  a 

1  1  '  '  1 

vapor  phase  chromatograph.  He  calibrated  the  system,  which  included 

flow  controls,  a  portion  of  tube  identical  with  that  used  in  the  1 

construction  of  the  shock  tube,  and  , the  handling  procedures  by  means 

of  the  arrangement  shown  in  Figure  6.  His  results  are  summarized  In 


1  1 
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Shock- front  Density  gradient  (Schlieren  Optical  Arrangement) 
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Figure  3,  both  for  CC^  -  1^0  and  CO2  -  D2O.  Comparisons  are  also 
made  in  this  figure  with  selected  data  taken  from  the  literature. 

These  experimental  results  for  the  relaxation  of  the  (010)  state  of 
CO2  with  H2O  and  D2O  as  collision  partners,  were  replotted  in  Figure  4 
for  comparison  with  predictions  based  on  the  theory  of  Sharma,  the 
expression  proposed  by  Shields  (for  t^O  and  D2O)  and  the  theory  of 
Widom-Bauer  for  ^0.  The  effect  of  isotopic  substitution  is  best 
reproduced  by  the  Shields  equation  as  is  also  the  general  trend  with 
temperature.  The  precise  magnitudes  for  the  relaxation  times  are 
subject  to  adjustment  with  the  a  parameter.  Similar  shock  front  density 
gradient  data  for  mixtures  of  CO2  with  CH^ ,  CD^  and  CH^F  were  obtained 
by  Mr.  Peter  Walsh;  these  are  reproduced  in  Figures  5,6,  and  7.  As 
an  overall  check  on  the  technique  Mr.  Walsh  also  determined  the 
vibrational  relaxation  time  in  pure  CH^F  as  a  function  of  temperature. 

As  can  be  seen  from  his  curve,  Figure  8  his  values  extrapolate  to 
the  room  temperature  value  report  by  Lambert,  based  on  sound  dispersion 
measurements.  Finally,  these  results  along  with  appropriately  selected 
values  from  the  literature  were  summarized  in  a  series  of  plots.  For 
each  catalyst  molecule  the  de-excitation  probability  for 
C(>2  (010)  +  M  -*■  C02  (000)  +  M  and  the  vibrational  de-excitation 

C02  (001)  +  M  -*■  C02(m,n,0)  +  M  were  compared  (i)  to  ascertain  which 
collision  partner  would  be  most  suitable  for  facilitating  laser  operations, 
such  that  the  lower  lasing  state  is  preferentially  de-exclted,  and 
(11)  to  call  attention  to  intrinsic  molecular  parameters,  such  as  moments 
of  inertia,  dipole  moments,  etc.,  for  the  various  collision  partners 
which  would  provide  essential  parameters  for  any  theory  which  accounts 
for  the  different  efficiencies  for  energy  transfer. 

In  another  shock  tube  experiment  Mr.  R.  Jones  is  measuring  the 
(010)  relaxation  time  in  pure  CO2  as  a  function  of  the  translation 
temperature,  in  the  low  temperature  regime.  His  tube  can  be  cooled  to 
about  -40 °C.  By  now  the  serious  problems  of  reliable  temperature 
control  and  clean  operation  appear  to  have  been  solved;  a  series  of 
runs,  from  about  600°K  and  down,  are  being  made. 


The  discovery  by  Hr.  Miller  that  on  shock  heating  a  mixture  of 

N^O  with  CO  to  a  temperature  around  2500°K,  CC^  is  produced  in  an 

inverted  population,  as  demonstrated  by  gain  in  an  incident  beam  of 

10.6  p  radiation,  is  most  intriguing.  If  confirmed,  this  is  indeed 

a  very  novel  process.  To  date  three  experiments  under  varying 

conditions  indicate  some  degree  of  gain.  Also,  more  than  a  half 

dozen  experiments  have  been  run  with  ^0  and  CO2,  where  no  gain  was 

observed.  The  tentative  assumption  is  that  at  these  high  temperatures 

the  N^O  dissociates  rapidly  to  produce  vibrationally  excited  N2  and 

0(*D);  that  these  (^D  atoms  combine  with  CO  to  produce  C0«(*B_),  which 

1  +  t  l 

then  must  get  de~excited  to  the  C0.(  E  )  by  collision  with  another 
molecule.  The  latter  process  leads  to  preferential  population  of 
the  (001)  state.  This  is  consistent  with  a  paper  published  in  1967  by 
Clerc  and  Barat  on  the  reaction  between  CO  and  0(^D)  for  which  they 
report  a  very  large  association  rate  constant.  The  surprising  feature, 
of  course,  is  that  the  product  COg  is  in  a  non-Boltzmannian  distribution. 
Further  work  is  in  progress. 

Very  briefly  I  wish  to  mention  our  pre-occupation  with  chemical 
effects  Induced  by  laser  irradiation.  The  experiments  to  date  on  the 
reaction  between  NO^  and  CO  to  produce  NO  -I-  CO^  have  led  nowhere.  At 
present  we  are  convinced  that  no  augmentation  of  the  reaction  by  C02 
laser  radiation  occurs.  This  ends  (I  hope)  a  long  list  of  inconclusive 
experiments,  where  we  finally  traced  the  difficulties  to  variations  in 
the  method  of  analysis.  We  may  return  to  it  after  we  set  up  a  laser  which 
will  provide  a  much  closer  coincidence  between  the  incident  and  the 
absorption  spectrum.  In  order  to  select  appropriate  systems  for  such 
studies  we  set  up  a  spectrophone,  as  shown  in  Figure  13.  This  will 
permit  us  to  determine  relatively  easily  which  substances  have  strong 
absorptions  for  particular  lines  in  the  C02  spectrum.  Following  such 
scans  we  should  be  in  a  good  position  to  select  those  species  which  are 
good  candidates  for  selectively  induced  reactions  based  on  their  high 
absorption  and  long  vibrational  translation  relaxation  times. 
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RELAXATION  TIMES  AT  ROOM  TEMPERATURE  FOR, 
C02(001)  +  AB(v=0)  -*  CO2(001)  +  AB(v=l) 

J.  C.  Stephenson,  AVCO  Corporation  • 
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ATOM  SOURCES 

Y.  S.  Liu,  R.  A.  McFarlane  and 
G.  J.  Wolga 

I 

To  provide  atomic  species  for  use  in  laser  systems  employing  exothermic 
reactions  of  the  type 

A  +  EC— »AB*  +  C 

the  molecules  is  usually  dissociated  in  an  electrical  discharge  source  con¬ 
sisting  of  one  of 

|  / 

(a)  C.  W.  direct  current  discharge  between  internal  metal  electrodes. 

(b)  High  current  capacitor  discharge. 

(c)  Microwave  discharge  in  a  resonant  cavity. 

I 

i  (d)  Electrodeless  radio- frequency  discharge. 

Where  relatively  small  flow  rates  of  atomic  hydrogen,  nitrogen,  oxygen  or  the 
several  halides  are  required  any  of  these  sources  is  adequate. 

In  general  however,  d.c.  discharges  provide  a  smaller  percentage  of  dis¬ 
sociation  and  can  give  rise  to  gas  contamination  because  of  the  presence  of 
metal  electrodes.  The  high  electric  fields  employed  in  capacitor  discharge 
sources  cab  be  beneficial  for  dissociating  molecular  species  requiring  larger 
amounts  of  energy  such  as  N2  but  the  flow  of  atoms  is  not  continuous  in  fast 
flow  systems  except  under  high  repetition  rate  conditions. 

f 

The  efficiency  of  atom  production  can  be  very  high  for  radio-frequency 
sources,  and  generators  operating  at  10  -  30  MHz  with  kilowatt  power  levels 
have  been  applied  with  success.  Efficient  coupling  of  the  available  energy 

to  the  gas  column  is  at  best  difficult,  particularly  under  varying  conditions 

! 

of  gas  pressure  and  excitation  level.  Impedance  mismatch  between  the  r.f. 

I 

generator  and  the  gas  column  can  give  rise  to  substantial  radiation  of  energy 
at  the  operating  wavelength  and  high  intensity  r.f.  fields  can  seriously  in- 

i 

fluence  the  qpepation  of  optical  detectors  and  electronic  signal  processing 

I 

apparatus  being  used  in  the  vicinity. 

I 

Similarly  high  efficiencies  of  atom  production  can  be  achieved  with  micro- 

i  , 

wave  excitation  and  the  Evenson-Broida  cavity  driven  by  an  80Wc.w.  S  band 
magnetron  has  been  widely  applied.  This  type  of  source  has  the  advantage  of 
inherently  clean  operation  —  a  1  cm.  diameter  quartz  tube  carrying  the  gas 
passes  through  the  cavity  —  and  is  capable  of  operating  over  a  wide  pressure 
range.  Early  measurements  which  we  made  on  such  a  source  for  atomic  oxygen 
production  indicated  that  in  a  fast  flow  system  the  rate  of  atom  production 


i 


// 


was  continuing  to  rise  at  the  maximum  available  microwave  operating  power. 

The  availability  of  relatively  inexpensive  C.W.  magnetrons  capable  of 
producing  powers  in  excess  of  1  kw.  at  2450MHz  led  to  the  discharge  source 
shown  in  Figure  1.  The  microwave  power  tube  is  a  Litton  Type  L-5001  normally 
employed  for  microwave  heating  and  cooking.  This  tube  is  E-probe  coupled 
to  a  section  of  3"  X  bg"  C  band  waveguide  and  an  adjustable  shorting  plunger 
is  arranged  to  present  the  appropriate  high  impedance  plane  at  the  probe 
position.  The  waveguide  size  employed  is  smaller  than  normally  used  at  2,453  MHz 
but  the  cutoff  wavelength  of  15.2  cms.  is  more  than  adequate  for  an  operating 
wavelength  of  12.2  cms. 

A  100  c.f.m.  blower  provides  cooling  for  the  anode  fin  structure  and 
cathode  terminals.  The  required  heater  power  depends  upon  the  operating  power 
level  because  of  back  bombardment  of  the  cathode.  We  have  found  that  for  D.C. 
anode  currents  between  500  ma.  and  750  ma.  a  voltage  of  4.3  volts  as  measured 
at  the  magnetron  terminals  and  a  heater  current  of  18  amperes  is  adequate. 

High  voltage  is  supplied  by  a  full  wave  three  phase  bridge  rectifier  power  supply 
and  no  filtering  is  required  as  the  inherent  ripple  on  such  a  system  is  only 
about  4$.  The  rated  power  output  to  a  flat  load  for  anode  currents  of  500- 
750  ma.  is  1300-1700  watts.  No  direct  measurements  have  been  made  of  the  power 
delivered  by  the  source  to  the  gas  column  during  operation.  The  absence  of 
a  substantial  temperature  rise  of  the  fin  structure  suggests  that  70-80$  of 
the  available  power  is  deposited  in  the  gas. 

For  initial  studies  of  this  source  a  quartz  discharge  tube  2  cm.  I.D.  was 
passed  through  the  broad  face  of  the  waveguide  section  at  an  angle  of  10°  to 
the  propagation  direction.  A  microwave  noise  source  waveguide  section  was 
used  for  these  experiments  and  the  quartz  tube  was  simply  substituted  for  the 
fluorescent  lamp  noise  generator.  The  waveguide  section  was  closed  by  a 
second  adjustable  shorting  plunger. 

During  operation  the  high  voltage  is  increased  until  the  indicated  anode 
current  is  within  the  range  discussed  above.  For  gas  pressures  in  the  range 
of  1  Torr  a  discharge  will  have  been  initiated  by  this  time  in  the  quartz  tube, 
the  exact  point  of  onset  being  somewhat  dependent  upon  the  position  of  the  second 
shorting  plunger.  This  occurs  of  course  because  in  the  absence  of  gas  break¬ 
down  the  waveguide  structure  can  be  tuned  as  a  resonant  cavity  near  to  the 
operating  frequency.  The  attendant  large  E  fields  then  lead  to  breakdown 


FIGURE 
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in  the  quartz  tube.  The  plunger  position  is  not  critical  and  once  the  discharge 
is  operating,  variation  in  its  position  has  no  effect.  It  is  concluded  there¬ 
fore  that  little  or  no  microwave  energy  is  transmitted  beyond  the  discharge 
tube  after  breakdown  takes  place.  The  angle  of  10°  is  sufficiently  small  that 
microwave  power  reflected  from  the  discharge  column  back  toward  the  magnetron 
should  be  unimportant.  The  only  adjustment  required  therefore  is  the  position 
of  the  first  shorting  plunger  near  the  magnetron  and  it  is  simply  set  for 
maximum  gas  column  brightness.  Once  set  it  needs  no  further  adjustment. 

Measurements  were  made  to  determine  if  leakage  of  microwave  energy  from 
such  a  high  power  source  could  present  a  hazard  to  any  nearby  experimenter. 

Where  the  discharge  tube  leaves  the  waveguide  structure,  radiation  levels  of 
10  mw/cm^  were  detected,  falling  rapidly  to  less  than  1  mw/cm2  at  a  distance 
of  six  to  eight  inches.  A  second  waveguide  structure  described  below  employs 
cylindrical  metal  tube  sections  beyond  cutoff  surrounding  the  quartz  tube  and 
even  lower  levels  of  radiation  were  found  for  that  system.  It  is  concluded 
that  no  hazard  exists  for  normal  operating  procedures  but  the  experimenter  should 

O 

be  aware  that  10  mw/cm  is  an  upper  limit  to  avoid  cornea  damage  and  some 
caution  should  be  exercised  when  viewing  the  discharge  column  at  very  close 
distances. 

To  evaluate  the  performance  of  the  discharge  source,  measurements  were 
made  of  the  dissociation  of  molecular  oxygen  at  pressures  up  to  2.3  Torr  using 
initially  a  relatively  low  capacity  pumping  system.  Gases  were  introduced 
into  the  discharge  tube  via  Ma  the  son  Type  602  flowmeters  and  exhausted  by  a 
Precision  Scientific  pump  with  a  rated  speed  of  6  litres/sec  at  1  micron. 

Atomic  oxygen  concentration  was  measured  at  a  point  50  cm.  from  the  discharge 
region  by  titration  with  NOg.  The  atomic  oxygen  flow  rate  is  equal  to  the  NOg 
flow  rate  required  to  cause  the  disappearance  of  the  green  fluorescence  pro¬ 
duced  as  a  result  of  the  second  reaction  below  - 

0  +  NOg — >  NO  +  0g  (1) 

NO  +  0  — >N0g— *N0g  +  hv  (2) 

At  that  time  reaction  (1)  has  removed  all  atomic  oxygen  and  none  is  available 
for  reaction  (2). 

It  is  known  that  the  presence  of  helium  or  trace  amounts  of  nitrogen  can 
strongly  influence  the  rate  of  atomic  oxygen  production  and  this  is  demonstrated 


in  Figure  2  where  the  production  rate  of  atomic  oxygen  is  shown  as  a  function 
of  the  flow  rate  of  molecular  oxygen.  Figure  3  shows  that  the  pumping  rate 
was  essentially  constant  over  the  pressure  range  employed  and  Figure  4  shows 
the  percent  dissociation  observed  as  a  function  of  molecular  oxygen  flow  rate. 

The  maximum  atom  production  rate  was  1.4  X  10  atoms/sec.  at  a  molecular 
ojQrgen  pressure  of  2.3  Torr  and  approximately  4  litres/sec.  punning  speed. 

For  a  fixed  pressure  this  will  scale  directly  with  pumping  speed  and  a  two 
order  of  magnitude  increase  should  easily  be  obtained.  The  percentage  of 
dissociation  approached  20$  at  the  higher  pressures  used.  The  excitation 
source  will  operate  up  to  pressures  of  10  Torr  and  additional  increase  of  atom 
production  rate  should  be  possible  at  higher  pressures  than  shown  in  Figure  3* 
However  recombination  effects  will  have  to  be  considered  as  the  pressure  is 
increased  and  it  is  to  be  expected  that  an  optimum  pressure  will  be  determined 
for  any  particular  system. 
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NEW  LASER  SYSTEM  STUDIES  -  +  0 

Y.  S.  Liu,  R.  A.  McFarlane  and 
G.  J.  Wolga 


The  capability  of  producing  large  flow  rates  of  atomic  oxygen,  makes 
possible  the  exploration  of  reactions  which  are  sufficiently  exothermic  to 
allow  the  formation  of  vibrationally  excited  products  which  correspond  to 
very  high  vibrational  temperatures.  If  the  rotational  to  vibrational  tem¬ 
perature  ratio  can  be  maintained  adequately  small,  partial  inversion  and 
laser  oscillation  can  be  established  with  the  excited  products.  A  par¬ 
ticularly  attractive  candidate  for  the  production  of  vibrationally  excited 
carbon  monoxide  is  the  reaction  of  oxygen  with  acetylene.  Recent  measure- 

(l  g] 

ments  have  been  reported  on  infrared  chemiluminescence  for  this  system'  * 
at  relatively  low  pressures  (20-300  microns).  The  primary  reactions  are 

0  +  CgHg . ■-■>  CO  +  CHg  +  51  kcal./mole  (1) 

0  +  CHg - )C0  +  2H  +  71  kcal./mole  (2) 

for  low  fuel  flows,  and  excitation  of  CO  up  to  v'  =  l4  is  observed.  With 
increased  fuel  flow  or  the  addition  of  H  atoms,  changes  in  vibrational  distri¬ 
bution  above  v'  =  8  occur  and  emission  corresponding  to  v'  =33  has  been 
identified.  ^  This  is  ascribed  to  the  reaction 

0  +  CH - »C0  +  H  +  176  kcal./mole  (3) 

and  is  clearly  an  additional  potential  source  process  for  populating  higher 
vibrational  levels  for  utilization  as  initial  laser  states. 

Measurements  of  the  relative  populations  of  the  vibrational  levels  of 

(2) 

the  CO  molecule  which  have  been  reported  in  the  literature'  '  indicate 
that  at  low  pressures  (25  microns)  the  population  of  levels  up  to  v'  =  13 
can  be  characterized  by  a  Boltzmann  distribution  with  a  vibrational  temp¬ 
erature  of  approximately  7500°K.  Above  v'  =  13  the  population  fell  off 
more  rapidly  with  vibrational  quantum  number  but  could  be  enhanced  by 
operation  at  higher  pressures  (150  microns).  The  result  can  be  under¬ 
stood  in  terms  of  a  preferential  V-V  pimping  mechanism  examined  by  Rich, 

Watt  and  Thomson  and  shown  to  be  of  crucial  importance  in  electrically 
excited  CO  lasers.  a  striking  increase  in  the  population  of  higher  vib- 


rational  levels  was  observed  by  Thrush  and  co-workers ^  with  fuel  rich  mix¬ 
tures  and  this  can  be  interpreted  as  a  direct  consequence  of  V-V  pumping  at 
increased  levels  of  CO*  concentration.  The  data  for  fuel  lean  and  fuel  rich 
operation  were  analysed  by  Thrush  to  show  that  in  the  latter  case  a  non 
equilibrium  distribution  of  population  was  present  in  CO*  for  vibrational 
levels  higher  than  v'  =  13 •  This  corresponds  directly  to  the  appearance 

of  the  plateau  in  the  N  ,  versus  v'  curve  in  the  analysis  of  Rich,  Watt 
(3)  V 

and  Thomson. 


If  it  were  possible  to  obtain  a  sufficiently  high  "local"  vibrational 

temperature  characterizing  adjacent  vibrational  levels  while  maintaining 

a  low  rotational  temperature,  laser  oscillation  could  be  obtained  in  the 

chemically  formed  CO  species.  Vibrational  temperatures  near  8000°K  for 

(3) 

v'  =8,9  have  been  determined  for  the  electrical  CO  laserVJ'  increasing  to 
greater  than  20, 000°K  for  vibrational  levels  above  v'  =  13.  For  laser 
oscillation  near  5*2  microns  (v1  =  8,9)  the  required  vibrational  temper¬ 
atures  are  apparently  available  from  the  0  +  CgHg  reaction.  There  are  how¬ 
ever  other  temperatures  which  must  be  examined  and  controlled. 


The  V-V  pumping  rates  decrease  with  increasing  translational  temper- 
(3) 

ature  and  losses  due  to  V-T  processes  increase. It  will  be  necessary  there¬ 
fore  to  operate  at  the  lowest  possible  translational  temperatures  in  such  a 
laser  and  it  is  expected  that  the  addition  of  helium  can  be  used  to  lower 
the  translational  and  rotational  temperatures  by  increasing  the  rate  of  heat 
conduction  to  the  tube  walls.  While  C0*-He  collisions  will  contribute  addi¬ 
tional  V-T  losses  the  effect  is  expected  to  be  small  in  comparison  with  the 
enhancement  due  to  the  lowering  of  T,  and  T  , 

uZ*culS«  1*0 Xf  • 

A  series  of  measurements  has  been  undertaken  to  examine  the  0  + 


reaction  in  a  laser  configuration  as  shown  in  Figure  1.  Atomic  oxygen  is 
formed  in  the  source  described  above  using  a  3  cm.  I.D.  quartz  discharge 
tube  and  the  acetylene  is  injected  radially  into  the  flow  in  a  teflon  reaction 
tube  18"  long  and  with  a  §•"  central  bore.  A  laser  cavity  is  formed  by  two 
gold  mirrors  and  output  coupling  is  provided  by  a  1  ram  central  hole  in  one 
of  the  mirrors.  Slots  in  the  teflon  tube  permit  observation  of  the  1R  emis¬ 
sion  from  the  chemically  formed  excited  species.  The  system  is  pumped  with 
a  100  cfm  blower  backed  by  a  Model  1397B  Welch  forepump. 


Measurements  of  atomic  oxygen  production  rate  in  this  system  were 
essentially  in  agreement  with  observations  on  the  smaller  source.  To  avoid 


SCHEMATIC  DIAGRAM  OF  THE  EXPERIMENTAL  SET  UP: 


4  7 


< 


o  < 


u 

ce 

3 

O 

Ul 


i  / 

■f 

r  ■  ■  ; 

,  i 

I 

the  influence  of  vibrationally  excited  nitrogen  on  the  chemiluminescence 

■  ■  I  »  t  f  j 

studies,  trace  amounts  of  hydrogen  were  substituted  for  the  nitrogen  pre¬ 
viously  employed  to  enhanpe  the  dissociation  of  Og.  Figures  2  and  3  sum- 
i  marize  the  performance  with  hydrogen.  .  The  long , narrow-bore  reaction 

.  i  1  ■  i  1 

channel  limited  the  'pumping  sjieed  to  about  one- third  the  capacity  of  the 

|  :  I  I 

pump  system  and  rates  of  13  litres/sec  were  observed  at  the  location  of 
i  the  pressure  tap  in  the  r^actibn  tube.  This  corresponds  to  a  flow  veloc¬ 
ity  of  83  meters/second. 

1  '  1  1 

1  1 

The  optical  detection, system  consists  of  a'  Jarrell-Ash  one-half 

I  > 

meter  Ebert  monochromator  employing  several  detectors— a  room  temperature 

i  ■  '  ' 

1  lead  sulphide  photoconductive  detector,  a  liquid  cooled  InSb  photo- 
1  voltaic  detector  and  a  liquid  Np  cooled  Au:Ge  photoconductive  detector. 

1  ■  1  1  *  1  1  1  j 

Radiation  from  the  reaction  tube  is  focused  with  a  silicon  lens  through  a 
400  c.p>.  s.  chopper  onto  the  entrance' slit  of  the  monochromator  and  the 
detector  output  is  amplified  by  a  Keithley  Model  103  Low  Noise  Amplifier  i 
The  signal  is  synchronously  detected  using  the  chopper  as  reference  with  a 

1  J.  1  1  i  1  1  1 ' 

Princeton  Applied  Research  Model  JB-5(  Lock-in  amplifier. 

The^  spectrometer  employs  a  grating  witii  105  lines/mm  blazed  at  8.6  p 
in  first  order.  It  was  intended  to  observe  th,e  CO  fundamental  spectral 

j  j  . 

region1  in  first  and  second  order,  isolated  by  a  4  micron  long-pass  filter. 
The  measured  spectrum  and  a  60Q°C  blackbody  response  is,  shown  in  Fig.  4. 

No  well  def ined . features  attributable  to  the  CO  fundamental  are1  apparent. 
These  would  be  expected  to  fall  in  the  spectral  regioh  4.6 -4.9  microns. 

The  spectrometer  slit  width  of  1.5  mm.  should  permit,  the  resolution  of 
features  approximately  .025  V  apart.  A  similar  failure  to  observe  a 
1  structural  spectral  feature  due  1  to  the  CO  fundamental  was  experienced  when 
,  operating  in  second  order  of  the  spectrometer  and  in  first  order  at  one-: 

half  the  initial  slit  width.  It  is  concluded  that  in  addition  to  any  CO 

1 

fundamental  emission  present,  there  was  also  present  a  broadband  emission 
of  substantial  intensity  throughput  this  region.  Ihe  4  micron,  long  pass 

*  1  1  (  ' 

‘  filter  distorted  the  apparent  intensity  at  short  wavelengths  on  the  scan 
1  ,  .  !  ,  ,  1  .  1  / 

shown.  Radiation  is  present  however  out  to  a+  least  0.5  microns  and  other 

scans  have  shown  similar  emission  at  wavelengths  down  to  below  3*5  microns. 

As  dry  nitrogen  was  used  to  flush  the  spectrometer  and  all  external  optical 

i  ‘  ,  ,  f 

paths,  little  or  no  water  vapour  absorption  is  apparent  on  the  blackbody 
scan.  There  does  appear  however,  against  the  broadband  radiation  from  the 
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o  02  +  He;  PHe=l.6torr 

□  02  +  He  +  H2 ;  PHe  =  l.6torr;  PHz=O.OI  torr 

A  02+He+H2;  PHe  =  l.6torr;  PH  =0.02torr. 
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0  +  CgHg  reaction,  sharp  absorption  features  identifiable  as  due  to  water 
formed  in  the  reaction.  These  features  do  exhibit  the  resolution  estimated 
above.  The  source  of  this  excess  radiation  has  not  been  identified  for  the 
present  experiment  but  makes  extremely  difficult  any  determination  of  vibra¬ 
tional  population  distribution  based  on  observations  of  CO  fundamental  emis¬ 
sion.  It  is  also  true  that  significant  absorption  due  to  water  could 
drastically  modify  the  spectrum  of  any  laser  which  might  be  operated  based 
on  the  0  +  C^Hg  reaction. 

A  series  of  spectral  scans  were  made  in  the  CO  first  overtone  region 
in  the  vicinity  of  2.5  microns  using  the  spectrometer  in  third  order.  The 
spectral  response  was  limited  by  a  2.3  micron  long  pass  filter  and  the 
rapid  decrease  in  sensitivity  of  the  lead  sulphide  detector  at  wavelengths 
greater  than  2.8  microns.  This  decrease  will  be  seen  to  be  unfortunate 
for  the  observation  of  the  higher  vibrational  levels  but  the  operation  of 
the  monochromator  in  third  order  required  a  long  wavelength  cutoff  for  un¬ 
ambiguous  measurements. 

Figure  5  shows  a  spectral  scan  for  conditions  where  the  oxygen  disso- 
iation  was  only  about  10$,  providing  a  ratio  of  oxygen  atoms  to  fuel  mole¬ 
cules  of  2.3:1*  The  band  centers  for  the  indicated  v  ■  2  transition  are 
shown  above  the  wavelength  scale.  Emission  from  v'  =  l4  is  apparent  al¬ 
though  its  relative  intensity  is  distorted  by  the  falling  detector  response. 

A  computer  program  is  being  set  up  with  temperature  and  population  parameters 
to  establish  the  relative  distribution  of  CO*  over  the  vibrational  levels 
by  comparing  observed  and  synthetic  spectra.  At  the  present  time  this  is 
not  yet  running  but  nevertheless  a  qualitative  comparison  of  the  overtone 
emission  under  a  variety  of  conditions  is  significant. 

The  addition  of  a  trace  of  Hg  to  the  oxygen  source  increased  the  dissoc¬ 
iation  to  above  30$  and  the  O/CgHg  ratio  to  8:1.  A  twenty  fold  increase  in 
overall  signal  level  was  observed  as  shown  in  Figure  6.  There  is  a  water 
vapour  absorption  structure  in  the  vicinity  of  2.58  and  2.68  microns  but  the 
disappearance  of  vibrational  structure  compared  to  that  in  Figure  5  indicates 
a  substantial  increase  in  rotational  temperature,  certainly  in  excess  of  600°K. 
Lacking  the  synthetic  spectra  it  is  not  possible  to  establish  exactly  any 
changes  in  the  population  of  higher  vibrational  levels.  Such  changes  in 
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relative  population  do  not  appear  to  be  large  as  judged  from  a  comparison  of 
Figures  5  and  6  at  the  longer  wavelengths. 

Figure  7  shows  that  an  increase  in  He  pressure  from  1.6  Torr  to  3*1  Torr 
lowers  the  rotational  temperature  as  ejected,  with  little  effect  on  the 
vibrational  distribution.  The  overall  signal  level  has  however  decreased 
an  order  of  magnitude  compared  to  Figure  6. 

Figure  8  shows  that  an  increase  in  fuel  flow  by  a  factor  of  three  over 
the  conditions  in  Figure  7  increases  the  overall  signal  level  by  a  factor  of 
2  compared  to  Figure  6  and  increases  slightly  the  rotational  temperature. 

More  significant  is  the  increase  in  relative  emission  from  vibrational  levels 
with  v'  =  9  and  greater.  The  increase  is  substantial  for  v'  >  12  and  is 
indicative  of  the  plateau  in  the  N^,  versus  v'  dependence  mentioned  above. 
Such  a  distribution  is  exactly  that  found  for  electrically  excited  CO  lasers. 
It  is  to  be  expected  that  optical  gain  should  be  present  under  operating  con¬ 
ditions  similar  to  those  employed  in  the  experiment  of  Figure  7  if  a  further 
decrease  of  Ttrans  and  Trot  could  be  brought  about.  This  might  be  accom¬ 
plished  by  liquid  N2  cooling  of  the  reaction  chamber  or  more  directly  by 
employing  a  supersonic  expansion  of  the  reaction  products.  Calculations 
will  be  initiated  shortly  to  examine  this  latter  possibility. 
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MECHANISMS  IN  THE  CS2  +  0  REACTION 
S.  H.  Bauer,  S.  Tsuchiya,  and  N.  Nielsen 

A  longitudinally  excited  pulse  discharge  has  been  set  up  in 
a  1"  I.D.  glass  tube,  about  one  meter  in  length,  for  the  study  of  the 
CS2  +  0,  reaction  in  a  flow  system.  The  laser  output  is  controlled 
with  a  grating  with  sufficient  resolution  so  that  individual  rotation- 
vibration  lines  can  be  observed.  The  block  diagram  for  the  set  up 
is  given  in  Figure  1.  Dr.  Tsuchlya  and  Mr.  Nielsen  measured  the 
dependence  of  total  power  produced  on  the  following  parameters: 

CS2,  02  and  He  pressures  for  a  range  of  exciting  voltages.  The  half 
widths  of  the  pulses  range  from  0.3  to  4.0  p  sec.  The  current  is 
measured  by  the  voltage  drop  across  a  resistor  of  0.27-40  ohms, 
depending  on  the  current  levels,  which  range  from  5  to  500  amps. 

Plots  that  show  the  effect  of  oxygen  pressure  on  the  integrated 
power  for  various  voltages  are  shown  in  Figure  2.  Similar  graphs 
were  obtained  by  varying  the  CS2  and  He  pressures.  For  example,  the 
laser  power  increases  monotonically  with  CS2  content  for  voltages 
higher  than  about  9  Kv.  The  laser  pulses  appear  20-250  y  sec.  after 
the  discharge.  In  the  02  rich  case  the.  power  peaks  rapidly  and  then 
decreases  in  an  oscillatory  manner.  These  oscillations  have  been 
correlated  with  sound  speeds  appropriate  for  the  gas  mixture,  for  waves 
propagating  in  a  direction  perpendicular  to  the  tube  axis.  The 
distribution  of  laser  light  intensity  over  the  various  vibrational 
transitions  is  shown  in  Figure  3,  for  three  peak  voltages.  The 
intensity  was  integrated  over  a  period  of  about  10  seconds  in  these 
measurements. 

About  100  lines  have  been  observed  and  assigned  to  their  respective 
vibration- rot at ion  transitions  covering  the  range  Av(14-13)  to  (2-1). 

Of  most  interest  are  the  measurements  of  delay  in  lasing  subsequent 
to  pulse  initiation,  as  dependent  on  the  vibrational  transition  (for 
three  rotational  lines);  this  is  shown  in  Figure  4.  It  is  anticipated 
that  these  data  plus  the  relative  intensities  of  chemiluminescent 
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radiation  in  the  infrared  (non-lasing  condition) ,  which  are  now  being 
measured  as  a  function  of  the  time  after  initiation,  will  provide 
the  necessary  Information  to  permit  the  unambiguous  formulation  of 
a  mechanism  and  values  for  the  rate  constants.  At  present  we  are 
operating  under  the  assumption  that  the  high  voltage  discharges 
dissociate  both  the  -*■  20  and  CS2  CS  +  S;  and  that  the 
subsequent  reaction  CS  +  0  -*•  C0^  +  S  leads  to  the  vibrationally 
excited  CO.  However,  in  the  low  voltage  discharges  a  chain  mechanism  must 
be  present.  This  is  represented  as  follows: 


pulse 

discharge 


20 


X  BXB  w  no  x  qr>  SO  +  O  ^  ~  V  SO  +  O 

+  °  19  kcal^-CS  80  J 2  12. G  kcal  2 

S’ 

♦“SK?  ♦  8 . - 8  +  °2  i^r-'80  +? 
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Whether  the  four  center  reaction: 


CS  +  o2 


exo  . 
81  kcal 


CO  -f  SO 


plays  a  role  is  yet  to  be  determined;  this  is  one  of  the  important 

questions  we  hope  to  answer.  After  completing  recording  of  the 

chemiluminescence  in  the  infrared,  we  plan  to  measure  the  absorption 

by  CS  (in  the  UV) ,  both  under  lasing  and  non-lasing  conditions,  and 

thus  to  correlate  the  production  of  CS  with  its  destruction  during 
+ 

CO  lasing. 
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EPR  CHEMICAL  KINETIC  RATE  DETERMINATION 


The  kinetics  of  chemical  reaction  rates  may  be  studied  and  absolute 

reaction  rates  determined  by  use  of  Electron  Paramagnetic  Resonance  (EPR) 

techniques.  The  types  of  reactions  that  may  be  studied  with  this  method 

are  those  in  which  free  atoms  are  present  as  either  reactants  or  products 

or  both  since  most  free  atoms  in  their  ground  states  possess  a  magnetic 

moment  and  are  therefore  paramagnetic.  Certain  molecules  such  as  02*  NO 

are  also  paramagnetic  and  reactions  involving  them  can  also  be  studied . 

In  addition,  when  atoms  or  paramagnetic  molecules  are  present  in  excited 
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electronic  states  with  sufficient  number  density  (  n>  10  cm  )  they 
may  be  distinguished  from  their  ground  state  counterparts  since  their 
magnetic  properties  are  nearly  always  different.  Many  reactions  of 
interest  to  chemical  lasers  can  be  studied  with  EPR  techniques  and  we 
have  begun  a  program  of  absolute  reactions  rate  determination  for  a  series 
of  these  reactions  whose  rates  are  unknown  or  in  doubts 

An  EPR  spectrometer  defects  the  absorption  of  microwave  power  by 
the  paramagnetic  species  in  a  microwave  cavity.  Maximum  power  absorption 
occurs  for  fixed  microwave  frequency  when  the  relevant  energy  levels  are 
split  by  the  Zeeman  effect  such  that  the  microwave  photon  energy  equals 
the  Zeeman  energy  splitting.  If  several  paramagnetic  species  are  simul¬ 
taneously  present  in  the  microwave  cavity  they  are  separately  distinguish¬ 
able  since  each  will  display  a  absorption  resonance  at  a  different 

magnetic  field  value.  A  good  EPR  spectrometer  has  high  sensitivity  and 
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S/N  >1  can  be  observed  at  specie  densities  of  10  -10  cm  at  pressures 
of  approximately  1  Torr  and  T-300*K.  Moreover,  by  comparing  the  absorption 
strength  of  the  sample  specie  with  that  of  a  reference  specie  (frequently 
taken  to  be  02)  absolute  density  determinations  are  routinely  made. 

Chemical  reactions  are  studied  in  the  following  way  for  a  reaction 
of  the  type  A+BC  ^  AB+C.  Atoms  A  are  produced  upstream  of  the  cavity 
in  a  rapidly  flowing  gas  stream  and  the'r  absorption  strength  S  at  the 
cavity  is  determined.  The  signal  S  is  directly  proportional  to  the 
concentration  of  A,  [A].  Molecules  BC  are  rapidly  mixed  into  the  flowing 
stream  cf  A  at  different  points  upstream  of  the  cavity  and  for  each 
location  of  gas  mixing  [A]  or  S  is  determined  at  the  cavity.  The  farther 


upstream  the  injection  occurs  the  longer  the  reaction  can  proceed 
before  [A]  is  determined.  The  gas  flow  velocity  is  measured  so  that 
distance  upstream  is  equivalent  to  time  of  flight.  If  [A]  -  aS  then 
^1-  -1<2  [A]  [BC]  and  S  "-^oe  “*2*  where  k2  -  k2(BC] 

and  the  reaction  is  run  under  conditions  that  [BC]  is  much  greater  than 
[A]  and  thus  can  be  considered  essentially  constant.  A  plot  of  S  vs. 
distance  upstream  on  semi  log  paper  yields  a  straight  line  (for  the 
assumptions  stated)  whose  slope  gives  k^  and  therefore  k2  .  Under  the 
assumptions  stated,  only  relative  concentration  determinations  are 
required.  More  sophisticated  variations  upon  the  above  method  will  be 
discussed  later. 

He  have  assembled  a  high  sensitivity  EPR  spectrometer  (np  >  lO^can 
be  measured)  from  existing  parts  and  coupled  it  to  a  400  cfm  pumping 
facility.  Means  for  gas  injection  at  continuously  variable  upstream 
locations  have  been  established  and  calibration  procedures  for  gas  flow 
and  paramagnetic  specie  density  have  been  developed.  He  are  equipped 
to  handle  most  gases  including  reactive  species  such  as  NO,  MF  and  the 
halogens.  Atoms  are  produced  by  molecular  dissociation  in  a  2.45  GHz 
microwave  discharge  energized  by  a  80  watt  CH  magnetron.  The  following 
paramagnetic  species  have  been  observed  with  excellent  signal  to  noise: 

H,  N,  0,  Cl,  F,  NO,  02>  The  reactions  chosen  for  initial  study  were 

h+f2  !*2»af+f,  (1) 

H+ClJ^HCl+Cl,  (2) 

NO+F^2>  NOF+F.  (3) 

Reactions  (1)  and  (2)  are  of  course  half  of  a  chain  reaction  with 
corresponding  chain  partners 

F  +  H2  ^HF  +  H  (1*) 

Cl  +  H2  ^.HCi  +  H  (2') 

Our  Initial  studies  attempted  to  use  the  method  of  data  analysis  obtained 
above  in  which  reactions  (l' )  and  (2'  )  were  ignored  since  under  the 
assumptions  previously  discussed  the  data  analysis  is  particularly  simple. 
Linear  semi  log  plots  of  S  vs.  distance  were  indeed  obtained  under 
typical  flow  conditions,  for  reaction  (1),  for  example: 
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QF2  ■  7.32x10  cin  torr  sec"1 
QH  ■  4.72x10^0^  torr  sec"1 

C  M 

QH2  -  23.2cmJ  torr  sec"1 
v  -  61.8  m/sec 
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[ f 2 1  ■  4.85x10  molecules  cm 

The  hydrogen  flow  was  chosen  very  much  less  than  that  of  F2  to  help 

validate  the  assumptions  required  to  neglect  (1*).  The  average  of  many 

determinations  analyzed  with  this  method  are: 

-13  -1  3  -1 

Reaction  (1)  :  k?  -  4.69x10  molecules  cm  sec 

c  _i  o  -1  3  -1 

Reaction  (2)  :  k2  -  3.27x10  molecules  cm  sec 

These  rates  are  somewhat  slower  than  those  quoted  by  Albright ,  et  al 

J.  Chem.  Phys.  50,  3632  (1969)  who  find 

Reaction  (1)  :  k2  ■  3.56xl0~^  molecules'^cm^sec”^ 

Reaction  (2)  :  k2  ■  3.04x10”^  molecules" ^cm^sec”^ 

We  have  therefore  undertaken  to  refine  our  data  analysis  procedure  by 

taking  Into  account  both  reactions  (1)  and  (1*)  and  (2)  and  (2')  in  the 

following  manner.  We  can  find  operating  conditions  such  that  *-k[H] 

ie  a  linear  plot  of  S  vs.  x  on  semi  log  paper  is  obtained.  Using  (1)  and 
(l')  we  have  k  -  We  shall  determine  [A2]  and 

{F21  from  flow  measurements  and  {F]  and  [A]  from  absolute  EPR  density 
determinations.  Knowing  K  and  S  at  two  points  on  the  ssml  log  plot 
we  can  determine  both  k2  and  k^.  We  are  currently  completing  this 
series  of  measurements. 

We  are  also  partway  into  the  determinations  of  k^  for  reaction  (3). 

By  operating  at  low  [NO]  we  can  avoid  consideration  of  the  competing 
reaction  NO+F  -*• OHF  since  both  [NO]  and  [F]  can  be  made  small  compared  te 
[F2J.  While  we  hive  observed  F  production  by  reaction  (3)  and  have 
observed  paramagnetic  absorption  by  NO  we  are  currently  fabricating 
another  microwave  cavity  in  which  a  larger  electric  field  component  will 
Interact  with  the  gas  as  compared  with  our  present  cavity  that  Is 
optimum  for  magnetic  dipole  transitions.  Resonant  obsorption  by  NO  can 
be  observed  with  greater  sensitivity  using  electric  dipole  transitions 
and  thus  operation  at  smaller  [NO]  will  be  possible. 

After  completion  of  these  studies  we  would  like  to  study  the  0+CS2 
reaction  in  which  0,  s  and  SO  are  observable  paramagnetic  species.  The  rates 
for  this  reactlon(s)  are  not  known. 
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STUDY  OF  RAPID  GAS  MIXING 
Donald  L.  Turcotte 


A  mixing  wind  tunnel  has  been  built  to  study  turbulent  diffusion 
at  low  pressures  and  to  test  mixing  devices.  The  mixing  tunnel  is 
illustrated  in  Pig.  1.  The  plenum  is  divided  into  two  parts  by  a  splitter 
plate.  The  two  ga6es  that  are  to  be  mixed  are  metered  into  the  two  plenums. 
For  the  experiments  reported  here  the  two  gases  used  were  helium  and  argon. 
The  splitter  plate  ends  with  a  sharp  edge  at  the  entrance  to  the  mixing 
section.  The  mixing  section  has  a  3/8"  x  3/4"  area  and  is  4"  long.  The 
mixing  section  exhausts  into  a  plenum.  A  removable  section  at  the 
entrance  to  the  mixing  section  holds  screens  or  grids  which  induce  turbulence 
in  the  mixing  section. 

Samples  of  the  mixed  gases  are  obtained  using  a  probe.  This  probe 
can  be  translated  across  the  mixing  region  and  measurements  can  be  obtained 
at  three  axial  positions.  Gas  samples  extracted  from  the  mixing  region 
pass  through  a  thermoconductivity  cell.  This  cell  is  illustrated  in 
Pig.  2.  The  cell  contains  four  heated  thin  wire  resistance  elements.  The 
sample  gas  passes  over  two  of  the  wires.  A  reference  gas  (in  these 
experiments  argon  from  the  argon  plenum  was  used  as  the  reference  gas) 
passes  over  the  two  remaining  wires.  The  four  wires  comprise  a  bridge 
circuit  which  measures  the  difference  in  voltage  drop  between  a  wire  in  the 
reference  gas  and  a  wire  in  the  sample  gas.  This  voltage  difference  is 
related  to  the  thermal  conductivities  of  the  two  gases  by 
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In  order  to  evaluate  the  constant  C  ,  the  probe  is  placed  in  pure  helium 
outside  the  mixing  region.  With  the  sample  ga  helium  and  the  reference  gas 
argon,  the  two  thermal  conductivities  are  known,  the  AV  is  obtained 
and  C  is  evaluated.  Using  this  C  (  C  is  obtained  for  each  run) 

the  AV  is  obtained  at  several  points  within  the  mixing  zone.  The 
thermal  conductivity  as  a  function  of  the  per  cent  of  helium  in  the  helium- 


Balance 


FIGURE  2 


a  PR  on  mixture  is  given  in  Ref.  1.  The  thermoconductivity  cell  is 
sensitive  to  concentrations  of  helium  in  argon  as  low  as  1  part  in  1000. 

A  complete  schematic  diagram  of  the  apparatus  is  shown  in  Fig.  3. 

In  order  to  calibrate  the  mixing  tunnel  the  laminar  mixing  of  helium 
and  argon  has  been  studied.  Assuming  a  uniform  velocity  U.  ,  the 
molal  concentration  of  helium  Cnm  satisfies  the  equation 
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where  x  1  is  measured  along  the  channel  and  y  is  measured  across  it. 
Dfofi  is  the  binary  diffusion  coefficient  for  the  mixing  of  helium  and 
argon.  Since  mixing  occurs  across  only  a  fraction  of  the  channel,  the  ap¬ 
propriate  boundary  conditions  for  Eq.  (l)  are  C^e  0  as  y  — *  — 
and  Cna  — >  /  as  y  — >  +  oo  .  The  solution  of  Eq.  (l) 

that  satisfies  these  boundary  conditions  is 


CH.  =  0.f(l  +  erf  {tIc£x]  *}) 


Measurements  of  the  laminar  mixing  profile  for  helium  and  argon  have  been 
carried  out  for  a  range  of  velocities  and  pressures  with  X  =  0.87". 

The  results  are  given  in  Fig.  4.  Good  agreement  between  theory  and  experiment 
la  obtained  if  the  binary  diffusion  coefficient  Dh*a  (reduced  to 

p 

standard  conditions,  300°K  and  1  atm)  is  taken  to  be  0.18  cm  /sec.  This 
is  a  factor  of  four  less  than  the  value  quoted  in  the  literature  (Ref.  1). 

Reasons  for  this  discrepancy  are  being  investigated,' 

‘  ■'  I 

In  order  to  study  the  influence  of  turbulence  induced  by  screens  on  the 
mixing  process,  concentration  profiles  were  obtained  for  two  different 
screens  (with  Wire  diameter  of  0.012"  and  0.018").  The  results  are 

'  i  1 

given  in  Fig.  5  along  with  the  profile  for  laminar  diffusion  (no  screen). 

For  the  finer  screen  the  concentration  in  helium-rich  mixtures  was 

unaffected  but  a  turbulent  enhancement  of  mixing  is  noted  for  argon- rich  mixtures. 

,  ,  < 

This  difference  can  be  explained  by  the  different  densities  of  argon  and 
helium  and  the  resultant  differences  in  Reynolds  number  based  on  screen  wire 
diameter.  The  Reynolds  number  for  the  screen  in  argon  is  a  factor  of  ten 

!  i 

larger  thhn  the  Reynolds  number  for  the  screen  in  helium.  For  the  coarser 
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FIGURE  5 


S’? 


screen  no  additional  mixing  is  observed  in  the  argon  but  the  mixing  in  the 
helium  is  now  as  effective  as  in  the  argon.  Apparently  the  transition 
Reynolds  number  for  the  screen  in  helium  has  new  been  exceeded.  Initial 
calculations  indicate  that  screens  are  effective  in  generating  turbulence 
at  a  Reynolds  number  (based  on  screen  wire  diameter)  as  low  as  10.  An 
understanding  of  this  phenomenon  at  such  low  Reynolds  numbers  requires 
further  study. 

Initial  studies  have  been  carried  out  with  a  coarse  grid  (l/l8" 
diameter  wire)  for  generation  of  turbulence.  Results  indicate  that 
substantial  turbulence  is  induced  and  mixing  is  increased  by  up  to  an  order 
of  magnitude. 

The  results  of  a  test  of  a  typical  chemical  laser  injectGr  is 
shown  in  Fig.  6.  The  injector  is  a  1/8"  diameter  cylinder  with  0.007" 
diameter  holes  spaced  3/16"  apart.  Injection  takes  place  perpendicular  to 
the  flow.  Helium  was  injected  at  sonic  velocity  into  an  argon  flow  with 
a  speed  of  150  ft/sec  and  a  pressure  of  30  torr.  A  concentration  profile 
for  the  helium  2"  downstream  of  the  injector  is  shown  in  Fig.  6.  Also  in¬ 
cluded  is  the  profile  expected  for  purely  laminar  mixing.  It  is  seen  that 
the  agreement  between  the  measured  profile  and  the  profile  expected  for 
purely  laminar  mixing  is  quite  good.  This  indicates  that  any  turbulence 
induced  by  the  injector  has  not  been  effective  in  speeding  up  the  mixing. 
This  is  consistent  with  other  results  which  indicate  that  single  obstacles 
at  the  centerline  of  mixing  are  not  effective  in  generating  turbulence. 
Rather,  a  coarse  grid  of  turbulence  producing  elements  is  required. 

The  apparatus  is  now  available  for  measuring  the  mixing  from  any 
proposed  injector  for  a  wide  range  of  pressures  and  subsonic  velocities. 
Complete  mixing  profiles  can  be  obtained  1",  2f"  and  3|"  downstream  of 
the  injector. 

References : 

1.  S.  Chapman  and  T.  G.  Cowling,  The  Mathematical  Theory  of  Non-Uniform 
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HIGH  PRESSURE  VIBRATIONAL  RELAXATION 


A.  Introduction 


A  program  of  experimental  determination  of  high  pressure  vibrational 
relaxation  rates  measured  over  a  wide  temperature  range  was  initiated  this 
past  year.  At  the  present  time,  such  data  does  not  exist  since  prior 
determinations  (in  CO^  for  example)  using  the  spectrophone  or  ultrasonic 
dispersion  measurements  were  not  accurate  and  were  taken  only  up  to  1  atm. 
and  over  a  very  limited  temperature  range.  Our  plan  is  to  measure  vibra¬ 
tional  relaxation  from  selected  molecular  levels  in  laser  molecules  (starting 
with  CO2 (00°1))  in  the  pressure  range  1-100  atm.  and  temperature  range 
77°K  to  greater  than  600°K  with  a  variety  of  molecular  collision  partners. 

The  need  for  reliable  experimental  determination  of  high  pressure  vibra¬ 
tional  relaxation  times  In  molecular  laser  molecules  stems,  of  course, 
from  the  desire  to  operate  molecular  lasers  at  high  pressures.  Already, 

CO2  lasers  operating  at  1  atm.  are  commonplace  in  the  laboratory  and  the 
operation  of  lasers  at  substantially  higher  pressures  is  being  contemplated. 
The  motivation  for  high  pressure  operation  is,  of  course,  to  Increase  the 
stored  energy  density  and  permissible  repetition  rate  in  molecular  laser 
devices . 

B.  Attempts  to  Operate  a  4.3  pm  CO 2  Laser 


The  method  we  have  chosen  to  obtain  vibrational  relaxation  data  is 
the  high  pressure  adaptation  of  the  laser  excited  fluorescence  technique 
that  has  been  successfully  applied  to  CO2,  ^0,  HF,  etc.,  at  pressures  of 
one  to  ten  Torr.  The  principal  requirement  for  such  experiments  to  be 

possible  is  to  radiatively  populate  the  molecular  level  whose  vibrational 
relaxation  is  to  be  studied  in  a  time  short  compared  to  the  relaxation 
time.  In  the  CC^OO0!)  experiment  a  10  n.  sec.  pulse  at  4.3  ym  and 
substantially  resonant  with  C02(00°0  -  00°1)  vibration-rotation  absorption 
is  required.  Our  initial  plan  was  to  obtain  the  4.3  ym  laser  from  the  00®2  - 
00*1  transition  in  C02  in  accordance  with  a  report  by  Javan  et  al1.  In  this 
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work  the  00*2  -  00°1  Inversion  was  supposedly  obtained  In  a  two  step 

process!  First,  a  pulsed  discharge  In  a  C02  laser  mixture  produced  a 

high  vibrational  temperature  In  the  (00*n)  system  of  CO^ ;  second  an  Intense 

pulse  of  00*1  -  10*0  laser  radiation  (obtained  from  a  Q  switched  CO^ 

laser)  was  coupled  through  the  medium  to  saturate  the  00*1  -  10*0 

transition  thereby  inverting  00*2  relative  to  00*1.  Our  initial  plan 

was  to  enhance  operation  of  this  4.3  ym  laser  by  using  a  TEA  laser  pulse 

of  higher  peak  power  to  saturate  the  00*1  -  10*0  transition  in  the  medium 

to  be  inverted  for  4.3  ym  operation.  Accordingly,  we  built  a  TEA,  three 

electrode  structure  laser  that  provided  1-2  y  sec.,  1/2  joule  pulses  of 

00*1  -  10*0  radiation.  With  a  1  meter  gain  cell  and  suitable  mirrors  for 

4.3  ym  we  were  unable  to  detect  any  4.3  y  lasing  after  a  careful  search 

over  a  broad  range  of  operating  conditions.  That  sufficient  00*1  -  10*0 

saturating  power  was  available  was  confirmed  by  lasing  the  gain  tube  at 

10.6  ym  and  observing  the  termination  of  lasing  upon  coupling  in  the  TEA 

2 

pulse.  Subsequent  communication  with  Professor  G.  Flynn  revealed  that  the 

4.3  ym  oscillation  was  extremely  weak,  had  been  observed  in  a  3  m.  gain  tube, 

and  was  probably  a  weak  transition  in  a  combination  band  10*2  -  10*1.  This 

3 

tentative  assignment  is  discussal  briefly  in  a  later  publication  .  It  is 
appropriate  to  mention  that  the  4.3  ym  CO^  source  would  not  have  given 
10  n.  sec.  pulses  and  that  we  have  designed  and  assembled  all  components 
of  a  GaAs  electroptic  shutter  designed  to  pass  an  appropriate  portion  of  the 

4.3  ym  pulse. 

C.  Operation  of  a  4.3  ym  HBr  Laser 

In  order  to  secure  the  required  4.3  ym  laser  source  our  attention  turned 

to  the  HBr  chemical  laser  for  which  4.3  ym  laser  output  hud  been  observed 
4 

by  Deutsch  .  Deutsch's  observations  indicated  that  several  coincidences  between 
HBr  laser  lines  and  C02  could  be  expected  to  within  0.1  cm  *,  and  thus  pressure 
broadening  vould  not  have  to  be  relied  upon  to  achieve  appreciable  absorption 
in  C02(00*0  -  00*1).  A  1  meter  long,  1"  ID,  longitudinally  pulsed,  low 
pressure  discharge  in  Hj-r  Br2  was  made  to  lase  on  chemically  produced,  excited 
HBr.  Lasing  was  achieved  in  both  the  free  running  (2  Au  mirrors)  and  grating 
controlled  modes.  Free  running,  peak  powers  of  10  watts  inside  the  cavity 
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were  observed  at  a  total  operating  pressure  of  about  5  torr.  Pulse  length 
vas  about  5  y  sec.  The  absorption  of  this  radiation  by  CO^  and  NjO  was 
checked.  Of  the  17  lines  made  to  lase  in  this  system,  4  were  abosrbed  so 
heavily  in  NjO  that  an  a  of  1/meter  was  observed  at  pressures  of  less 
than  5  torr.  Later  studies  with  the  TEA  HBr  lanar  to  be  described  revealed 
2  or  more  lines  in  CO2  with  correspondingly  good  coincidences.  The  other 
lines  required  pressures  of  100  to  500  torr  to  give  a  -  1/meter. 

D.  The  High  Pressure,  Transverse  Discharge  HBr  Laser 

1.  Experiments  with  a  TEA,  HBr  Laser 

The  rather  modest  output  power  of  the  longitudinally  pulsed  HBr  laser 
was  not  sufficient  for  purposes  of  the  fluorescence  experiment.  Consequently, 
efforts  were  made  to  raise  the  output.  Experiments  at  much  higher  pressures, 

10  to  400  torr,  resulted  in  factors  of  100  improvement  over  the  old  design. 

The  laser  consisted  of  1  meter  long  TEA  pin  laser  using  1000  ohm  resistors 
fed  from  a  0.1  yf  capacator  charged  to  20  KV.  This  multi-pin  geometry  had 
to  be  reverted  to  because  the  gases  used  gave  exceedingly  poor  performance 
in  the  3-electrode  pre-ionlzing  structure.  The  discharge  contracted  into  one 
bright  arc  for  any  pressures  in  excess  of  about  40  torr. 

A  mixture  of  Hj  and  Brj  in  this  multi-pin  laser  gave  powers  on  the  order 
of  1  kw  inside  the  cavity  at  pressures  of  30  torr.  The  pulse  length  depended 
on  the  pressure,  but  was  about  2  y  sec.  under  conditions  of  highest  power. 

Both  grating  and  2  mirror  operation  was  obtained  (the  powers  quoted  in  both 
the  TEA  and  the  long  tube  experiments  were  measured  with  a  2  mirror  cavity) . 

It  should  be  noted  that  more  lines  could  be  made  to  lase  in  the  TEA  than 
with  the  longitudinal  discharge  tube. 

A  series  of  quantitative  tests  were  made  to  determine  the  effects  of 
a  variety  of  chemical  species  on  the  performance.  These  included  the  usual 
gases  (Nj,  He,  02,  etc.)  as  well  as  a  series  of  bromlnated  hydrocarbons  in 
analogy  to  the  use  of  hydrocarbons  in  the  HF  laser  as  sources  of  the  F  or  H 
atoms,  or  both.  In  total,  19  species  in  addition  to  Hj  end  were  introduced. 

The  testing  scheme  used  on  any  given  substance  consisted  of  establishing 
lasing  on  a  good  coincidence  line  of  HBr  with  NjO  in  a  Hj/B^  mixture  at 
low  pressure  (about  5  torr)  and  then  progressively  introducing  more  of  the 
substance  until  the  power  output  changed  appreciably.  Up  to  approximately 
5  to  10  times  this  amount  was  then  introduced  to  see  if  the  trend  was  monotonic. 
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and  Br  +  RH  -*•  RBr  +  H  -  (9) 

or  Br  +  RH  -*■  HBr  +  R  -  (10) 

As  we  are  Interested  in  short  times*  we  have  only  shown  above  reactions 
involving  reactants  not  products.  The  reverse  of  reactions  (1)  ...  (10) 
above  will  be  present  but  be  proportional  to  [intermediate]  [produce]  and 
thus  be  slow.  Reactions  Involving  reactants  and  products  will  be  one 
of  two  forms. 

AB  +  CD  -*■  ABC  +  D  (11) 

or 

AB  +  CD  -*■  AC  +  BD  (12) 

Reaction  (11)*  if  ABC  exists*  will  contribute  to  the  number  of  free  atoms 
present  and  will*  in  general*  help  matters.  (12)  is  a  our  center  exchange 
and  can  be  expected  to  have  a  high  activation  energy  and  consequently  to  be 
slow  at  room  temperature.  For  the  same  reason*  reactions  between  the  rearjb^nts 
e«g.  * 

RH  -r  Br  X  -V  RBr  +  HX  (13) 

RH  +  Br  X  R  X  +  HBr  (14) 

will  not  be  considered.  This  contention  is  supported  by  experimental 
evidence  against  H2  +  Br2  -*•  2HBr5. 

Deactivation  of  the  HBr*  produced  in  (1>  by  any  of  the  other  species 
present*  particularly  the  Br  X  and  the  RH,  will  limit  lasing,  but  is  hard 
to  estimate. 

In  general,  we  are  looking  for  reactions  that  are  easy  to  initiate 
with  both  chain  propagating  steps  that  are  fast;  with  competing  reactions 
that  are  slow*  or  if  fast,  feedback  into  the  chain  by  (7)  and  (9),  not  by 
(8)  and  (10).  If  any  of  these  steps  is  unusually  exothermic*  the  possibility 
of  V-V  pumping  from  excited  reactants  is  a  hopeful  possibility. 

An  analysis  of  the  thermochemistry  of  reactions  with  various  additions 
was  made  with  the  following  results: 

a.  Reactions  with  hydrocarbons  or  brominated  hydrocarbons  result  in 
highly  endothermic  chain  propagation  steps  or  less  exothermic  not 
than  the  direct  reaction  H  +  Br2  -*■  HBr  +  Br.  Thus  hydrocarbons 
will  not  help  in  the  performance  of  the  HBr  laser*  either  as 
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This  was  repeated  starting  with  a  high  pressure  (about  100  torr)  of 
the  H2/Br2  nix.  If  the  substance  contained  Br,  the  opposite  was  done. 

For  those  substances  which  increased  the  output ,  various  mixtures  of 
these  and  H2  and  Br2  were  tried  to  optimize  the  peak  power. 

The  results  are  summarized  below.  The  symbols  +f,  +,  0,  — ,  and 

t  indicate  a  strong  increase,  a  mild  increase,  no  appreciable  effect,  a 
mild  decrease,  a  strong  decrease  and  termination  of  the  lasing. 


Air:  + 

He :+  (High-press);  -  (low  press) 
Ar:  +f 


COjlt 

h^:-f  (high-press) ;  -(low  press) 

V:C 

0:+(high  press);  +f  (low  press) 
HBrt  t 
HI:— 

2*  Chemistry  of  the  HBr  Laser 


Chain  prop  steps: 

H  +  Br  X  ->  HBr*  +  X 

— - 

(1) 

X  +  RH  ->  RX  +  H 

— 

(2) 

Initiation 

RH  +  m  -»■  R  +  H  +  m 

— 

(3) 

BrX  +  m-*Br  +  X  +  m 

— 

(4) 

Competing 

H  +  Br  X  ->  HX  +  Br 

— 

(5) 

X  +  RH  ->  HX  +  R 

— 

(6) 

R  +  Br  enter  chain  by  R  +  Br  X  ->  RBr  +  X  —  (7) 

R  +  Br  X  ->  RX  +  Br  -  (8) 


CH4:- 


ch3 

Br  :0 

ch2 

Br2:- 

CHj 

BrCl: 

CH 

Br3:  — 

C2 

v- 

C2 

H5  Br: 

or 


sources  of  H  or  Br  atoms. 

b.  Use  of  the  non-equilibrium  abundance  of  atomic  hydrogen  present  in 
some  other  chain  reaction  to  feed  into  the  hot  H  +  Br2  reaction 

is  suggested.  The  addition  of  02  to  the  discharge  should  be 
beneficial  if  the  H2  +  02  reaction  can  be  used  as  a  source  of 
atomic  hydrogen. 

c.  The  brute  force  technique  of  using  some  hydride  with  very  weak  bonds 
was  investigated  as  a  useful  additive.  As  an  example,  Silane 
(SiH^)  decomposes  immediately  in  an  electric  discharge  or  at 

T  >  400°C  and  can  be  considered  as  a  source  of  atomic  hydrogen 
for  the  H  +  Br2  ■+■  HBr  +  Br  reaction.  In  addition  since  the 
Si-X  (X  is  a  halogen)  bonds  are  much  stronger  than  the  Si-H 
bonds.  Silane  might  be  effective  as  a  hydrogen  atom  source  via 
the  chain  reaction 

H  +  Br2  -*>  HBr  +  Br  ,  AH0  *  -  41.4  Real 

Br  +  SiH^  +  SiH3Br  +  H  ,  AH®  *»  +  0.4  +  5  Real 

in  which  the  second  reaction  is  much  more  favorable  than  the 
corresponding  reaction 

Br  +  H2  HBr  +  H  ,  AH°  ■  +  16.7  Real 

3.  Correlation  with  Experimental  Results 

The  addition  of  any  of  the  hydrocarbons  or  brominated  hydrocarbons  tried 
did  not  increase  the  laser  power  or  support  lasing  by  serving  either  as  a 
source  of  hydrogen  atoms  or  of  bromine  atoms.  The  enrgetics  are  unfavorable. 
At  best,  they  deactivate  the  HBr. 

The  addition  of  oxygen  increased  the  peak  power  by  a  factor  of  5  at 
low  pressures  and  somewhat  less  as  the  pressure  of  the  ^  +  Br2  was  increased. 
The  pressure  dependence  is  not  understood  at  present,  but  the  increase  is 
most  certainly  due  to  a  chemical  reaction  rather  than  intermolecular  V-V 
pumping  as  irt  the  case  of  N2  -  C02>  The  addition  of  N2  increased  the  peak 
power  by  at  most  a  factor  of  2  at  low  H,,/Br2  pressures.  At  higher  pressures 
it  decreased  the  peak  power.  N_  vibrations  are  nearly  resonant  with  HBr 
vibrations  (Av  <,  ~  100  cm  )  and  thus  the  increase  at  low  pressures  could  be 

explained  in  terms  of  resonant  pumping  as  in  the  case  of  C07.  However,  the 

-1  /‘  -1 
vibration  frequency  of  02;  ~  1550  cm  .  The  mismatch  with  HBr  is*  ~  800  cm  . 


t  I 
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Considering  the  size  of  kT  this  Is  hardly  resonant.  Since  the  masses  and 
the  Leonard- Jones  parameters  are  roughly  comparable  for  the  two  molecules, 

i  I  i 

It  Is  hard  to  attribute  the  Increased  power  to  anything  but  the  chemistry 
of  the  mixture. 

The  effects  of  the  other  gases  added  to  the  H^/B^  mixture. hold  few  , 
surprises.  The  addition  of  C02  and  ^0  terminated  lasing  due  to  their 
strong  absorption  at  the  frequency  to  which  the  laser  was  tuned  (with  the 

l  '  ' 

grating).  He  and  Ar  Increased  power.  This  was  concurrent  with  the  Increased 
diffuseness  of  the  discharge  upon  the  addition  of  these  gases  and  was  ' 

'  '  I 

probably  due  to  the  usual  Increased  electron  mobility  In  these  gases  (vs. 

electronegative  gases).  Cold  HBr  terminated  lasing  due  to  absorption  or 

‘  ’  1  / 

resonant  transfer.  Absorption  Is  less  likely  since  the  ,laslng  transition 
was  either  4-3  or  a  3-2.  HBr  Introduced  externally  could  not  be  made  to 
lase.  The  exothermicity  was  needed.  Monobromoethape  terminated  lasing, 
whereas  ordinary  etnane  did  not.  This  remains  a  mystery.  /  , 

i  i  < 

He  are , currently  designing  a  larger  TEA,  HBr  laser  and  seating  up 
to  experimentally . study  the  effect  of  Silane  addition  to  laser  performance. 

With  the  addition  of  the  beneficial  additives  discussed  above  the  performance 
of  our  tlBr  laser  Is  entirely  sufficient  to  proceed  with  the  vibrational 

i  .  !  1  I  1 

relaxation  experiments  and  further  beneficial  additives  will  simply  Improve 

l  • 

an  already  satisfactory  signal  to  noise  ratip.  He  performed  a  crude 
vibrational  relaxation  experiment  using  C02  in  a  cell  at  1  atm.  and  the 
simple,  1  meter,  TEA  pin  HBr  laser  without;  additives  described  above. 

I  • 

A  lmV  fluorescent  signal  Has  observed  with  100  watts  of  HBr  laser  excitation. 

i  1  1 '  1 

Since  we  conservatively  expect  1  Kwatt  HBr  output  with  additives  we  calculate 
that  at  100  atm.  C07  pressure  we  will  have  10  mV  of  signal  after  signal 
averaging.  This  is  a  minimum  figure  also,  sinee  significant  improvements 

'  1  i  /  ! 

can,  be  expected  in  fluorescent  collector  efficiency  with  proper  cell  design. 

4 .  Other  Aspects  of  Experimental  Design  > 

‘He  are  currently  completing  drawings  for  the  construction  of  the 

1  .  i 

optical  cell  and  the  high  pressure  generating  equipment.  Vacuum  and  gas 

handling  facilities  have  been  Installed.  Optical  and  signal  processing 

equipment  are  on  hand.  He  expect  final  experimentation  to  begin  within 

60  days.  In  the  interim  time  we  shall  optimize  the  final  HBr  laser,  study 

1  1  ,  ,  /  < 

its  output  spectrum  and  determine  the  degree  of  coincidence  in  all  favorable 

cases  with  CO2  and  N2<>.  ,  ,  ( 

1  1  . 
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Uiser-iiiduccd  vibrational  fluorescence  experiments  have  been  accomplished  with  a  newly  developed  cw  chemical 
laser.  The  design  and  operating  characteristics  of  this  small-scale  transverse  flow  chemical  laser  source  are  pre¬ 
sented.  The  fluorescence  measurement  technique  and  apparatus  provide  several  advantages  over  the  use  of  pulsed 
chemical  lasers  for  similar  measurements.  The  present  laser  is  also  useful  as  a  selectively  tunable  probe  for  single 
pass  gain  measurements  in  chemical  lasers. 


INTRODUCTION 

EASUREMENTS  of  the  fluorescence  from  molecular 
vibrational  bands  selectively  excited  by  the  absorp¬ 
tion  of  a  modulated  laser  beam  have  provided  a  unique 
body  of  knowledge  concerning  vibrational  energy  transfer 
processes  in  gases.1-4  At  present,  pending  the  development 
of  continuously  tunable  laser  sources  of  sufficiently  narrow 
linewidth  in  the  near  infrared,  fluorescence  measurements 
are  limited  to  vibrational  energy  transfer  processes  in¬ 
volving  molecules  possessing  absorption  bauds  matching 
existing  laser  wavelengths. 

The  exceptionally  rapid  vibrational  energy  transfer 
processes  associated  with  the  hydrogen  and  deuterium 
halides  have  heretofore  been  studied  by  monitoring  the 
fluorescence  decay  following  excitation  by  a  pulse  of 
radiation  from  HC1  or  HE  chemical  lasers.4  6  Experimental 
difficulties  associated  with  the  use  of  such  pulsed  chemical 
lasers  include  a  long  pulse  duration,  an  inherently  low 
pulse  repetition  rate  (only  single  pulses  have  been  used  to 
date),  a  large  electrical  noise  associated  with  flash  lamps  or 
electrical  discharges  used  to  initiate  laser  action,  and 
variability  in  pulse  characteristics  from  shut  to  shot. 

The  recent  development  of  small  scale  cw  chemical 
lasers  for  laboratory  research  has  made  available  con¬ 
venient  sources  of  laser  radiation  from  HE,  DF,  HC1,  and 


Fig.  1.  Apparatus  for  laser  induced  fluorescence  measurements. 


CO4-*  (presumably  Hlfr,  DBr,  DC1,  and  OH  can  be 
added  to  this  list  with  some  further  investigation10).  Such 
lasers  typically  provide  multiline  power  outputs  ranging 
from  0.1  to  2.0  W  for  conv  niently  small  scale  devices 
with  moderate  vacuum  pumping  requirements  (25-150 
liters/sec  at.  pressures  ~2-20  Torr).  Selective  single 
frequency  operation  at  somewhat  lower  power  output  is 
possible  with  use  of  a  dispersive  resonator.* 

We  wish  to  discuss  here  the  methods  and  apparatus  we 
have  recently  developed  for  successful  measurements  of 
vibrational  energy  transfer  rates  of  importance  in  the 
HF,  DF,  HF-COj,  and  DF-COj  chemical  laser  systems.11 
The  use  and  construction  of  a  small  cw  transverse  flow 
chemical  laser  are  describe.'  which  provides  an  alternative, 
superior  in  many  respects,  to  the  use  of  pulsed  chemical 
lasers  for  such  studies. 

I.  EXPERIMENTAL  APPARATUS 

Figure  1  shows  a  schematic  diagram  of  the  cw  chemical 
laser,  absorption  cell,  and  associated  apparatus  used  for 
the  present  measurements.  The  laser  operated  by  means  of 
the  chemical  reactions 

F+H.— »HF(r=«)+H,  (1) 

where 

»  =  0, 1,  2,  or  3  for  HF  and  m=*0,  1,  2,  3,  or4  for  DF. 

An  rf  discharge  was  employ  ed  to  produce  requisite  con¬ 
centrations  of  F  atoms  by  dissociation  of  SF»  by  electron 
impact.  Design  features  of  the  laser  test  section  included  a 
subsonic  flow  and  an  optical  axis  transverse  to  the  flow 
direction. 

About  90%  of  the  useful  laser  output  was  coupled  from 
the  optical  cavity  by  the  partially  transmitting  mirror 
Mi  and  entered  the  absorption  cell  after  being  focused, 
modulated,  and  then  recollimated  by  use  of  lenses  Lj  and 
I.j  and  the  mechanical  chopper.  The  additional  10%  of 
the  laser  power  output  was  coupled  from  the  cavity  by 
partially  transmitting  mirror  to  provide  a  determina¬ 
tion  of  the  laser  spectral  output  characteristics  with  a 
McPherson  model  216  monochromator-spectrograph.  All 
external  light  paths  were  purged  with  dry  nitrogen  to 
prevent  absorption  by  atmospheric  gases.  A  helium  purge 
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|j<i.  2.  l‘lioio"r.i|ili  of  civ  chemical  laser, 
(aj  U|ilical  avis;  (h)— laser  lesl  section; 
(c> — rf  discharge  luhe;  (il) -wlotvnslrcnm 
cltclr.ele;  (<•)-•  exhaust;  (f)  —purged  lij>ht 
housing  between  mirror  M i  and  the  Iwx 
(g)  which  contained  the  chopper  and  lenses 
l.j  and  i.j;  (h) — location  of  internally 
mounted  mirror  Nl(;  (it  -  H a  (or  Dj)  injec¬ 
tion  manifold;  (j)  -Nat.'l  Brewster  window; 
(k)— inlet  for  helium  pur^c ;  (I)- -pressure 
lap;  (m)-  -5  cm  i  d.  plenum  chamlicr. 


was  used  in  ihe  laser  test  section  to  prevent  absorption  by 
ground  stale  HF  or  DF  in  the  space  separating  the  inner 
optical  surfaces  from  the  main  flow. 

The  desired  lluorcscence  signal  was  provided  by  an 
InSb  photovoltaic  detector  which  was  exposed  to  radiation 
sampled  from  the  absorption  cell  through  a  narrow  hand 
interference  filter  by  means  of  lens  Lp 

The  requisite  gas  mixtures  for  the  absorption  cell  were 
carefully  prepared  and  purified  before  being  admitted  to 
the  cell.  A  continuous  llow  through  the  cell  was  provided 
to  maintain  sample  purity  despite  inevitable  surface 
reaction  and  contamination.  Absorption  cell  pressures 
were  measured  with  a  capacitance  manometer.  Cell  tem¬ 
peratures  were  monitored  with  a  Chromcl-  Aluniel  thermo¬ 
couple  with  an  insulated  stainless  steel  outer  sheath. 

The  essential  components  of  Fig.  I  arc  discussed  in  more 
detail  in  the  remainder  of  this  section. 

A.  Continuous  Wave  Chemical  Laser 

Figure  2  shows  a  photograph  of  the  cw  chemical  laser 
employed  for  these  studies.  Identifiable  in  the  photograph 
is  the  rf  discharge  tube  ~sed  to  accomplish  initial  F  atom 
production  by  partial  dissociation  of  SF»  in  a  premixed 
flow  of  SF»  and  helium  (see  Fig.  1).  The  capacitivelyr 
coupled  discharge  was  maintained  by  means  of  a  pair  of 
aluminum  foil  strips  of  5  cm  width  with  a  5  cm  space 
separating  them,  tightly  wrapped  around  a  16  nun  i.d. 
quartz  tube.  A  one-half  wavelength  coax’al  cable  con¬ 
nected  the  discharge  electrodes  to  the  tank  circuit  of  a 
2500  W,  27.12  MHz  rf  oscillator  amplifier  of  our  own 
construction.  Under  normal  operating  conditions  the  sup¬ 


ply  input  power  was  maintained  at  !4(X)  W,  of  which 
perhaps  about  GOO  \Y  was  delivered  to  the  gas  mixture. 
Forced  air  cooling  of  the  discharge  tube  was  provided  from 
a  fan  which  fed  a  Pyrex  duct  that  enveloped  the  discharge 
tube. 

Upon  leaving  the  discharge  and  flowing  20  cm  from  Ihe 
downstream  electrode,  the  primary  flow  expanded  into  a 
plenum  chamber  of  5  cm  inner  diameter  before  entering 
the  upstream  end  of  the  transverse  flow  chemical  laser. 
The  laser  test  section  (sec  Fig.  5)  was  machined  from  a 
single  piece  of  aluminum.  A  sonic  throat  of  5X4.5  cm 
cross  section  terminated  a  5  cm  length  of  constant  area 
channel  of  5X4.5  cm  cross  section  to  provide  subsonic 
llow  within  the  laser  test  section.  The  piimary  llow  entered 
the  test  section  through  an  inlet  of  6.4X5  cm  cross  section 
and  was  rapidly, mixed  with  a  secondary  flow  of  Hj(or  Dj) 
injected  by  means  of  two  opposing  arvavs  of  supersonic 
jets.  The  two  rows  of  jets  were  aligned  transverse  to  the 
primary  flow  with  each  row  formed  of  24  orifices  of  0.(155 
nun  diameter  spaced  2  mm  apart.  The  orifice  rows  were 
located  in  the  top  and  bottom  test  section  walls  such  that 
the  jets  intersected  the  primary  flow  at  right  angles  with 
the  lop  and  bottom  jet  axes  forming  a  staggered  array 
(sec  Fig.  3). 

The  transverse  optical  axis  was  located  6  mm  down¬ 
stream  from  the  injection  orifices  to  allow  for  initia’  mixing 
and  chemical  reaction.  The  original  design  had  provision 
to  position  the  optical  axis  at  various  locations  further 
downstream;  however,  Fig.  3  shows  a  fixed  optical  axis 
located  to  give  maximum  power  output.  Not  shown  in 
Fig.  3  arc  the  H<  injection  manifold  covers  visible  in  the 
photograph  of  Fig.  2. 
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It  was  found  necessary  to  provide  a  small  purge  How 
of  helium  to  eliminate  recirculation  of  ground  state  HF 
(or  DF)  molecules  within  the  region  separating  the  main 
flow  from  the  inner  surfaces  of  the  Brewster  window  and 
mirror  Mi  (sec  Fig.  A).  Without  the  purge  (low  the  power 
output  of  I  — >  U  hand  transitions  was  markedly  reduced 
hy  absorption  from  the  cold  HF  (or  DF). 

The  partial  flow  rates  found  to  provide  a  maximum 
power  output  for  HF  laser  operation  were:  primary  He, 
25m.l//sec ;  He  purge  flow,  0.4 m A/ /sec ;  SF«,  0,8 7m. 4// 
sec  ;  Hj,  2.1mA/  see.  'Hie  plenum  pressure  was  10.5  Torr; 
the  pressure  at  the  optical  axis  location  was  9.8  Torr. 
The  test  section  static  temperature  and  mean  flow  velocity 
were  about  515  K  and  4.1X104  cm/sec  at  the  optical  axis. 

For  DF  laser  operation,  optimum  power  output  condi¬ 
tions  were:  primary  He,  llmAf./scc;  He  purge  flow, 
0.4mA//sec;  SF«,  0. 2fim.il/ /sec ;  D»,  1.4m M/sec.  The 
plenum  pressure  was  5.1  Torr,  the  pressure  at  the  optical 
axis  was  4.65  Torr.  The  test  section  static  temperature 
and  mean  flow  velocity  were  about  480  K  and  A.75X104 
cm/'sec  at  the  optical  axis. 

The  laser  output  spectrum  for  I1F  consisted  of  the 
/,(4),  /'(5),  and  /'(6)  2—*  1  anu  1  — >0  band  transitions. 
The  total  output  power  from  a  5%  transmitting  mirror 
Mi  was  1.5  W  for  all  lines.  The  1— >0  band  transitions 
made  up  about  50%  of  this  power.  For  the  HF  laser  case 
an  additional  approximately  1.5  W  was  lost  because  of 
mirror  absorption. 

For  the  DF  laser  case  the  laser  output  spectrum  con¬ 
sisted  of  the  /’(5),  /’( 6),  and  l’(7)  transitions  on  the 
3— » 2,  2  — *  l,  and  1—  >tl  bands.  The  total  output  power 
from  a  2%  transmitting  mirror  Mi  was  0.5  W  for  all  lines, 
of  which  about  25%  was  represented  by  1  — >  0  b&nd  transi¬ 


tions.  The  power  loss  because  of  mirror  absorption  was 
negligible. 

It  should  be  noted  that  both  the  IIF  and  DF  lasers  were 
operated  at  somewhat  higher  power  outputs  than  those 
reported  here  when  F-  was  substituted  for  SF»  as  a  source 
of  F  atoms,  to  enable  initiation  of  the  chain  reaction  be¬ 
tween  Fs  and  H;  or  Dj.  Because  of  the  greater  precautions 
necessary  in  the  handling  of  Fs  it  was,  however,  much  more 
convenient  to  use  SF». 

The  configuration  of  Fig.  3  required  a  vacuum  pump 
capable  of  a  minimum  pumping  rate  of  180  liters -'sec 'at  a 
pin  nr  del  pressure  ol  about  2  Torr.  Operation  at  smaller 
flow  i...cs  and  correspondingly  smaller  power  outputs 
could  be  accomplished  by  reducing  the  sonic  throat  height 
and/or  transverse  flow  dimension  from  the  values  of  Fig.  3. 

B.  Optical  Cavity 

The  cavity  mirrors  of  Fig.  1  were  of  3  mm  thick  ger¬ 
manium  with  dielectric  film  coatings  supplied  by  Laser 
Optics  Incorporated  to  our  specifications.  The  external 
mirror  M«  had  a  4  m  radius  and  a  transmission  of  0.5%. 
A  sodium  chloride  Brewster  window  sealed  with  a  Viton 
O-ring  (see  Fig.  2)  was  necessary  to  enable  niirror  Mj  to 
be  positioned  externally.  Any  of  several  mirrors  (Mi) 
with  transmissions  of  2,  5,  and  10%  and  of  10  in  radit.s 
could  be  mounted  internally  (sec  Fig.  3).  Mirror  absorp¬ 
tion  could  be  neglected  for  the  DF  case,  but  the  II F 
mirrors  had  absorption  values  approximately  equal  to 
their  transmissions  because  of  the  presence  of  II30  in  the 
dielectric  coating.  A  mirror  separation  of  95  cm  was  used 
to  insure  the  presence  of  a  longitudinal  cavity  mode  near 
the  peak  of  the  unsaturated  gain  profile  for  a  given  laser 
transition. 
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The  flow  channel  height  of  5  mm  was  selected  to  insure 
a  negligible  loss  from  the  TEMm  mode.  In  practice  the 
laser  tended  to  oscillate  in  cither  the  TEMu  or  TEMu 
imxles  (x  axis  along  flow  direction)  unless  the  beam  was 
apertured  slightly,  which  suggests  a  somewhat  reduced 
flow  channel  height  would  be  passible. 

C.  External  Optical  System 

Lenses  L2  and  La  had  focal  lengths  of  5  cm ;  lens  Li  had 
a  2.5  cm  focus.  Lens  Li  was  made  of  Irtran  II;  lenses  L« 
and  Li  were  of  calcium  fluoride.  All  external  optical  paths 
were  purged  with  dry  nitrogen  to  prevent  atmospheric 
absorption. 

The  mechanical  chopper  was  driven  by  a  Globe  syn¬ 
chronous  motor  model  75AKXJ3-2  at  25  000  rpm  with  an 
audio  oscillator  and  amplifier.  The  chopper  wheel  consisted 
of  six  slits  of  0.25  mm  width  and  5  nun  length  aligned 
radially  and  equally  spaced  around  the  periphery  of  a 
7  cm  diam  aluminum  disk  of  2  mm  thickness.  The  slit 
edges  were  machined  to  a' 0.3  mm  thickness  by  milling 
away  material  in  the  vicinity  of  the  slits  on  the  surface 
of  the  wheel  Lacing  away  from  the  incident  beam.  The 
rccollimatcd  light  pulses  produced  with  this  chopper  and 
lens  arrangement  were  approximately  Gaussian  with  an 
FWIIM  of  4  /usee,  with  peak  powers  equal  to  78%  of  the 
cw  laser  output  from  mirror  Mi.  The  standard  deviation 
in  pulse  height  was  less  than  5%.  Narrower  pulses  would 
be  difficult  to  achieve  with  the  present  arrangement  be¬ 
cause  of  the  observed  inability  to  focus  the  3  nun  diam 
lacer  beam  to  a  spot  much  smaller  than  (he  0.25  mm  slit- 
widths  employed.  The  approximately  400/rsec  duration 
between  pulses  provided  a  convenient  time  interval  during 
which  to  observe  the  fluorescence  decay  under  appropriate 
absorption  cell  gas  compositions  and  pressures. 

A  rotating  mirror  Q-switch  arrangement  was  initially 
attempted,  instead  of  the  present  arrangement,  by  rotating 
mirror  Mj  (in  this  case  flat  and  gold  coated)  at  150-200 
cps.  The  pulse  widths  were  shorter  (about  1  jisec),  and 
the  peak  pulse  powers  were  about  2.5  times  the  cw  power 
level.  This  very  modest  increase  in  peak  Q-switchcd  power 
over  the  cw  power  level  reflects  the  small  energy  storage 
capacity  of  the  HE  system  compared  with  other  molecular 


fiG.  5.  Signal  processing  instrumentation. 


lasers  such  as  the  N*--COj  system.  The  slight  advantages 
in  peak  power  and  pulse  width  were  overcome  by  the  dis¬ 
advantages  of  piilsc-to-pulsc  fluctuations  and  relatively 
low  pulse  repetition  rate  compared  to  the  present  chopping 
wheel  arrangement. 

Though  not  usually  necessary,  it  was  advantageous 
under  conditions  of  very  low  fluorescence  intensity  to 
employ  lens  Li  to  focus  fluorescence  at  the  InSb  detector 
element.  The  filters  employed  were  manufactured  by 
Optical  Coating  Laboratory,  Inc.  For  HF  fluorescence  the 
filters  \V -02430-5  were  used;  the  filters  31*31-551  -045  and 
N -044 10-9  were  employed  for  the  DF  and  CO*  fluores¬ 
cences,  respectively.  The  HF  and  I)F  fillers  were  specially 
designed  to  pass  /{-branch  emission  on  the  1  — >0  bands, 
but  rejected  the  scattered  /’-branch  emission  from  the 
laser  pulse.  Usually  two  filters  of  a  given  type  slacked  in 
series  (at  a  slight  angular  displacement  to  prevent  multiple 
reflection  between  fihtrs)  were  necessary  to  provide  ade¬ 
quate  rejection  of  the  scattered  HF  or  1)F  /'-branch  laser 
emissions. 

D.  Absorption  Cell 

Figure  4  shows  I  lie  essential  features  of  the  absorption 
cell,  machined  from  a  single  piece  of  aluminum.  The  cell 
operated  at  total  pressures  ~  10  Torr  and  could  be  operated 
at  temperatures  from  300  to  450  K  with  a  temperature 
regulated  Nichromc  heater.  A  typical  total  gas  flow  of 
(M)2mJ//sec  was  maintained  through  the  cell  to  minimize 
the  effects  of  surface  reactions  in  causing  absorption  or 
contamination  of  the  HF  or  DF  gas  samples.  Calcium 
fluoride  windows  sealed  with  Yiton  O-rings  (see  Fig.  4) 
wiO u jwI  iuU'  nil"' at  jperiiltoft  «JR*itdfug 
over  several  months  with  no  visible  coating  or  deteriora¬ 
tion  caused  by  HF  or  DF.  Ik  cause  of  ihc  large  absorption 
coefficient  of  IIF,  the  detector  was  located  as  close  as 
possible  to  the  Brewster  window  (see  Fig.  4). 

The  gas  handling  and  purification  procedures  described 

Iwi.*  itL  i. '  Misl  BTj  {.■  (■&  -  j  -  l.iLL  'iJiuii  4 1 ii  I  ..ill  mif  a 

Ha  ii  I »  -jj  nll'I  T  tifcvt  m  rw  ti  l|£>C  Writ  iirt  nn  *  j 
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tion  by  IM).  'I  lie  total  cell  pressure  and  the  partial  pres¬ 
sure  of  each  gaseous  component  wit  bin  the  cell  were  mea¬ 
sured  with  an  MRS  model  1-1-11  l-IOSl*  Har.U ron  capaci¬ 
tance  manometer  anil  the  associated  model  170M-5  elec¬ 
tronics  unit.  Cell  temperatures  were  monitored  with  an 
Omet'.a  model  SCASS-t)2t)C.-6  Chromel-Almucl  insulated 
thermocouple  before  and  after  each  run. 

E.  Instrumentation  for  Fluorescence  Detection 

A  block  diagram  of  the  electronics  instrumentation  is 
show'n  in  big.  S.  The  lluorosccnce  signal  was  detected  by  a 
l'hilco-Ford  I  SC-30 1C  l»Sb  l’V  detector  with  a  £>*(500, 
767,  1)  of  1.45X10I0cm  (Ifz)>/\V.  The  element  was  2.25 
mm  in  diameter  and  was  cooled  to  77  K.  The  detector  was 
back-biased  to  a  aero  voltage  and  operated  with  a  10000 
load  resistor.  The  detector  RC  time  constant  with  this  load 
impedance  was  1.2/<scc. 

The  voltage  appearing  across  the  load  resistor  was  fed 
into  a  PAR  model  1 1.3  low  noise  preamplifier  which  was 
used  with  a  pass  band  from  0.03  Hz  to  3(X)  kHz  and  an 
amplification  factor  of  5000.  The  preamplifier  output 
signal  was  further  amplified  50  limes  by  the  internal  cir¬ 
cuits  of  a  Tektronix  type  532-S7  oscilloscope  and  type 
53/54D  plug-in  unit  which  had  a  bandwidth  of  2  MHz. 
The  over-all  time  response  of  the  detector',  preamplifier, 
and  oscilloscope  amplifier  was  2  usee. 

The  vertical  signal  output  of  the  oscilloscope  was  used 
to  drive  a  PAR  model  CW-l  “boxcar"  integrator.  The 
scanning  circuits  of  both  the  boxcar  integrator  and  the 
oscilloscope  were  triggered  by  a  photolransisior  syn¬ 
chronized  to  the  laser  pulse.  The  photolransisior  received 
a  pulse  of  light  from  a  light  bulb  located  on  the  opposite 
side  of  the  chopping  w.ictl  slot  through  which  the  focused 
laser  beam  passed. 

The  boxcar  integrator  output  was  recorded  on  one 
channel  of  a  dual  channel  strip  chart  recorder.  The  second 
channel  was  used  'or  display  of  a  signal  proportional  to 
the  total  cw  laser  output  provided  by  a  thin  film  bolometer 


Fio.  6.  Fluorescence  of  the  COj(00°l)  -  >  C0,(00*0)  band  followirg 
absorption  of  a  single  pulse  of  I  IF  laser  radiation  ir  a  mixture  of  Hr 
and  COi. 


Fig.  7,  Fluorescence  of  the  CO,(00°1)  —  CO, ((WO)  band  after 
signal  averaging  for  Ihe  conditions  of  Fig.  6. 


and  amplifier  which  detected  a  fraction  of  the  laser  output 
through  mirror  Mg.  This  provided  a  check  on  laser  output 
stability  while  data  were  being  recorded. 

II.  PERFORMANCE  OF  LASER-EXCITED 
FLUORESCENCE  APPARATUS 

Figures  6-8  demonstrate  typical  operation  with  the 
present  apparatus.  These  data  were  taken  at  an  absorp¬ 
tion  cell  temperature  of  295  K  and  cell  partial  pressures  of 
CO;,  0.27  Torr;  HF,  0.18  Torr;  and  argon,  9.88  Torr 
(argon  was  present  as  an  inert  buffer  to  reduce  the  rate  of 
diffusion  to  the  tube  walls).  The  fluorescence  signal  from 
either  the  HF(r*»l)  — » HF(t'*0)  band  or  from  the 
CO;((X)°l)  — ♦  C0;(00°0)  band  was  monitored  to  examine 
vibrational  energy  transfer  processes  following  the  periodic 
production  of  HF  molecules  in  the  first  vibrational  level 
by  the  induced  absorption  of  the  i>=‘.  — +0  laser  pulses.11 

Figure  6  shows  an  oscilloscope  trace  of  the  fluorescence 
signal  from  the  CO;(00°1)  C'O*(00°0)  band  excited 

with  a  single  laser  pulse.  It  is  evident  that  the  fluorescence 
signal-to-noisc  ratio  is  no  better  than  unity  for  a  single 
pulse.  This  is  markedly  inferior  to  the  performance  one 
would  expect  with  the  usual  pulsed  chemical  laser4  6  which 
is  capable  of  millijoulc  pu  ses  in  contrast  to  the  microjoulc 
pulses  achievable  with  the  present  device.  However,  the 
high  repetition  rate  and  uniformity  of  the  present  pulses 
permit  full  utilization  of  modern  signal  averaging  tech¬ 
niques  to  recover  weak  signals  front  background  noise  to 
provide  quite  good  fluorescence  traces  as  Figs.  7  and  8 
demonstrate. 

Figure  7  shows  a  chart  recording  of  the  output  of  the 
boxcar  integrator  (see  Fig.  5)  taken  under  the  same  condi¬ 
tions  under  which  the  single  pulse  of  Fig.  6  was  recorded. 
This  trace  was  taken  with  a  3  Msec  gatewidth,  a  400 //see 
time  base,  and  a  1  msec  time  constant.  The  entire  trace 
w..s  recorded  in  a  2  min  time  interval.  The  exponential 
decay  of  Fig.  7  yielded  a  value  of  the  rate  constant  for 
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Fig.  8.  The  initial  rise  of  the  fluorescence  signal  of  Fig.  7 
obtained  with  an  expanded  lime  scale. 

vibrational  deactivation  of  COj(00°l)  by  ground  state 
HF" 

In  order  to  examine  the  over-all  process  of  vibrational 
energy  transfer  between  HF(t'=l)  and  ground  slate  COt, 
it  was  necessary  to  examine  more  closely  the  rise  of 
fluorescence  signals  such  its  that  shown  in  Fig.  7."  Figure  8 
is  a  recording  taken  under  identical  conditions  with  a 
1  gsec  gate  width,  a  50  /iscc  time  base;  and  1  msec  time 
constant  in  order  to  closely  resolve  the  initial  portion  of 
the  fluorescence  signal.  Again,  the  entire  trace  was  ac¬ 
complished  in  a  2  min  time  interval. 

It  is  clear  from  these  results  that  much  weaker  signals 
could  be  meaningfully  examined  without  difficulty  by 
increasing  the  signal  averaging  time.  Our  measurements 
of  the  HF  (»=  1)  deactivation  rate  in  pure  HF  obtained  by 
monitoring  the  decay  of  HF(»«1)  —  >  HF(r=0)  band 
fluorescence"  agree  well  with  the  previous  measurement  of 
Airey  and  Fried*  who  employed  the  single-shot  pulsed 
chemical  laser  technique. 

III.  A  SELECTIVELY  TUNABLE  PROBE  LASER 

The  present  cw  chemical  laser  has  provided  a  convenient 
probe  laser  for  single-pass  gain  measurements  of  HF  and 


L»F  purely  chemical  lasers.’5  For  this  application,  the 
minor  M;  was  replaced  with  a  diffraction  grating  (Bansch 
&  Lonib  35-63-05-730),  blazed  at  2. 5 /i,  of  300  lines/nim, 
which  was  Lit  trow  mounted  so  that  the  first  order  dif¬ 
fracted  beam  sustained  laser  oscillation.  Single  frequency 
operation  on  any  of  the  individual  transitions  discussed 
in  Set.  11  was  possible  with  this  arrangement.  It  was 
important  to  insure  that  the  cavity  resonance  frequency 
for  oscillation  was  stabilized  near  the  center  frequency  of 
the  unsaturated  line  profile  to  obtain  meaningful  gain 
measurements.  For  this  purpose  the  grating  was  mounted 
upon  a  piezoelectric  ceramic  (Lansing  model  21.817) 
stressed  by  a  dc  voltage  controlled  by  a  feedback  network 
stabilized  at  maximum  laser  output  power  (Lansing 
model  80.210  lock-in  stabilizer). 
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The  optical  saturation  of  a  vibration- rotation  transition  in  a  molecular  gas  generally  cannot 
be  described  by  any  of  the  saturation  formulas  obtained  with  the  assumption  of  fixed  pump 
rates.  Its  explicit  form  is  derived  here  by  treating  rotational  relaxation  into  each  of  the  ra- 
diatively  interacting  levels  as  a  variable  pump  rate  dependent  upon  the  degree  of  saturation. 

The  formula  exhibits  a  collision-induced  power-broadening  feature  which  is  distinct  from  the 
well-known  homogeneous  power  broadening  in  form  as  well  as  in  physical  origin.  This  effect 
has  been  observed  indirectly  in  a  C02  laser. 


I.  INTRODUCTION 

The  optical  saturation  of  a  pair  of  atomic  levels 
with  fixed  pump  rates  has  been  considered  in  the 
homogeneously  and  inhomogeneously  broadened  lim¬ 
its  by  Gordon  et  al. 4  and  in  general  by  Close. 2  In 
the  analysis  of  a  vibration-rotation  transition  in  a 
molecular  gas,  however,  one  cannot  realistically 
assume  the  pump  rates  for  the  radiatively  interac¬ 
ting  rotational  levels  to  be  constant.  It  is  more 
reasonable  to  assume  fixed  pump  rates  for  the  as¬ 
sociated  vibrational  levels,  since  rotational  relaxa¬ 
tion  is  always  much  faster  than  vibrational  relaxa¬ 
tion.  Consequently,  the  saturation  formula  given 
by  Close  is  no  longer  applicable  in  the  case  of  a 
vibration-rotation  transition.  One  must  modify  the 
expression  specifically  to  reflect  the  additional 
complication  caused  by  the  variable  pump  rates. 

There  have  been  several  previous  attempts  at  ana¬ 
lyzing  the  saturation  of  a  vibration-rotation  transi¬ 
tion,3-5  but  none  of  them  succeeded  in  obtaining  an 
analytical  expression  for  the  saturated  incremen¬ 
tal  gain  which  is  valid  for  any  admixture  of  Doppler 
and  collision  broadening,  as  well  as  for  arbitrary 
field  intensities.  Furthermore,  these  theories  did 
not  examine  the  frequency  dependence  of  the  satu¬ 
rated  gain,  which,  as  we  shall  see  later,  brings 
into  evidence  an  entirely  new  kind  of  power  broad¬ 
ening. 

The  method  employed  in  arriving  at  the  saturation 
formula  consists  of  the  simultaneous  use  of  rate 
equations  and  the  Lamb  formalism  for  molecular 
interaction  with  an  electromagnetic  field*  The 
rate  equations  are  written  for  the  radiatively  inter¬ 
acting  rotational  levels  and  their  associated  vibra¬ 
tional  levels.  An  additional  relation  is  obtained  be¬ 
tween  the  populations  of  these  levels  and  the  inten¬ 
sity  of  the  saturating  field,  as  well  as  its  frequency. 
The  simultaneous  solution  of  these  equations  leads 
to  an  explicit  saturation  formula  for  a  vibration-ro¬ 
tation  transition  in  a  molecular  gas. 

For  sufficiently  low  pressures,  our  saturation  for¬ 
mula  reduces  to  that  given  by  Close,  with  the  rota¬ 
tional  relaxation  rate  representing  the  level  decay 
rate.  At  elevated  pressures,  the  half-width  of  the 
saturated  transition  increases  with  increasing  field 


intensity.  This  may  occur  even  when  the  homoge¬ 
neous  width  itself  (as  may  be  measured  from  the 
Lamb  dip)  is  not  significantly  power  broadened. 

To  distinguish  this  phenomenon  from  the  power 
broadening  of  the  homogeneous  width,  we  call  it 
collision- induced  power  broadening. 

II.  DERIVATION  OF  SATURATION  FORMULA 


In  a  model  of  rotational  relaxation  which  character¬ 
izes  the  coupling  between  one  particular  rotational 
level  and  all  the  remaining  ones  in  the  same  vibra¬ 
tional  level  by  a  single  rate  constant,  a  total  of 
four  rate  equations  is  needed  for  the  discussion  of 
the  optical  saturation  of  a  vibration-rotation  transi¬ 
tion  in  a  molecular  gas.  Denote  the  radiatively  in¬ 
teracting  rotational  states  in  the  upper  and  lower 
vibrational  levels  by  a  and  b  and  the  aggregates  of 
the  remaining  rotational  states  in  the  correspond¬ 
ing  vibrational  levels  by  A  and  B  and  their  respec¬ 
tive  population  densities  by  Na,  Nb,  NA,  and  Ng. 
Molecules  in  individual  rotational  states  can  decay 
via  rotational  relaxation  to  other  states  in  the  same 
vibrational  level,  or  by  vibrational  relaxation  (ei¬ 
ther  radiatively  or  collisionally)  into  a  different  vi¬ 
brational  level,  or  by  diffusion  to  the  wall  and  deac¬ 
tivation  thereupon.  Therefore,  so  long  as  the  mean 
free  path  of  the  molecules  is  much  smaller  than  the 
dimensions  of  the  apparatus,  we  can  write  the  fol¬ 
lowing  rate  equations: 

-yA -  Ra(Na  -  K) + DavzNa+s, 


RanI 

I 


I 

Rbnb 


r.  No 


s 


Wb 


FIG.  1.  Sche¬ 
matic  of  mo¬ 
lecular  pro¬ 
cesses  on 
which  the  rate 
equations  are 
based. 
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>'IG.  2.  Saturation  of  the  C02  gain  medium  at  two  differ¬ 
ent  frequencies  of  a  10.4-ji  band  transition. 


fj§*  =  -rA  +  ^Nb  -  Rb{Nb  -  ND  +  DbV*Nb -S, 

=  rA  -  A,  -  K) + DA*NA, 

-I*NB  -  Rb(Na  -Nl)  +  n^Nt, 
at  nb 

where  ya  and  yb  are  the  rotational  relaxation  rates, 
R„  and  Rb  are  the  vibrational  relaxation  rates,  and 
Da  =  DA  and  Db  =  Da  are  the  diffusion  constants  for 
the  upper  and  lower  levels,  respectively.  The  in¬ 
teraction  with  radiation  Is  represented  by  S  and  the 
superscript  e  denotes  equilibrium  populations  (which 
are  not  thermal  for  a  gain  medium) .  The  fractions 
l/w„  and  l/w„  are  equal  to  N*a/N\  and  N%/N%,  respec¬ 
tively. 

For  cylindrical  geometry,  the  diffusion  terms  can 
be  replaced  by  effective  relaxation  rates  with  radi¬ 
al  dependence 


These  pump  terms  will  be  shown  to  be  variables  de¬ 
pendent  upon  the  intensity,  as  well  as  the  frequency 
of  the  saturating  field.  This  dependence  Is,  In  turn, 
reflected  in  the  saturated  gain  for  the  transition 
connecting  a  and  b . 

In  order  to  make  use  of  the  results  in  Ref.  2,  we 
have  to  assume  that  a  and  b  are  nondegenerate, 
which  admittedly  is  unrealistic  for  vibration-rota¬ 
tion  transitions  because  of  the  existence  of  the  mag¬ 
netic  sublevels.  But  for  transitions  where  the  mean 
deviation  of  the  dipole  moments  of  the  degenerate 
sublevels  from  some  average  value  Is  small,  the 
approximation  should  be  a  good  one.  Such  Is  the 
case  for  the  P-  and  fl -branch  transitions,7  to  which 
we  shall  restrict  our  attention  in  what  follows. 


In  this  connection,  it  may  be  helpful  to  point  out 
that  for  linearly  polarized  light,  the  number  of  ac¬ 
tive  sublevels  (those  with  nonzero  matrix  elements) 
In  the  upper  rotational  level  Is  the  same  as  in  the 
1-vwer  whether  Ja  ~Jb  -  1  or  Ja  =Jb  + 1 .  It  Is  equal  to 
the  smaller  of  ga  =  2 Ja  +  1  and gb-2Jb  + 1,  and  we 
shall  label  It  g.  It  should  be  noted  that  Na  and  Nb 
us'd  in  the  rate  equations  above  consist  of  only  the 
active  sublevels.  That  Is,  they  are  both  #-fold  de¬ 
generate. 


With  this  understanding,  then,  we  can  use  the  basic 
results  from  Ref.  2.  The  saturated  incremental 
gain  for  an  electromagnetic  wave  traveling  in  the  z 
direction  with  Intensity  7  was  given  in  rationalized 
mks  units  as 


1  dl  Jnp1  (na  _ N^ReWix  +  iy) 
I  dz  €0^VW«  ”*»/(  1+/A>)U’ 

with 


(5) 


\iHx+iy)  =  re~ix’,y)i  /”  e~*dt< 


b’S) ■ 

Furthermore,  under  the  assumption  of  large  appa¬ 
ratus  dimension,  we  always  have  R„,  r„  «y„  and  Rb, 
r„«yb.  Defining  RA(r)  =  rA(r)+R,  and  RB(r)  =  rB(r) 
+  Rb  then,  we  simplify  the  above  equations  to  obtain 

(1) 

dNb  yb ..  c 

-^=-rA+-Wfl-S. 

(2) 

= rA  -£na  -RA(r)(NA  -  n*a), 

(3) 

=  VA  ~Ns  -  HB(r)(NB  -  Na), 

(4) 

which  are  illustrated  schematically  in  Fig.  1. 

It  Is  clear  now  that  and  y^fB/nb  are  the 

pump  terms  for  the  radiatively  interacting  levels  a 
and  b,  respectively. 


where 

u  =(2fe7'/M)1/®, 
x  =  (c/k«)(m„-  «), 
y=(cY'/uu)(l‘.!/rn)l/t, 
and 

^o=  (2f0c^//>I)y4y6y,/(ya +y5). 

Here  u  Is  the  angular  frequency  of  the  saturating 
light  beam,  w0  is  the  angular  frequency  of  the  ab¬ 
sorbing  transition  at  line  center,  y'  is  the  dephas¬ 
ing  collision  rate  [not  necessarily  equal  to  \(ya  +y»)], 
and,  according  to  the  discussion  above,  we  have 

g  i.i j.i 

where  a,  and  bt  are  the  magnetic  sublevels  of  a  and 
ft,  respectively.  In  the  limit  of  y  —  0,  R eW(x  +  iy) 
becomes  the  familiar  Gaussian  function. 

We  note  In  particular  that  NA  and  NB  In  (5)  are  varl- 
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ables.  In  order  to  write  the  saturation  formula  in 
terms  of  known  parameters,  we  need  an  additional 
relation  expressing  the  rotational  disequilibrium 
within  the  vibrational  levels.  For  that  purpose,  we 
rewrite  Eq.  (36)  in  Ref.  2  in  terms  of  our  own  sym¬ 
bols  as  follows: 


where  Na(v)  and  Nb(v)  are  population  densities  of 
molecules  with  velocity  v  and  U(v)  is  the  velocity 
distribution  function 

U(v)  =  (fiu)-1  exp (-//i/)  . 

Integrating  (6a)  over  velocity,  we  find 

..  ..  /Na  Nr\  r,  y'fic  7RelF(x+»yj| 

"•  -"•■(sr-s;)  L ‘-—r.  (T^pJ- (6b) 

Equation  (6b)  shows  precisely  how  rotational  equi¬ 
librium  is  disturbed  by  the  saturating  optical  field. 
The  simultaneous  solution  of  Eqs.  (1)  -  (4)  in  steady 
state  and  (6b)  yields 

Na_Nb^(n^_N^\  (.  fir/  I  ReW(x  +iyh~\ 
n„  ~  nb  \wa  nb  J  \  +  kuy  /„  (l  +  I/I0)[/i/ 

where 

y  =  2(l/ya  +  l/y„)-1 

and 


y  =  2(l/n„RA  +■  1/wjRb)'1. 

We  thus  arrive  at  the  following  explicit  saturation 
formula  for  vibration-rotation  transitions: 

1  rf7  fi?(K  _N|\ _  ReW(x+iy) _ 

I  dz  =  tohu  V*  “  n J  (1  +  7/70)1/2 +K(7/70)ReW'(*+  iy) ' 

where  ^ 

K = fir/ /kuy. 


We  first  investigate  the  low-pressure  behavior  of 
this  formula.  The  rates  y  and  y'  are  directly  pro¬ 
portional  to  pressure,  whereas  a  diffusion-domi¬ 
nated  y  (below  ~0. 1  Torr)  is  inversely  proportional 
to  pressure.  Therefore,  the  parameter  K  ap¬ 
proaches  zero  with  a  cubic  dependence  on  pres¬ 
sure.  However,  it  must  be  remembered  that  Eq. 
(7)  was  derived  with  the  assumption  that  the  rota¬ 
tional  relaxation  rate  is  much  faster  than  the  diffu¬ 
sion  rate  (ra,  r6«ya,  ya).  This  condition  implies 
pressures  in  excess  of  10“ 2  Torr  for  typical  mole¬ 
cules.  But  it  is  still  possible  to  have  K  « 1  for 
pressures  greater  than  this  lower  limit,  because 
the  numbers  Na  and  Nb  are  generally  quite  large. 

In  this  case  (7C«  1)  the  saturation  formula  reduces 

t0  1  dl  M.  (&M\  Rew(x+iy)  m 

I  dz  t0tiu  \»a  nb  J  (1+7/fo)1" 


Interpreting  the  rotational  relaxation  rate  in  70  as 
the  level  decay  rate,  we  see  that  Eq.  (8)  is  identi¬ 
cal  to  the  formula  given  by  Close  for  atomic  levels 
with  constant  pump  rates.  This  is  to  be  expected, 
since  in  the  pressure  regime  for  which  Eq.  (8)  is 
valid  the  levels  a  and  b  are  coupled  mainly  to  .4 
and  B,  while  the  latter  are  coupled  mainly  to  the 
region  outside  of  the  saturating  beam.  Conse¬ 
quently,  the  vibrational  populations  will  remain  in 
equilibrium  regardless  of  the  intensity  of  the  satu¬ 
rating  field.  It  is  not  surprising,  then,  that  we  ob¬ 
tain  the  result  corresponding  to  constant  pump 
rates  in  this  limit. 

For  K  on  the  order  of  unity,  a  comparison  of  Eqs. 
(7)  and  (8)  show  that  saturation  sets  in  for  a  smaller 
7/70  for  the  former  at  all  frequencies.  Further¬ 
more,  due  to  the  presence  of  the  ReW(.r  +  fj>)  term 
in  the  denominator  of  Eq.  (7),  the  transition  is 
less  readily  saturated  away  from  line  center.  This 
kind  of  frequency  dependence  gives  rise  to  the  phe¬ 
nomenon  of  collision-induced  power  broadening, 
which  we  now  discuss. 


HI.  COLLISION  -  INDUCED  POWER  BROADENING 

For  pressures  above  10-1  Torr  or  so,  y  becomes 
dominated  by  volume  quenching,  and  K  begins  to 
have  a  linear  dependence  on  pressure.  For  K  »  1, 
the  saturation  formula  becomes 

1  dl  _  fip?  fa  N|\  KelV(-r+  iy) 

I  dz  e0hu  yna  "  nb  )  l  +  (l/l0)KReW(x+iy)  ’  l  ' 

For  the  purpose  of  comparison,  we  write  down  the 
frequency  and  intensity  dependence  of  the  saturated 
linear  absorption  coefficient  for  a  collision  broad¬ 
ened  transition8: 


lrf/_  y2fyz-t-(tu-  &>n)2]' 

I  dz  1  +  ( I /!,)■/[■/■  +  (cJ  -  u>0)2]-‘  ’ 


where 


/«=€0c(y2/?V/>*). 

The  power  broadening  of  such  a  homogeneously 
broadened  line  is  well  known.  For  the  same  field 
intensity,  saturation  sets  in  more  slowly  as  the 
frequency  of  the  field  moves  away  from  line  center 
simply  because  the  small-signal  absorption  coeffi¬ 
cient  decreases  away  from  line  center.  In  other 
words,  the  saturating  field  becomes  less  effective 
in  distorting  the  level  populations  from  equilibrium 
as  its  frequency  moves  towards  the  wings  of  the 
Lorentzian  line  shape. 

Superficially,  Eq.  (9)  is  just  Eq.  (10)  with  a  Gaus- 
sian-like  function  replacing  the  Lorentzian  func¬ 
tion.  But  in  physical  origin,  the  power  broadening 
exhibited  by  Eq.  (9)  is  quite  different.  According 
to  Eq.  (9),  power  broadening  can  take  place  even 
for  a  Doppler-broadened  transition.  This  substan¬ 
tially  modifies  the  customary  notion  that  power 
broadening  is  to  be  associated  with  a  homogene- 
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ously  broadened  line.  For  vibration- rotation  tran¬ 
sitions,  power  broadening  (in  the  high-pressure 
regime)  results  from  the  varying  degrees  of  devia¬ 
tion  of  the  vibrational  level  populations  (and  hence 
the  pump  rates)  front  equilibrium  as  the  frequency 
of  a  fixed -intensity  saturating  field  changes,  for 
the  same  po*er  density,  the  pump  rates  for  the  ra- 
diatively  interacting  rotational  levels  are  altered 
the  most  for  a  saturating  wave  with  its  frequency 
equal  to  that  at  molecular  line  center.  Because  of 
the  obvious  connection  between  rotational  relaxation 
and  this  phenomenon,  we  refer  to  it  as  collision- 
induced  power  broadening. 

'  ,  ’  '  1  ' 

In  order  to  make  the  distinction  between  collision- 
induced  power  broadening  and  the  usual  potver 
broadening  of  a  homogeneous  transi  tion  absolutely 
clear,  let  us  examine  the  power  dependence  of  the 
width  of  the  hole  burned  in  the  Doppler-like  profile. 
It  can  be  shown8  that  the  hole  width  is  given  by 

yh=y'(l+I/lj1,z. 

Therefore,  the  homogeneous  component,  which  is 
Represented  by  t|ie  hole,  begins  to  broaden  for  I 
comparable  to  On  the  other  hand,  from  Eq.  (9) , 
we  see  that  collision-induced  power  broadening 
would  occur  for  /  of  the  order  of  1JK.  Since  K »  1, 
it  is  possible  to  have  power  broadening  of  the  over¬ 
all  molecular  line  shape  while  y„  remains  close  to 
y'.  This  analysis  shows  clearly  that  collision-in- 
Aaced  power  twbarfMtttffM  not  \  varAtoGt&Ucr  of 
<  the  power  broadening  of  the  homogeneous  compo¬ 
nents,.  i  ■  1 

r  •  I 

IV.  EXPERIMENT  IN  A  COa  LASER  , 

A  simple  experiment  was  performed  using  a  COa  , 
laser  to  demonstrate  the  collision-induced  power- 
hroadenin^  effect.  Consider  a  symmetrically  con¬ 
structed  laser  operating  on  a  single  mode.  The 
holes  burned  in  the  gain  profile  in  such  a  laser  are  ' 
symmetric  about  line  Center.  If  the  overlap  of 
these  holes  is  small,  they  can  be  assumed  to  satu¬ 
rate  the  transition  independently. .  If  /  is  the  inten¬ 
sity  of  the  wave  traveling  in  either  direction  inside 
the  cavity,  it  can  be  shown8  that  Eq.  (7)  would  ap¬ 
ply  now  with  2K  replacing  K.  Thus,  for  1  «IQ  and 
7  «ku,  we  have  1  .  , 

1  ^L-  «nexpl-(0)  -  OinlVtftM)4] 

/  dz  i  +2A(///0) exp[-(w  - 

i  '• 

where 


The  output  of  the  laser  from  either  end  mirror  is 
the  intensity  of  one  of  the  traveling  waves  inside 
the  cavity  multiplied  oy  the  transmTUTVity  ot  The 
mirror.  Therefore,  the  field  intensity  inside  the 
cavity  can  be  determined  from  the  detected  output 
power  density.  Since  the  gain  inside  the  cavity  at 


the  frequency  of  oscillation  is  saturated  down  to 
the  loss  level,  we  can  write,  for  small  a^L 

ctnL  exp[-(o)  -  o»n)2/(fa<)2|  .  , 

1  +  mi/lS expHtx)  -  o>o )7(fcM)z] "  ’ 

f 

Where  L  is  the  length  of  the  active  medium  and  l  is 
the  cavity  loss  per  transit.  Rearranging,  we  ob¬ 
tain  1  • 

(afoL/1)  ekp[Mw  -  a>o)2/(few)2]  ‘ 

=  1  +  2K(I/Iq)  exp[-(o>  -  w0)z/(ku)z] .  (il) 

For  a  given  w,  Uen,  the  plot  of  (a^L/l)  exp[-(w 
-  Wq)2/ (few)2]  vs  /  is  a  straight  line  intercepting  the 
ordinate  at  unity.  The  slope  of  the  straight  line 
should  have  the  Doppler  dependence  on  in  as  re¬ 
quired  by  Eq.  (11),  The  ratio  oiqL/I  itself  is  easi¬ 
ly  calculated  from 

ctoL/l  =  exp[(iAw)7(fen)2],  , 

where  Aw  is  the  single-mode  oscillation1  width. 

The  experiment  employed  a  1.45-m-dong  cavity 
with  a  1-m  water-cooled  COz  gain  tube.  The  reso¬ 
nator  was  formed  by  a,  98%  reflecting  curved  mir¬ 
ror  (with  a  5. 7rm  radius  of  curvature)  and  a  150- 
line/mm  grating  operated  in  Littrow  (with  96%  effi¬ 
ciency),  Only  CO*  and  Nz  were  used  in  the  gain 
tube,  with  a  total  pressure  of  0. 58  Torr  (thus  en¬ 
suring  y'  <<  ku).  Gain  was  varied  by  adjusting  the 
discharge  current,  and  it  was  assumed  that  the 
proetdun.  tem  ttm  uavt*  param^t-rs  un¬ 
changed .  With  increasing  current,  the  tuning 
curves,  which  all  exhibited  the  Lamb  dip,  showed 
progressively  larger  oscillation  widths.  But  the 
Condition  I  «I0  was  fulfilled  (even  for  the  highest 
gain  used.  i  1  i 

Analysis  was  made  of  the  saturation  behavior  at 

0  t  71  Ua  wUr»  Uw  o»j:iap- 
ping  of  the  homogeneous  holes  is  indeed  insignifi¬ 
cant.  The  results  aRe  shown  in  Fig.  2,  where  the 
power  intensity  axis  is  approximately  scaled  to 
W cm-2.  When  the  slopes  of  the  straight  lines  are 
computed,  we  find  this  ratio  to  be  1.3,  in  good 
agreement  wifh  the  predicted  ratio  of  1. 37  ob¬ 
tained  from  (11).  Thus,  the  experiment  suggests 
strongly  that  the  effect  of  collision-induced  power 
broadening  does  exist  for  vibration- rotation  tran¬ 
smute  in 
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Repetitive  Passive  Q  Switching 
of  Single-Frequency  Lasers 

HOWARD  T.  POWELL  and  GEORGE  J.  WOLGA 


Abstract — A  almple  model  describing  repetitive  passive  Q  switch¬ 
ing  of  e  single-frequency  laser  with  an  inirac&vlty  saturable  absorber 
Is  presented.  Three  coupled  nonlinear  rate  equations  are  postulated 
to  describe  the  inversion  denr.lty  cf  the  amplifier,  the  inversion 
density  of  the  absorber,  and  the  photon  density  in  the  laser  cavity. 
Each  of  the  three  variables  Is  assumed  to  exhibit  a  simple,  exponen¬ 
tial  relaxation.  By  considering  the  linear  equations  describing  small 
deviations  about  the  expected  equilibrium  point,  a  condition  on  the 
various  parameters  is  found  that  describes  when  repetitive  pulsing 
will  occur.  Digital  computer  solutions  are  shown  to  verify  this 
condition.  Experiments  with  a  COj  laser  with  an  internal  COi 
absorber  provide  qualitative  agreement  with  the  model. 

I.  Introduction 

ATURABLE  absorbers  within  a  CO*  laser  cavity 
can  produce  spontaneous  repetitive  pulsing  of  the 
laser  intensity.  The  gases  SF#,  CaH«  (propylene), 
heated  CO*,  BC1*  plus  a  small  amount  of  NH„  and 
others  havo  been  used  to  provide  saturable  absorption 
(IH6J.  This  mode  of  operation  was  labeled  repetitive 
passive  Q  switching  by  its  first  observers,  Wood  and 
Schwarz  [1].  They  have  suggested  that  the  operating 
intensity  of  the  laser  must  produce  significant  saturation 
of  the  absorber  for  such  pulsing.  A  simple  model  that 
yields  a  more  detailed  condition  on  tho  laser  parameters 
delineating  when  repetitive  pulsing  will  occur  is  presented 
hero.  An  approximate  formula  for  the  pulsation  frequency 
is  also  given,  but  the  analysis  docs  not  predict  the  pulse 
height  or  width. 

A  clear  distinction  must  be  made  between  repetitive 
passive  Q  switching  and  pulsation  in  mode-locked  lasers. 
It  was  observed  that  passive  Q  switching  can  occur  in 
a  single-mode  single-frequcney  CO*  laser.  Hence  such 
Q  switching  docs  not  depend  on  the  interaction  of 
numerous  cavity  modes.  Multimode  lasers  with  an 
internal  saturable  absorber  ean  exhibit  passive  Q  switching 
and  mode  locking  concurrently  [6).  Mode  interactions 
are  neglected  and  thus  this  treatment  applies  best  to 
single-mode  lasers. 

The  phenomenon  of  output  spikes  in  solid-state  lasers 
has  been  widely  discussed.  Theoretical  considerations 
of  a  mathematical  system  of  two  rate  equations  governing 
tho  field  intensity  and  inverted  population  density  have 
predicted  that  only  damped  relaxation  oscillations  should 
occur  (7).  The  present  treatment  uses  an  independent 
third  rate  equation  to  describe  the  absorber.  This  ad- 
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ditional  degree  of  freedom  will  be  shown  to  allow  the 
possibility  ftf  repetitive  pulsation.  Shimoda  has  shown 
theoretically  that  a  supposed  cavity  loss,  whose  saturation 
follows  the  instantaneous  value  of  the  field,  can  produce 
undamped  relaxation  oscillations  [8].  A  more  general 
model  is  used  for  better  examination  of  the  system 
dynamics. 

■»  Intracavity  absorbers  are  commonly  used  in  the  pulsed 
ruby  laser  to  act  as  a  saturable  switch  that  has  the  effect 
of  compressing  the  optical  energy  into  a  single  giant 
pulse.  Theoretical  treatments  of  such  giant  pulses  (9), 
[10]  use  equations  resembling  those  applied  here.  How¬ 
ever,  since  such  studies  presumed  a  pulsed  pumping 
source,  the  need  to  consider  repetitive  pulsing  did  not 
arise. 

II.  Rats  Equations  and  Approximations 

In  the  simplest  use  of  phenomenological  rate  equations 
one  neglects  nny  spatial  dependence  and  deals  only  with 
inversion  and  photon  densities.  In  this  model  all  gain 
and  loss  parameters  are  considered  to  be  uniformly 
distributed  throughout  the  opticnl  cavity.  This  approxima¬ 
tion  is  only  valid  if  the  photon  transit  time  is  short 
comparod  to  the  time  scale  of  population  or  laser  intensity 
changes.  For  a  3-m  cavity  tho  use  of  sueh  rate  equations 
to  describe  ehanges  on  tho  microsecond  scale  is  well 
justified.  It  is  also  assumed  that  tho  entire  level  popular 
tion  interacts  with  the  radiation  field,  i.e.,  homogeneous 
broadening.  Since  this  assumption  is  not  satisfied  for 
Doppler-broadened  gaseous  absorbers,  the  model  must 
be  applied  with  some  discretion. 

A  two-level  approximation  for  both  the  amplification 
medium  and  the  absorption  medium  is  assumed.  For 
application  to  a  CO*  amplifier  the  levels  correspond  to 
specific  rotational  states  of  the  00°1  and  10°0  vibrational 
states.  Treating  the  CO*  laser  as  a  two-lcvcl  system 
despite  the  presence  of  strongly  coupled  rotational  levels 
is  certainly  an  approximation  and  is  discussed  in 
Section  VI. 

The  relaxation  of  a  population  difference  to  an  equi¬ 
librium  that  produces  cither  ait  amplification  or  absorption 
ean  often  bo  described  by  a  single  exponential  decay. 
For  many  continuously  lasing  media  tho  lower  level 
exhibits  a  much  foster  relaxation  (either  radiative  or 
collisional)  than  the  upper.  For  moderately  saturating 
intensities  the  upper  state  population  is  sizably  depressed 
while  the  lower  state  population  is  virtually  unchanged. 
The  time  dependence  of  the  population  inversion  is 
then  described  by  a  single  relaxation  rate,  which  is  the 
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relatively  slow  rate  of  the  upper  level.  Assume  that 
tho  population  dilTcrcnccs  describing  amplification  and 
Absorption  in  a  Q-swi  telling  system  arc  each  governed 
by  a  single  relaxation  rate.  These  assumptions  lead  to 
a  set  of  three  rate  equations: 

rfn, 


-jf  -  -  anxS  -  ?,(«,  -  »?) 


=  -firtaS  -  y2(n}  -  nS) 


dn i 

dl 

-j~  =  (am,  -  0nt)S  -  y.S 


(la) 


(lb) 


(lc) 


where 

», 

n, 

5 

a,  fi 


Ti.  7j 

„0  „0 

a,  i  k, 


equivalent  inversion  density  of  the  gain  medium 

if  distributed  over  the  cavity  length; 

negative  inversion  density  of  the  absorption 

medium  also  distributed; 

photon  number  density; 

field  interaction  coefficients; 

cX*/(8irfilw„)f7(x)  (0(»)  is  a  normalised  line  shapo 

factor  and  t,„„n  is  the  radiative  lifetime  between 

the  upper  and  lower  levels); 

relaxation  rates  for  the  gain  and  loss  media; 

cavity  losa  rate; 

inversion  densities  in  the  absence  of  a  laser 
field  but  with  the  presence  of  pumping  sources. 


n«i  gn, 
r/ a  '  r,/8 


Fig.  1.  Saturated  amplifier  Rain  minus  saturated  absorber  loss 
is  shown  as  a  function  of  CW  laser  intensity.  Possible  steady-state 
photon  densities  are  found  when  this  effective  saturated  Rain  rate 
equals  the  cavity  loss  rate. 


Notice  that  the  usual  pumping  terms  have  been  in¬ 
corporated  into  n”  and  nj;  hence  the  relaxation  returns 
the  population  difference  to  the  equilibrium  value  equally 
-  -  for  cither  positive  or  negative  deviations.  A  more  general 
model  with  rate  equations  for  both  the  upper  and  lower 
* 1  levels  of  each  medium  could  bo  postulated  giving  a  more 
complex  system  of  five  equations. 

By  digital  computer  solution  it  is  shown  that  the  set 
.  i  of  nonlinear  equations  (la)-(lc)  docs  exhibit  a  periodic 
pulsing  mode  for  some  sets  of  parameters.  Such  a  solution 
-f  is  visualised  as  a  closed  trajectory  or  limit  cycle  in  the 
n,,  na,  5  space.  Before  discussing  what  choices  of  pa¬ 
rameters  lead  to  a  limit  cycle  first  consider  the  expected 
steady  state  where  all  the  time  derivatives  of  (1)  are 
equal  to  zero. 

r 

mtn 

III.  Equiliiiimum  Solution  of  Rate  Equations 


I 

I 

I 
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Denote  the  expected  steady-state  values  by  fllt 
and  S.  Setting  the  time  derivatives  in  (1)  equal  to  zero, 
wo  obtain 


m,  = 


1  +  wS'/yi 


ft  2 


n2 


1  +  a.S'/y, 


(2a) 

(2b) 


an"S _ finaS 

1+mS/y,  1+fiS/y,' 


(2c) 


Equation  (2c)  determines  the  value  of  S  for  which  the 
total  saturated  gain  is  equal  to  the  cavity  loss.  A  graphical 
solution  of  (2c)  provides  considerable  insight  and  is 


shown  in  Fig.  1.  Fig.  1(a)  shows  a  monotonic  decrease 
of  the  effective  gain  with  the  steady-state  field.  Fig.  1(b) 
shows  an  initial  increase  in  the  effective  gain  due  to  the 
saturation  of  the  absorption  for  small  fields.  This  initial 
increase  can  occur  when  the  absorption  rate  fin"  divided 
by  its  saturation  parameter  y2/fi  is  greater  than  the 
corresponding  ratio  for  the  amplifier.  This  condition  is 
easily  found  by  examining  the  sign  of  the  linear  term 
in  S  in  the  expansion  of  the  effective  gain  (2c)  for  small 
fields. 

In  Fig.  1(b),  iS,  =  0  is  a  stable  solution  of  (2c)  since 
the  total  unsatumted  gain  is  less  than  the  cavity  loss. 
We  expect  S*  to  be  a  stable  solution  while  St  is  obviously 
not.  From  (lc)  notice  that  small  deviations  from  <S'j 
cause  the  photon  density  to  move  from  that  point  to 
cither  tho  stable  point  or  S*.  This  argument  is  es¬ 
sentially  that  detailed  in  [11].  Kazantsev  cl  al.  [11)  argue 
that  the  presence  of  two  stable  solutions  (S',  and  N,) 
can  cause  hysteresis  behavior  and  sudden  transitions 
from  ono  solution  to  the  other  as  laser  parameters  are 
slowly  changed.  A  simple  example  of  this  is  seen  by 
changing  the  cavity  loss  rate  y,.  For  Fig.  1(b)  there  is  a 
region  of  intensities  bordering  zero  that  arc  not  allowed. 
If  the  loss  is  decreased  from  a  large  value,  the  zero  in¬ 
tensity  solution  is  preferred  until  the  intensity  must 
jump  to  a  finite  value.  However,  if  the  loss  is  increased 
from  a  small  value,  intensity  solutions  like  arc  chosen 
until  the  field  must  drop  suddenly  to  zero.  Such  intensity 
jumping  and  hysteresis  have  been  observed  experimentally 
(12). 
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Ill  the  next  section  it  is  shown  tlmt  supposedly  stable 
solutions  considered  lierc,  nmy  give  pulsing  behavior. 

IV.  Stability  Requirements 

The  set  of  nonlinear  equations  (la)-(lc)  is  not  amenable 
to  analytic  solution.  However,  it  can  bo  determined  if 
a  presumed  equilibrium  point  gives  stable  operation. 
An  equilibrium  point  is  stable  if  and  only  if  the  linear 
system  of  equations  describing  small  deviations  about 
that  equilibrium  point,  exhibits  fundamental  modes  that 
aro  all  danced  [13].  It  will  be  shown  that  near  a  presumed 
equilibrium  point,  small  deviations  may  sometimes  grow 
in  an  oscillatory  fashion.  The  condition  for  this  happening 
is  identified  with  that  of  repetitive  passive  Q  switching. 

Strictly  speaking  an  oscillatory  instability  need  not 
load  to  a  limit  cycle,  but  our  computer  solutions  have 
shown  that  the  identification  is  valid  at  least  for  eases 
where  there  are  no  stable  equilibrium  points.  For  very 
special -eases  where  there  is  one  stable  equilibrium  point 
at  8  =  0  and  two  other  unstable  ones,  simply  studying 
small  deviations  is  not  adequate.  From  the  linear  equa¬ 
tions  the  condition  for  a  pulsing  intensity  to  occur  may 
be  correctly  identified,  but  no  information  is  obtained 
about  the  amplitude  of  the  expected  pulsation  or  its 
precise  time  dependence.  Our  numerical  solutions  show 
that  the  frequency  of  the  nonlinear  oscillation  is  ap¬ 
proximately  given  by  the  frequency  found  in  the  linear 
region. 

Following  the  general  prescription,  the  supposed 
equilibrium  is  expanded  to  first  order  in  the  deviations, 


equilibrium  solution  of  (1).  In  this  instance  nnc  must 
retain  all  the  terms  in  (4)  and  the  expected  values  of 
X  =  —ft,  X  =  —7*,  and  X  =  anj  —  f)nn1  —  yt  can  be 
easily  found.  For  amplifier  gair  minus  absorber  loss 
greater  than  the  cavity  loss,  8  =  0  is  not  a  stable  equi¬ 
librium  ns  expected.  » 

The  bracketed  coefficients  in  (fi)  are  all  functions  of 
8,  which  is  easily  varied  through  the  cavity  loss.  Presume 
that  the  eavity  loss  is  decreased  just  to  the  point  where  it 
equals  the  total  effective  gain.  The  system  roots  for  X 
are  then  —71,  —7*,  and  0.  As  »S  is  increased  from  zero  by 
decreasing  7.,  the  three  roots  must  move  smoothly  in 
the  complex  X  plane.  If  any  of  the  roots  enters  the  positive 
half-plane,  the  system  point  is  ruled  unstable. 

A  special  case  of  the  Hurwitz  algorithm  [14]  is 
applicable.  All  three  roots  of  the  cubic  polynomial 
X*  +  a,X*  +  a*X  +  is  (0|,  a*,  and  a,  real)  exhibit  negative 
real  parts  if  and  only  if  the  following  conditions  are 
satisfied:  at  >  0,  a,  >  0  and  0,0,  >  th.  For  polynomial 
(5)  the  first  condition  is  trivially  satisfied.  The  second 
condition  requires  for  stability  that 

a’n,  -  0%  ^  >  0.  (G) 

nj 

The  fact  that  certain  solutions  approaching  5  =  0  are 
not  stable  corresponds  to  the  root  at  the  origin  for  8  =  0 
moving  into  the  positive  half-plane  along  the  real  axis 
as  7.  is  decreased.  Indeed  these  unstable  solutions  can 
be  shown  to  bo  like  &  in  Fig.  1(b).  An  alternative  con¬ 
dition  for  this  type  of  instability  is 


n,(<)  =  H,  +  x(t) 
n,(t)  =  n*  -f  i/(<) 

5(0  =  S  +  s(l )  (3) 

=  —aSx  —  afl,8  —  y,x 
^  «  -0§y  -  -  7>2/ 

~  =  (ax  —  0y)S  +  (aft,  —  0n3  —  yt)s.  (4) 


Temporarily  assume  that  the  equilibrium  condition 
corresponds  to  nonzero  8  so  that  through  saturation 
afli  —  Pft*  —  7«  is  zero.  Assuming  a  solution  vector 
(x,  y,  s)  proportional  to  evi,  one  obtains  the  usual  deter¬ 
minant  condition  that  yields  a  ntbif  equation  for  X, 


sfa’A,  ^  -  d’n,  ^1  =  0. 

L  nt  m  _ 


Equations  (2a)  and  (2b)  have  been  used  to  reexpress 
such  terms  as  aS  +  7,.  The  case  <§  =  0  is  always  a  possible 


^(a»i  “  0n»)  U...11  >  0.  (7) 

One  may  easily  show  that  this  implies  the  stability 
condition  (6)  and  provides  a  nice  cross  check  of  the 
present  formalism. 

For  stability  the  third  condition  requires 

(7^  +  +  Stfyrf  _  fori)  >  0.  (8) 

It  is  asserted  that  passive  Q  switching  results  when  (8) 
is  not  satisfied.  This  third  condition  corresponds  to  one 
real  negative  root  and  two  complex  conjugate  roots 
whose  real  part  is  crossing  zero.  This  can  be  proved  by 
assuming  such  roots  and  expanding  to  find  the  corre¬ 
sponding  polynomial.  One  also  finds  that  the  square 
of  the  oscillation  frequency  is  given  by  a,  if  the  real  part 
of  the  complex  roots  is  zero.  Thus  condition  (8)  gives 
a  sinusoidally  growing  deviation  of  approximate  frequency 

a,8  ~  (^)(~)  +  S(a\  -  0%).  (9) 

Equation  (9)  is  exact  only  when  (8)  is  satisfied  os  an 
equality. 

Some  simple  predictions  can  be  mode  from  (8).  A 
necessary  requirement  for  pulsed  operation  is 

(fin°t)0y ,  >  (an?)a7„  (10) 
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The  imsatnrntcd  nm|>lificiition  and  absorption  rates 
have  been  enclosed  by !  parentheses.  We  see  from  (8) 
tliat  CW  operation  sliopld  always  be  achieved  for  8 
sufficiently  close  to  zero.  Since  the  first  term  of  (8)  con¬ 
tains  powers  up  to  there  is  also  a  maximum  8  that 
allows  pulsation. 

To  illustrate  the  typical  magnitude  of  the  terms  in  (8), 
consider  a  gaseous  molecular  medium  for  both  the  am¬ 
plifier  and  absorber.  The  effective  relaxation  rates  are 
likely  to  bo  between  10*  and  10*  s”1.  With,  for  example, 
a  3-m  cavity  and  an  average  coupling  loss  of  10  percent, 
y,  is  10r  s-1.  For  reasonable  amplification  and  absorption 
the  rates  an*  and  fin,  must  be  comparable  to  y,.  Further¬ 
more,  with  moderate  gain  saturation  the  stimulated 
emissiqn  rate  a8  is  comparable  to  the  relaxation  rate  yt. 
First  assume  that  the  inequality  (10)  is  substantially 
satisfied,  for  instance  by  a  factor  of  2.  At  typical  in¬ 
tensities  the  second  term  in  (8)  is  one  to  four  orders  of 
magnitude  greater  than  the  first  with  the  fastest  to 
slowest  relaxation  rates  mentioned.  Hence  the  upper 
and  lower  limits  for  8  bounding  the  pulsation  region 
may  bo  extremely  large  and  small,  respectively,  compared 
to  usual  laboratory  intensities.  When  (10)  is  barely 
satisfied  observable  transitions  between  pulsed  and  CW 
operation  can  occur  by  varying  8.  One  can  check  from 
(8)  that  if  the  relaxation  rate  of  cither  media  is  comparable 
to  y„  pulsing  is  difficult  to  achieve. 

Under  most  conditions,  (10)  determines  when  poising 
will  occur.  Condition  (10)  agrees  with  the  experimental 
observation  of  Karlov  et  al.,  (4J  that  passive  Q  switching 
is  promoted  by  increasing  the  relaxation  rate  of  the 
saturable  absorber.  Distinct  features  of  the  absorber 
favor  pulsation.  For  rapid  absorber  blenching,  which 
allows  a  pulse  to  rise  in  a  Q-switehed  manner,  the  absorber 
should  interact  more  strongly  with  the  field  thnn  the 
amplifier.  Then  as  the  intensity  starts  to  fall  the  absorber 
should  repump  its  population  more  rapidly  than  the 
amplifier  to  quench  that  pulse.  Eventually  the  system 
relaxes  to  a  point  where  the  pulse  notion  repeats.  Equa¬ 
tion  (10)  describes  the  balance  between  these  features. 

By  applying  the  saturated  gain  condition  (2c)  to  (9), 
it  is  observed  that  the  oscillation  frequency  of  small 
deviations  increases  with  S  provided  fit  a  is  less  than 
unity.  From  (!>)  <o3  is  given  at  very  low  intensities  by 
y,7j.  For  intensities  giving  moderate  saturation  of  the 
gain  medium  and  with  a  comparable  to  fi,  «*  is  approxi¬ 
mately  given  by  y,y e.  A  similar  range  of  nonlinear  pulsa¬ 
tion  frequencies  is  expected. 

V.  Digital  Computer  Solutions 

To  verify  the  general  conclusions  digital  computer 
solutions  of  the  nonlinear  system  (1)  were  obtained  for  a 
few  special  eases.  Dimensionless  variable  n,/n*,  n,/n\, 
and  aS/y,  were  used.  Complete  numerical  specification 
of  (1)  is  then  given  by  the  rates  «nj,  fin\,  y„  yt,  and  y. 
and  the  ratio  fi/a,  which  was  fixed  at  unity.  Primarily 
for  illustration  y,  was  fixed  at  5  X  10*  s"1,  y,  at  1  X  10* 
s'1,  an?  at  20  X  1(J)B ,  s-1  (corresponding  to  20  percent 


(b) 

Fig.  2.  (»)  Computer  solution  for  the  decay  of  an  initial  deviation 
to  a  stable-equilibrium  photon  density,  fb)  Growth  of  an  initial 
deviation  away  from  an  unstable-equilibrium  photon  density. 

gain/pass  for  a  3-m  cavity),  and  y.  at  5  X  10*  s~’.  For 
finl  greater  than  10  X  10*  s_l  (8)  predicts  pulsation  for 
nearly  all  8. 

The  initial  values  used  were  the  equilibrium  values 
fi,/n”,  nnd  a t8/y„  which  were  commuted  to  two 

significant  figures.  This  procedure  left  a  significant 
deviation  from  the  exact  equilibrium  to  test  for  stability. 
In  Fig.  2(a)  fin,  was  chosen  to  be  8  X  10*  s“',  and  thus 
the  system  is  predicted  stable  by  (8).  A  damped  sinusoid 
was  found.  In  Fig.  2(b)  fin0,  was  chosen  to  be  12  X  10*  s“' 
and  the  expected  exponentially  growing  sinusoid  was 
initially  generated.  Integrating  the  equatu  is  to  times 
beyond  1  ms  the  stable  limit  cycle  shown  in  Fig.  3  was 
eventually  reached.  Notice  that  the  pulses  generated 
look  very  much  like  reported  experimental  pulses  (2], 
(3).  The  computer-generated  pulse  amplitudes  arc  about 
10  times  the  expected  CW  level  for  our  choice  of  param¬ 
eters.  One  should  also  note  that  the  actual  pulse  spacing 
is  about  50  percent  greater  than  the  oscillation  period 
found  in  the  lineal  region. 

From  these  solutions  it  is  concluded  that  pulses  rather 
like  those  observed  experimentally  can  be  obtained  from 
(1)  using  reasonable  parameter  values.  The  assumption 
that  an  oscillatory  instability  leads  to  a  limit  cycle  is 
also  confirmed.  A  digital  computer  study  of  the  effects 
of  various  parameters  on  the  pulsation  was  viewed  ns 
impractical  due  to  the  lengthy  computation  required 
to  reach  the  stendy  state. 

VI.  Approximations  for  CO» 

In  these  experiments  a  heated  CO*  absorber  at  pressures 
less  than  5  torr  was  employed.  Since  Doppler  brondening 
applies  for  CO*  at  these  pressures  [15],  the  assumption 
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Mg.  3.  Computer  Holution  for  the  stable  pulsing  corresponding 
to  the  final  state  reached  from  Fig.  2  (b)  after  about  1  ms. 

of  homogeneous  broadening  in  our  rate-equation  model 
must  be  reconsidered.  Experiments  with  a  CO,  laser 
have  shown  that  even  when  the  Lorents  width  is  signifi¬ 
cantly  less  than  the  Doppler  width,  the  inversion  density 
versus  longitudinal  velocity  profile  in  the  presence  of 
a  strong  CW  laser  field  saturates  essentially  uniformly 
with  only  a  relatively  small  hole  (16).  This  type  of  homo¬ 
geneous  saturation  results  from  the  hole  filling  due  to 
elastic  collisions  (10).  For  a  line  saturating  in  this  uniform 
manner,  (1)  is  still  valid  provided  a  normalised  Doppler 
shape  is  used  for  g{v).  Despite  these  comments  at  least 
a  small  hole  is  burnt  in  the  CO,  inversion  at  pressures 
less  than  1  torr.  A  relatively  shallow  Lamb  dip  has  been 
observed  with  the  CO,  laser  (17),  [18]  and  in  the  next 
section  wo  report  the  observation  of  a  central  tuning 
peak  (inverted  Lamb  dip)  [19]  duo  to  a  low-pressure 
CO,  absorber. 

Tho  relaxation  of  a  CO,  laser  level  exhibits  at  least 
two  rates,  a  rapid  rotational  relaxation  rate  and  a  slower 
vibrational  relaxation  rate.  The  rotational  relaxation  rate 
is  about  107  s"'  torr"1  at  room  temperature  [20].  For  a 
typical  high-gain  CO,  amplifier  at  a  pressure  of  several 
lore  the  vibrational  relaxation  rate  is  between  10*  and 
10'  s'1  depending  on  tho  component  pressures  and  the 
discharge  conditions  [21].  Moreover  the  vibrational 
relaxation  of  the  gain  is  not  exactly  described  by  a  single 
exponential  [21].  For  a  thermal  equilibrium  CO,  absorber 
at  pressures  less  than  1  torr,  diffusion  of  molecules  out 
of  tho  laser  beam  dominates  the  vibrational  response 
[22],  [23]. 

In  agreement  with  [2]  pulse  rise  times  of  about  1  gs 
wore  observed.  This  is  considerably  longer  than  tho 
time  required  for  a  vibrational  band  of  the  amplifier 
to  adjust  to  the  field.  The  exact  deviation  of  the  set  of 
rotational  levels  from  equilibrium  depends  somewhat 
on  the  rate  of  transfer  by  inelastio  collision  between 
pairs  of  rotational  states.  If  we  assume  that  rotational 


thermalization  is  maintained,  (la)  is  modified  onl; 
slightly. 

Rotational  thermalization  requires  the  population  of  n 
rotational  level  to  be  some  constant  fraction  x>  determined 
by  the  rotational  temperature,  of  the  total  population  in 
that  vibrational  state.  If  in  (la)  the  density  n,  refers  to 
tho  total  population  in  all  rotationnl  levels,  —  (xwi)"*s' 
must  be  used  for  the  field  interaction  term.  Multiplying 
this  modified  equation  by  Xi  it  is  found  that  the  presence 
of  strongly  coupled  rotational  levels  merely  modifies  the 
rate  equation  of  a  single  rotational  level  by  reducing 
its  coupling  to  the  field  by  the  factor  x<  Thus  tho  measure*  1 
saturation  parameter  should  bo  x~*  larger  than  that 
expected  from  the  measured  vibrational  relaxation  rate 
and  interaction  coefficient  os  other  authors  have  remarked 
[22],  [24]. 

It  was  just  shown  that  for  tightly  coupled  levels  the 
saturation  caused  by  a  laser  field  is  reduced  by  simply 
sharing  tho  population  deviation.  The  calculation  of 
Section  IV  can  be  easily  redone  with  population-sharing 
factors  xi  and  x,  inserted  into  (la)  and  (lb).  Equation  (8/ 
is  changed  only  in  tho  second  term;  factors  of  xi  and  Xi 
enter  the  respective  products  of  amplifier  and  absorber 
parameters.  Thus  pulsing  is  enhanced  if  the  absorber 
is  tightly  coupled  to  fower  levels  than  the  amplifier. 

Complete  rotationnl  thermalization  of  a  CO,  absorber 
is  unlikely  at  typical  pressures  with  pulsed  operation 
(0.2-1.0  torr).  Hanst  et  al.  [2]  havo  suggested  that  only 
a  single  rotational  line  of  the  absorber  is  bleached  during 
a  pulse.  Lack  of  rotational  thermalization  of  the  absorber 
augments  tho  possibility  of  pulsation  in  two  ways.  First . 
a  single  absorber  lino  is  more  easily  saturated  by  the 
pulse  and  second,  a  rapid  rotationnl  component  is  added 
to  its  relaxation.  This  conclusion  is  substantially  the 
same  as  reached  in  tho  previous  paragraph  by  letting 
tho  absorber  be  thermalized  with  fewer  levels  than  the 
amplifier. 

If  a  equals  /?  os  is  applicable  in  these  experimental 
studies,  (10)  requires  y,  to  be  greater  than  y,  for  repetitive 
pulsing  to  occur.  In  this  work  the  absorber  was  often 
at  a  pressure  one-tenth  that  of  tho  amplifier.  Even 
accounting  for  an  increase  in  the  vibrational  relaxation 
duo  to  heating  the  absorption  cell  and  an  increase  from 
diffusion  ,\t  the'  lower  absorber  pressures,  it  is  unlikely 
that  the  vibrational  relaxation  rate  of  the  absorber 
exceeded  that  of  the  amplifier.  Lack  of  rotational 
thermalization  in  the  absorber  is  than  mandatory  for 
the  present  model  to  explain  tho  pulsing. 

VII.  Experimental  Observations 

Tho  repetitive  passive  Q  switching  of  a  CO,  lasci 
caused  by  a  heated  CO,  absorber  inside  the  laser  cavity 
was  studied  using  a  0.f)-m  gain  cell  and  a  1.5-m  absorption 
cell  separately  terminated  with  NaCl  Brewster  windows 
A  2.75-m  optical  cavity  was  formed  by  a  150-lines/mm 
grating  and  a  98  percent  reflectivity,  4-m-radius-of- 
curvature  Ge-output  mirror.  A  flowing  mixture  of  equal 
parts  CO,  and  N,  was  excited  by  a  do  voltage  to  provide 
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gain.  Several  torr  of  lie  were  often  added  for  increased 
amplification.  The  absorption  cell  could  bo  heated  to 
400°C  by  a  nichromc  wire  heater  coil.  The  laser  output 
was  monitored  by  the  dc  output  of  a  Cu:Go  detector, 
displayed  on  an  oscilloscope.  A  load  resistance  of  1  k U 
was  used  when  a  detector  response  time  of  less  than 
1  t  /us  was  desired.  The  output  mirror  was  mounted  to  a 
piezoelectric  translator  so  that  the  cavity  length  could 
be  swept  about  12  a  with  a  maximum  repetition  rato  of 
|;cn  times  per  second.  With  proper  adjustment  to  a  single 
inode  and  single  rotational  line,  distinctive  arch-shaped 
tuning  curves  separated  by  regions  of  zero  intensity 
were  produced  as  the  cavity  length  was  swept. 

Pulsed  operation  on  numerous  P  and  II  lines  could 
bo  obtained  for  absorber  pressures  between  0.2  and  1.0 
torr  and  tcmjicrnturcs  above  200°C.  The  laser  pulses 
resembled  those  pictured  by  Hanst  el  al.  (2J  and  were 
generally  separated  by  regions  of  zero  intensity.  By 
slowing  opening  an  aperture  (at  fixed  cavity  length) 
a  pulse  spacing  of  about  500  ps  near  threshold  was  con¬ 
tinuously  reduced  to  as  small  as  15  ps  as  the  average 
laser  intensity  was  increased.  The  pulse  heights  and 
pulsewidths  (1-5  gs)  were  relatively  constant  during 
this  variation.  At  the  highest  absorber  temperatures 
combined  with  the  lowest  absorber  pressures,  pulsing 
was  restricted  to  a  region  of  intensities  bordering  zero. 
By  increasing  the  intensity  further  the  pulses  came 
together  to  a  minimum  spacing,  diminished  in  amplitude, 
and  CW  oscillation  began.  At  the  lowest  absorber  tem¬ 
peratures  along  with  the  lowest  absorber  pressures, 
pulsation  was  found  at  intermediate  intensities  with  CW 
operation  occurring  at  both  higher  and  lower  intensities. 
For  less  marginal  conditions  pulsation  occurred  for  all 
intensities  achieved. 

The  dependence  on  the  average  intensity  of  both  the 
pulse  spacing  and  tho  transitions  between  pulsed  and  CW 
operation,  agree  with  the  qualitative  predictions  of  our 
model.  The  observed  range  of  pulsation  frequencies 
can  bo  used  together  with  the  approximate  values 
(YiY*  <  wJ  <  YiY.)  computed  in  Section  V  to  estimate 
the  phenomenological '  relaxation  rates  >i  and  y*.  With 
an  estimated  5  percent  cavity  loss  per  pass,  relaxation 
rates  of  KV-IO"  s_l  arc  predicted.  The  uncertainty  arises 
from  the  approximate  relation  between  «  and  the  actual 
pulsation  frequency.  These  computed  rates  arc  somewhat 
faster  than  vibrational  relaxation,  probably  indicating 
the  presence  of  a  rapid  rotational  component. 

The  dependence  of  the  pulsation  on  the  average  in¬ 
tensity  was  roughly  tho  same  if  the  cavity  length  rather 
than  tho  cavity  loss  was  changed  to  produce  the  intensity 
variation.  The  cavity-length  tuning  curve  in  Fig.  4 
Bhows  transitions  between  pulsed  and  CW  operation 
ae  mil  m  a  tuning  peak  liiaj  In  Ih®  CW  ngkuu  Hu* 
frequency  width  of  tho  peak  is  the  homogeneous  width 
of  the  CO*  absorber.  This  peak  was  only  observed  during 
CW  operation  over  a  very  limited  range  of  absorber 
temperatures  and  pressures.  Generally,  pulsing  occurred 
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Fig.  4.  Output,  intensity  vermin  cavity  length  showing  a  central 
peak  with  CW  operation.  The  intensity  hand  is  produced  by 

Kassive  Q  switching  on  a  much  shorter  time  scale.  The  pulse 
eights  are  RC  limited  by  tho  detector  load  resistance  and  some 
nonlinearity  of  tho  frequency  sweep  is  evident  (P.h,  =  0.3  torr, 
T*.  ~  300°C) 
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Fig.  5.  Intensity  versus  cavity -length  tuning  curve  showing  regions 
of  repetitive  passive  Q  switching  that  nro  well  separated  by  regions 
of  CW  operation.  The  model  presented  gives  a  qualitative  expla¬ 
nation  of  puking  domains  (P,b.  =  0.4  torr,  T. i,.  =  200*0). 

throughout  the  central  tuning  region  and  was  accompanied 
by  some  peaking  of  the  pulse  heights  about  line  center. 
The  pressure  dependence  of  the  half-maximum  width 
of  the  peaking  of  the  pulse  heights  had  a  slope  of  about. 
6  MHz/torr  and  a  zero  pressure  intercept  of  about 
2  MHz.  A  tuning  peak  produced  by  a  CO*  absorber  in  a 
pulsed  CO*  laser  has  been  previously  reported  (25). 

A  cavity-length  tuning  curve  displaying  marginal 
pulsing  is  shown  in  Fig.  5.  This  picture,  taken  near  the 
lowest  absorber  temperature  for  passive  Q  r  witching, 
shows  regions  of  pulsing  at  intermediate  intensities  with 
CW  operation  occurring  for  both  higher  and  lower  in¬ 
tensities.  Hie  central  band  of  pulsing  comes  where  a 
tuning  peak  might  otherwise  have  occurred.  The  satura¬ 
tion  of  the  absorber  is  heightened  for  frequencies  within 
the  homogeneous  width  of  its  line  center,  this  effect 
IcndlPB  bo  Uir  tuniufpudi  (JOJ  pi«»ctmrr  far  puking 
in  the  central  region  is  qualitatively  explained  by  the 
increased  absorption  saturation. 

For  CO*  absorber  pressures  above  2  torr  and  tem¬ 
peratures  below  200°C  sudden  jumps  between  zero  and 
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Fig.  6.  Intensity  versus  cavity-length  tuning  profile  showing  a 
region  of  unallowed  CW  intensities  bordering  zero.  The  jump 
distance  going  from  zero  to  a  finite  intensity  difTers  from  the 
distanco  going  from  a  finite  value  to  xcro  intensity  (P,b,  —  5.5 
torr,  T*.  “  50°C). 

finite  intensities  were  observed  by  slowly  varying  the 
cavity  loss.  This  indicates  that  the  absorber  loss  was 
saturating  more  rapidly  with  intensity  than  the  gain 
(Section  III).  Such  intensity  jumps  were  also  found 
by  sweeping  the  cavity  length  (Fig.  0).  Intensity  jumping 
and  pulsing  did  not  occur  together.  This  agrees  with  (6) 
and  (8)  that  imply  tliatyi  must  be  less  than  y,  for  intensity 
j  umps  but  greater  than  y,  for  pulsing  (provided  a  equals  $). 

The  addition  of  buffer  gases  to  cither  the  amplifier  or 
absorber  produced  interesting  effects.  Adding  larger 
amounts  of  He  to  the  amplifier  than  needed  to  produce 
the  best  gain,  enhanced  the  onset  of  pulsing.  However, 
by  adding  about  0.5  torr  of  either  He,  Nt,  or  H,0  to 
flic  absorber,  pulsation  was  stopped.  Both  observations 
aro  consistent  with  the  view  that  these  gases  aid  rotational 
tlicrmalisation.  In  the  amplifier,  rotational  thermalization 
augments  pulsation  but  in  the  absorber  it  hinders  pulsa¬ 
tion  (Section  VI). 

A  new  domain  of  passive  Q  switching  using  a  COt 
absorber  was  observed.  At  an  absorber  pressure  of  5  torr 
and  an  absorber  temperature  of  about  50°C,  the  addition 
of  0.3  torr  II20  vapor  caused  the  laser  to  pulse  regularly. 

•  'Rising  was  limited  to  a  narrow  region  of  intensities 
icar  threshold  operation.  The  pulses  were  about  500  pa 
aide,  separated  by  approximately  2  ms,  and  about  5 
l  imes  higher  than  the  expected  CW  intensity.  Wc  speculate 
i  hat  the  addition  of  II, 0  vapor  barely  caused  the  vibra¬ 
tional  relaxation  rate  of  the  absorption  to  exceed  that 
of  the  gain,  thus  allowing  weak  pulsing. 

VIII.  Summary 

A  model  has  been  given  that  explains  some  aspects  of 
repetitive  passive  Q  switching.  Equations  (8)  and  (10) 
could  bo  useful  in  the  selection  of  saturable  absorbers 
to  produce  pulsing  of  CO,  and  other  lasers.  It  is  likely 
that  repetitive  passive  Q  switching  can  be  based  on 
media  with  varied  relaxation  mechanisms. 
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Optical  Saturation  of  a  Single  Vibration-Rotation  Transition  in  the 

vi  Fundamental  of  SF«*  f 

N.  Djeu$  and  G.  J.  Wolga 
Cornell  University,  Ithaca,  New  York  14850 
(Received  22  April  1970) 

A  well-resolved  vibration-rotation  transition  in  the  i/j  infrared  fundamental  of  SF«  has  been  revealed 
at  200°K  by  using  the  tuned  laser  spectroscopy  technique.  The  frequency-dependent  saturation  behavior  of 
that  transition  has  been  studied  at  low  pressures,  and  the  results  fitted  to  theoretically  predicted  curves. 

From  the  pressure  dependence  of  the  saturation  parameter,  the  quantity  (y'/y)Aai,  which  is  the  product 
of  the  Einstein  coefficient  and  the  ratio  of  the  dephasing  collision  rate  and  the  rotational  relaxation  rate, 
was  determined  to  be  53  s:  c~‘.  This  is  considered  as  a  good  upper  bound  for  the  Einstein  Aab  coefficient 
itself. 


INTRODUCTION 

In  the  past  two  years,  a  large  number  of  experiments 
have  been  performed  on  SF«  in  connection  with  its 
absorption  at  the  10.6-/Z  wavelength  of  the  C02  laser. 
They  range  from  the  observations  of  self-induced  trans¬ 
parency,1  of  photon  echo, 2,3  and  of  transient  nutation4  6 
to  investigations  of  the  inverted  Lamb  dip67  and  studies 
of  steady-state  saturation.8  9  But  in  spite  of  the  ex¬ 
tensive  literature,  our  knowledge  of  the  physical 
properties  of  that  molecule  is  far  from  being  complete. 
As  a  case  in  point,  the  figures  reported  for  the  Einstein 
coefficient  for  the  vs  band,  which  is  responsible  for 
the  absorption  of  radiation  near  10.6  jt,  vary  by  as  much 
as  two  orders  of  magnitude.  Further  experiments  on 
that  gas  were  therefore  deemed  desirable. 

Quite  apart  from  the  above-mentioned  developments, 
we  have  derived  elsewhere10  a  saturation  formula  for 
molecular  gases  which  predicts  that  in  a  limited  pressure 
region  the  optical  saturation  of  a  vibration-rotation 
transition  should  be  identical  to  that  of  a  two-level 
atomic  transition.  Optical  saturation  measurements 
carried  out  in  that  pressure  region  ought  to  yield 
directly  the  product  (y'/y)Aab,  where  7'  is  the  de¬ 
phasing  collision  rate  and  7  the  rotational  relaxation 
rate.  Since  7'  is  expected  to  be  only  slightly  greater 


than  7,  one  can  obtain  a  good  upper  bound  for  Aab  from 
such  saturation  studies. 

Tuned  laser  spectroscopy  revealed10  that  at  200°K 
and  f<  r  sufficiently  low  pressures  there  lies  within  the 
Doppler  width  of  the  P{  18)  transition  of  the  C02 
10.4-ji  band  a  well-resolved  SF«  vibration-rotation 
transition  belonging  to  the  vs  fundamental.  It  was  thus 
possible  to  test  the  validity  of  the  saturation  formula  in 
SFe  over  a  restricted  pressure  range.  Furthermore,  the 
saturation  parameter  measured  as  a  function  of  pressure 
determined  {y'  /y)Aab  unambiguously. 

The  saturation  studies  reported  here  differ  from  the 
earlier  ones8,9  in  that  we  operated  under  temperature 
and  pressure  conditions  which  insured  single-line  ab¬ 
sorption  in  SFe  and  that  the  C02  laser  radiation  em¬ 
ployed  was  frequency  tuned.  This  procedure  permitted 
the  saturation  formula  to  be  tested  in  detail. 

EXPERIMENTAL  ARRANGEMENT 

A  schematic  of  the  apparatus  used  for  making  small 
signal  and  saturated  absorption  measurements  is  shown 
in  Fig.  1.  The  1-m  water-cooled  active  laser  discharge 
column  operated  on  a  flowing  mixture  of  12.5%  C02, 
25%  N,,  and  62.5%  He  at  a  total  pressure  of  3  torr. 
The  discharge  current  could  be  adjusted  to  give  a 
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Fig.  1.  Experimental  arrangement 
for  absorption  measurements. 
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variable  power  output.  The  cavity  was  formed  by  a 
curved  mirror  Mi  (5.7-nt  radius  of  curvature,  98% 
reflecting),  a  grating  G  (150-lines/mm  gold  replica 
blazed  for  8  n)  which  permitted  selective  oscillation  on 
any  one  particular  vibration-rotation  transition,  and  a 
flat  mirror  M2  (40%  transmitting).  The  total  cavity 
length  was  1.60  m,  and  fine  tuning  was  accomplished 
by  mounting  Mi  on  a  piezoelectric  element  driven  by  a 
voltage  ramp  generator  capable  of  sweeping  with 
various  ramp  slopes.  The  angle  formed  by  Mi,  G,  and 
M2  was  so  chosen  that  •  -hen  oscillation  took  place  on 
any  of  the  9.4-  or  10.4-/1  band  transitions  with  Mi  and 
Mi  as  end  mirrors  (the  grating  was  not  used  in  Littrow), 
Mi  and  G  alone  did  not  form  a  high  Q  cavity  for  any 
of  the  laser  transitions.  With  the  use  of  an  iris  near 
Mi,  it  was  possible  to  limit  oscillation  to  only  the 
fundamental  transverse  mode.  The  output  beam,  after 
going  through  an  attenuator  A,  was  deflected  by  a  flat 
gold  mirror  Mi,  mod  lated  by  a  chopper  wheel  C,  and 
detected  by  a  gold-doped  germanium  detector  D 
cooled  to  liquid-nitrogen  temperature.  The  output  of 
the  detector,  after  being  amplified  100  times,  was 
displayed  on  an  oscilloscope  for  absolute  measurements. 
It  was  also  sent  into  a  tuned  amplifier,  whose  output, 
after  rectification,  was  fed  into  a  chart  recorder. 

A  1-in.  i.d.  Pyrex  pipe  of  82  cm  over-all  length  with 
NaCl  Brewster  windows  served  as  the  absorption  cell. 
An  open  jacket  around  the  cell  could  be  filled  with  dry 
ice  to  cool  the  entire  cell,  except  for  4  cm  on  each  end, 
to  200°K.  The  power  detected  with  the  empty  cell  in 
and  out  of  the  path  of  the  laser  beam  agreed  to  within 
2%  or  3%.  Integrated  loss  was  deduced  by  comparing 
the  unattenuated  signal  with  the  transmitted  signal 
when  the  cell,  filled  to  a  given  pressure,  was  intercepting 
the  beam. 


SMALL  SIGNAL  ABSORPTION 

The  full  width  between  1/e  points  for  a  Doppler 
broadened  vibration-rotation  line  in  SF#  near  10.6  ju 
at  200°K  is  29.4  MHz.  On  the  other  band,  the  single¬ 
frequency  tuning  range  for  the  C02  laser  used  was  94 
MHz.  Therefore,  SF#  transitions  whose  line  centers  are 
only  a  few  tens  of  MHz  from  C02  laser  line  centers 
could  be  studied  in  detail  by  the  tuned  laser  spectros¬ 
copy  technique.  The  unattenuated  single-frequency 
tuned  laser  output  was  compared  with  the  transmitted 
signal  after  the  beam  traversed  the  absorption  cell. 
Provided  the  power  level  of  the  incident  radiation  was 
much  lower  than  that  needed  for  saturation,  taking 
the  logarithm  of  the  ratio  of  the  two  signals  at  various 
frequencies  yielded  directly  the  line  shape  of  the  ab¬ 
sorbing  SF#  transition. 

The  upper  curve  in  Fig.  2  is  the  incident  intensity  at 
the  minimum  noise-free  level  (about  1  mW/cm2)  as  a 
function  of  frequency  of  the  C02  laser  operating  on  the 
P(18)  transition  of  the  10.4-ju  band.  The  asymmetric 
lower  curve  is  the  transmitted  signal  after  passage 
through  the  absorption  cell  containing  a  mixture  of 
0.007  torr  of  SF#  and  0.11  torr  of  He  at  200CK.  The 
partial  pressure  of  He,  while  not  high  enough  to  cause 
any  appreciable  collision  broadening,  was  sufficient  to 
insure  unsaturated  absorption.10  The  vertical  lines  are 
drawn  in  to  match  frequency  points  on  the  incident 
power  curve  to  corresponding  ones  on  the  transmitted 
signal.  They  are  evenly  spaced  by  3.75  MHz  but  have 
no  fixed  reference.  The  SF#  line  center  is  seen  to  be 
displaced  from  the  C02  line  center  by  about  10  MHz,  in 
agreement  with  the  Lamb  dip  spectroscopy  result.7 

When  the  logarithm  of  the  ratio  of  the  incident  signal 
over  the  absorbed  signal  for  different  points  along  the 
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Fig.  2.  Small  signal  incident  and  transmitted  power  curves  as 
functions  of  frequency. 


frequency  axis  is  taken,  we  obtain  the  data  displayed 
in  Fig.  3.  This  can  be  fitted  very  closely  to  a  Gaussian 
curve  of  30-MHz  full  width  between  1/ e  points  if  the 
base  line  for  that  curve  is  taken  slightly  above  the  zero 
level  of  the  actual  absorption  data.  The  displacement 
is  interpreted  as  an  uniform  background  absorption 
attributable  to  the  wings  of  adjacent  lines.  It  is  about 
8%  of  the  total  absorption  at  line  center  under  the 
experimental  conditions  given  above.  The  single  line 
absorption  coefficient  a0  was  found  to  be  3.6  cm-1  •  torr-1 
for  SF,  at  200°K. 

A  similar  fit  for  data  taken  with  the  same  mixture 
but  at  room  temperature  revealed  the  background  to  be 


Fig.  3.  Experimental  line  shape  fitted  to  a  Gaussian  curve. 


35%.  The  additional  overlapping  is  attributed  to  the 
rising  importance  of  hot  band  absorption  at  room 
temperature.  Furthermore,  when  the  He  partial 
pressure  was  increased  to  1  torr,  the  single-line  feature 
was  almost  totally  obscured.  These  effects  are  com¬ 
patible  with  the  finding  of  Rabinowitz  et  aV  that  the 
typical  separation  between  vibration-rotation  transi¬ 
tions  near  10.6  y  is  of  the  order  of  the  Doppler  line- 
width. 

SATURATION  RESULTS 

We  have  shown  elsewhere10  that  the  saturated  in¬ 
cremental  loss  for  a  traveling  wave  interacting  with  a 
vibration-rotation  transition  in  a  molecular  gas  can  be 


Fig.  4.  Integrated  loss  at  line  center  vs  input  intensity  com¬ 
pared  to  computer  generated  theoretical  curve  (P =0.01 7  torr, 

r  =  200°K). 


written  in  Gaussian  units  as 

1  dl  _  _ otp  RefF(x-fty) _ 

7  Jz  ~  {\+I/Io)w+K{l/h)  RefF(x-fiy) 

with 

W(x+iy)  =  2r~m  exp[—  (x+iy)2] 

X  [  exp (-Fjdf, 

where 

x=  (McV2*r)‘«[(  «,-«)/«], 
y=  (Jlfc2/2^r)1'2[l+(///o)]1/27'/w, 
/o=m'fi2/4x/'2, 

k  =  {rMcyikTy’W  ( r0+ r»)/2«wr.r»]. 

Here  I  is  the  field  intensity,  70  the  saturation  parameter, 
z  the  direction  of  the  traveling  light  beam,  w  its  angular 
frequency,  op  the  unsaturated  (small  signal)  incre¬ 
mental  loss,  wp  the  angular  frequency  of  the  absorbing 
transition  at  line  center,  T0  and  rV  the  relaxation  rates 
of  the  upper  and  lower  vibrational  levels,  respectively, 
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Fig.  5.  Integrated  loss  at  *=0.5  vs  input  intensity  com-  1 
pared  to  Computer  generated  theoretical  curve  (P=0.017  torr, 

T  =i200°K) . 


And  1/ n  the  ratio  of  the  absorbing  rotational  population 
to  the  associated  vibrational  population  at  thermal 
equilibrium.  The  average  dipole  moment  squared  is 
given  by 

£  £  I  I  P.  I  hi)  I*. 

«-l  i- 1 

where  and  bj  art  magnetic  sublevels  of  the  upper  and 
lower  absorbing  rotational  states,  and  g„  and  gb  their 
respective  degeneracies.  The  degenerate  levels  to  be 
summed'  over  are  those  with  different  M  and  K  quan¬ 
tum  numbers  but  the  same  J.  If  Coriolis  and  other 
interactions  remove  the  K  degeneracy  partially,  then 
the  summation  is  over  those  sublevels  which  remain 
unsplit. 

The  saturation  formula  given  above  is  a  modified 
version  of  the  one  derived  by  Close.11  The  modification 
becomes  necessary  for  vibration-rotation  transitions 
.because  rotational  relaxation  effectively  constitute 
variable  pump  rates  for  the  particular  rotational  levels 
which  are  interacting  with  the  radiation.  In  the  low- 


Fig.  6.  Integrated  loss  at  *  =  1.0  vs  input  intensity  compared 
to  computer  generated  theoretical  curve  (/>=0.017  torr,  7"= 
200°K). 


i  , 

pressure  limit  (less  than  0.1  torr  for  SF„  at  200°K), 
k«l  and  the  saturation  formula  reduces  to12 


1  dl  ao  RelF  (x-\-iy) 

id*  =  C i+tW" ' 

Since  we  know  a0  from  our  own  measurements  and  7' 
from  Ref.  2,  the  quantity 

(V/tM<*=  (64irV/AX37)/>2 

can  be  deduced  from  the  experimentally  determined 
integrated  loss  as  a  function  of  incident  laser  intensity. 

Saturation  measurements  were  made  in  pure  SF8  at 
200°K  for  pressures  of  0.009,  0.0’3,  0.017,  and  0.022 
torr.  Input  and  output  power  curves  similar  to  those 
shown  in  Fig.  2  were  obtained.  The  integrated  loss  at 
line  center  as  well  as  at  othdr  frequencies  was  readily 
found  by  plotting  the  ratio  of  the  two  curves  at  various 
points  along  the  abscissa.  Also,  for  each  given  pressure, 
a  computer  was  used  to  generate  the  integrated  loss  at 
line  center  as  a  function  of  incident  power  intensity  in 
units  of  70.  A  fit  of  the  experimental  data  to  the  theo¬ 
retical  curve  then  ascribed  a  unique  value  to!  7o  at  each 
pressure.  Such,  a  fit  for  0.017  torr  of  SF»  is  exhibited  in 
Fig.  4.  1  1  1 

Using  the  70  determined  in  this  fashion,  saturation 
curves  were  alsq  generated  (now  as  a  function  of  in¬ 
tensity  in  units  of  milliwatts  per  square  centimeters) 
for  x=  0.5  and  1.0  at  0.017  torr.  They  are  shown  in  Figs. 
5  and  6,  together  with  the  experimental  results  cor- 


Table  I.  The  Einstein  A  oh  coefficient  as  given  (or  implied)  by 
various  authors. 


A  at 

-  -  -  -i 

Ref. 

67  sec-1 
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•  C.  K.  Rhodes  and  A.  Stoke,  PhyS.  Rev.  184,  25  (1969). 
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rected  for  background  absorption.  The  agreement  sug: 
gests  the  correctness  of  the  reduced  saturation  formula 
in  its  frequency  dependence. 

When  the  parameter  70  is  plotted  as  a  function  of 
pressure  squared  in  Fig.  7,  the  points  fall  on  a  straight 
line  passing  through  the  origin.  This  is  to  be;  expected, 
because  both  y  and  y'  are  proportional  to  the  pressure. 
From  the  slope  of  the  lipe,  we  calculate  (y'/y)Aab  to 
be  53  sec-1.  Since  nearly  all  phase-interrupting  col¬ 
lisions  involve  a  change  in  the  rotational  quantum 
number  in  a  molecular  gas,  the  ratio  y'/y  should  be 
greater  than  but  very  close  to  unity.  Therefore,  the 
rate  given  above  is  regarded  as  a  good  upper  bound 
for  Aabr 

Previous  appraisals  of  the  transition  moment,  lead  to 
the  values  for  the  Einstein  A a»  coefficient  given  in 
Table  I.  The  first  was  estimated  from  pqlse  trans¬ 
mission,  while  the  second  and  third  were  obtained  from 
integrated  band  intensity.  The  fourth  figure  came  from 
the  transient  nutation  effect  and  the  last  from  the  self- 
induced  transparency  experiment.  The  large  dis¬ 
crepancies  between  these  last  two  and  our  own  upper 
bound  require  some  comments. 

The  inverse  Lamb-dip  experiment,7  as  well  as  a 
recent  high-resolution  spectroscopic  study  of  §F« 
using  a  tunable  diode  laser,13  shows  quite  conclusively 
that  the  K  degeneracy  is  indeed  removed  by  the  Coriolis 
interaction.  With  this  assumption  and  the  knowledge 
of  the  statistical  weights  of  the  rotational  sublevels 
derived  from  symmetry  considerations,  one  readily 
calculates  the  population  densities  of  these  sublevels. 
The  results  show  that  if  the  absorption  of  the  P(18) 
COj  line  is  due  to  one  of  the  sublevels  with  the  largest 
populations,  the  measured  small  signal  absorption 
coefficient  would  lead  to  an  Aab  coefficient  in  agree¬ 
ment  with  that  deduced  from  our  saturation  experi¬ 
ment.  Thus,  i  the  small  signal  absorption  coefficient 
alone  sets  a  lower  bound  on  Aab  which  is  about  two 
orders  df  magnitude  larger  than  those  given  by  the 
self-induced  transparency  and  the  transient  nutation 
experiments.  We  therefore  hold  the  last  two  numbers  in 
Table  I  to  be  highly  suspect. 


i 


CONCLUSION 


We  coi'dude  by  mentioning  some  deficiencies  as  well 
as  me  its  of  studying  molecular  systems  by  means  of 
saturation  measurements.  The  two  major  drawbacks 
are  that  one  must  know  the  dephasir.g  collision  rate 
before  an  analysis  of  the  data  can  be  made  and  that  the 
procedure  yields  not  one  single  physical  constant,  but 
the  product  ( y'/y)Aab .  This  restriction,  however,  also 
proves  to  be  an  asset  in  that  it  introduces  a  degree  of 
flexibility.  The  knowledge  of  either  y  or  A,*  enables  a 
calculation  of  the  other.  In  general,  Aab  is  more  ac¬ 
cessible,  and  we  see  our  technique  as  an  accurate  way  of 
determining  the  rotational  relaxation  rate  in  compli¬ 
cated  molecular  systems.  Furthermore,  ( y'/y )  Aab  is  a 
very  good  upper  bound  for  Aab  itself. 
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Abstract — Saturated  gain  profile  measurements  in  COi  at  pres¬ 
sures  where  a  Lamb  dip  appears  in  the  laser  intensity  reveal  that 
individual  rotation-vibration  transitions  saturate  with  a  strong 
Doppler  component  in  addition  to  the  Lorentzian  holes  expected  for  a 
Doppler  broadened  line.  This  mixed  inhomogeneous-homogeneous 
behavior  represents  cross-relaxation  effects  induced  by  collisions 
that  change  the  molecule’!!  velocity.  A  rate  equation  formulation  of 
collision  effects  elucidates  the  role  of  elastic  and  rotational  thermaliz- 
ing  collisions  on  saturation  and  Lamb  dip  formation  in  molecular 
lasers. 

Introduction 

IN  THE  absence  of  collisions  a  Doppler  broadened 
transition  that  is  subjected  to  a  resonant  monochro¬ 
matic  field  will  saturate  by  the  mechanism  of  spectral 
hole  burning  [1].  This  saturation  behavior  results  in  the 
appearance  of  the  Lamb  dip  [2]  in  the  intensity  of  a 
standing  wave  laser  as  its  frequency  is  tuned  through 
the  line  center  frequency.  Since  the  dip  width  is  the 
radiative  interaction  width  of  individual  molecules 
exclusive  of  Doppler  broadening,  information  on  collision 
effects  can  be  obtained  by  investigating  the  pressure 
dependent  behaviour  of  the  dip. 

Collision  broadening  studies  in  atomic  Ne  utilizing  the 
Lamb  dip  technique  have  been  based  upon  phenomeno¬ 
logical  modifications  of  Lamb’s  collisionless  equation  for 
the  laser  intensity.  The  simple  procedure  of  introducing 
a  density  dependence  into  the  radiative  line  width 
parameter  y'  (total  phase  interruption  rate)  was  used 
by  Smith  [3]  to  study  the  3s,  — »  2p«  transition  in  Ne  at 
X  =  0.6328  ii.  In  a  different  approach,  Sz&ke  and  Javan  [4] 
used  two  collision  parameters,  the  hard  and  soft  line 
widths,  to  characterize  the  Lamb  dip  in  the  2aa  — *  2 p4 
transition  in  Ne  at  X  =  1.15  ft. 

The  Lamb  dip  has  also  been  observed  in  rotation- 
vibration  transitions  of  the  H,0  [5]  and  CO,  [6]  molecular 
laser  systems.  Attempts  to  derive  collisional  information 
for  these  molecules  from  the  Lamb  dip  have  involved  the 
use  of  both  the  one  and  two  parameter  expressions  for 
the  laser  intensity  developed  for  atomic  transitions 
mentioned  above.  Despite  the  outward  similarity  of  the 
Lamb  dip  phenomena  in  atomic  and  molecular  lasers, 
the  collision  processes  that  affect  saturation  are  not  the 
same  for  atoms  and  molecules  and  furthermore,  the  extent 
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to  which  collisions  influence  the  saturation  behavior  in 
atomic  and  molecular  transitions  are  different  due  to 
large  differences  in  level  lifetimes.  The  description  of 
collisions  in  molecular  systems  by  the  use  of  atomic 
formulas  in  [5]  and  [6]  is  therefore  not  a  prior i  justifiable. 

We  present  in  this  paper  a  study  of  the  COa  laser 
that  elucidates  the  role  of  collisions  in  saturation  and 
Lamb  dip  formation  in  molecular  lasers.  The  experimental 
study  focuses  attention  on  the  changes  induced  in  the 
velocity  distribution  of  the  levels  undergoing  saturation. 
This  information  is  obtained  by  measuring  the  frequency 
profile  of  the  saturated  gain  [7]  on  individual  transitions 
by  the  simultaneous  use  of  two  monochromatic  fields.  An 
intense  fixed-frequency  field  saturates  the  transition  while 
a  weak  frequency-tunable  field  measures  the  change  in 
the  profile  of  the  optical  gain.  A  rate  equation  theory  is 
developed  to  describe  the  collision  effects. 

Saturation  of  Molecular  Transitions 

The  fact  that  Doppler  broadened  transitions  can 
saturate  in  a  manner  that  is  different  from  the  simple  hole 
burning  process  is  demonstrated  by  the  results,  to  be 
presented  shortly,  of  the  experiment  shown  in  Fig.  1 
where  the  saturation-induced  change  in  the  frequency 
profile  of  the  gain  in  laser  L.  is  measured  by  a  single  pass 
gain  probe  of  the  lasing  medium. 

Both  lasers  were  single  frequency  oscillators  operated 
in  the  TEM0o,  mode.  The  CO,  +  N,  mixture  in  L,  was 
set  at  pressures  below  1  torr  where  the  Lamb  dip  is 
observable  [8],  while  the  gas  mixture  in  L,  was  set  at 
pressures  above  1  torr  to  produce  a  smooth  tuning  curve. 
The  probe  beam  was  injected  into  L,  by  reflection  off  an 
intracavity  NaCl  Brewster  angle  plate  and  its  return 
back  into  L,  was  prevented  by  operating  the  two  lasers 
in  orthogonal  polarizations  in  conjunction  with  a  pair 
of  Ge  Brewster  plates  within  the  L,  optical  cavity.  A 
second  pair  of  Ge  Brewster  plates  placed  before  the 
detector  served  to  reject  the  small  amount  of  L,  radiation 
scattered  from  the  internal  plates.  The  differential  gain 
(DG),  defined  as  the  unsaturated  gain  minus  the  saturated 
gain,  was  obtained  directly  by  periodic  interruption  of 
the  L,  lasing,  and  phase  sensitively  demodulating  the 
amplified  probe  signal.  The  saturating  frequency  oi,  can 
be  held  at  any  desired  point  on  the  transition  by  the 
L,  PZT  ceramic.  The  probe  frequency  is  tuned  across 
the  transition  by  a  linear  ramp  voltage  applied  to  the 
L,  ceramic. 

While  the  simple  hole-burning  picture  seems  to  ade¬ 
quately  describe  the  saturation  behaviour  in  the  2p«  level 
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Fig.  1.  Apparatus  for  probe  of  saturated  gain.  B — NaCl  beam 
splitter,  C — oollimation  optics.  D — deflecting  mirrors,  0 — 
diffraction  grating,  M — spherical  mirror,  P — Ge  Brewster  angle 
plates. 
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of  Ne  [7],  the  saturated  gain  in  CO*  on  the  P(18)  (00°1)- 
(10°0)  transition  at  a  pressure  of  0.4  torr  reveals  a  different 
behavior.  See  Fig.  2(a)--(e).  In  addition  to  the  Lorentzian 
holes  expected  of  an  inhomogeneously  broadened  line,  a 
Doppler  saturation  in  the  form  of  a  much  broader  and 
stronger  component  that  displays  the  Doppler  width  of 
the  transition  occurs.  Whereas  an  inhomogeneously 
broadened  line  responds  to  a  monochromatic  field  over  a 
frequency  width  that  is  small  compared  to  the  spectral 
line  width  and  a  homogeneously  broadened  line  responds 
to  a  resonant  field  over  the  entire  spectral  width,  the  CO, 
transitions  (Fig.  3)  behave  in  a  mixed  inhomogeneous- 
homogeneous  manner  despite  the  fact  that  the  radiative 
interaction  width  y'  is  small  compared  to  the  Doppler 
width.  The  condition  y'  «  Doppler  width  is  therefore 
not  a  sufficient  condition  for  the  spectral  line  to  saturate 
inhomogeneously  [9].  A  little  thought  will  show  that  the 
additional  Doppler  saturation  is  the  result  of  spectral 
cross-relaxation  processes  induced  by  collisions  that 
change  the  molecular  velocity.  These  collisions  cause 
diffusion  within  the  Doppler  distribution  and  act  to 
maintain  the  equilibrium  line  shape  against  radiative 
distortions. 

This  can  be  understood  in  the  following  way.  In  the 
absence  of  velocity  diffusion,  only  those  molecules  whose 
Doppler  shifted  emissions  fall  within  a  width  2  y'  about 
the  frequency  <■>,  can  interact  appreciably  with  a  field 
at  frequency  The  DG  of  such  a  nondiffusing  distribu¬ 
tion  will  therefore  consist  of  Lorentzian  holes.  When 
velocity  diffusion  occurs,  molecules  located  further  than 
a  width  y'  from  «.  now  possess  a  probability  of  increasing 
their  interaction  with  the  field  through  collisions,  which 
move  them  closer  to  w,.  The  probability  that  a  molecule 
will  augment  its  interaction  with  the  field  in  this  manner 
can  be  seen  to  be  an  increasing  function  of  the  ratio 
between  the  diffusion  rate  (inverse  velocity  diffusion  time) 
across  the  distribution  to  the  rate  at  which  molecules 
leave  the  distribution  (inverse  level  lifetime).  If  the  level 
lifetime  is  short  compared  to  the  diffusion  time  then  a 


Fig.  2.  (a)-(d)  Frequency  dependence  of  the  DG  in  L,  at  0.4 

torr  pressure  (top  tracings)  with  both  lasers  oscillating  on  the 
P(18)  (00°1)-(1040)  transition.  Bottom  tracings  show  the  position 
of  the  fixed  L,  frequency  relative  to  its  tuning  curve  for  which  the 
corresponding  DG  tracings  were  made.  Vertical  arrows  locate 
the  Lorentzian  holes  whose  asymmetric  strength  is  due  to  loss 
anisotropies,  (e)  DG  on  the  P(20)  (00°l>-(10o0)  transition  with 
L,  oscillating  on  the  P(18)  transition. 

molecule  will  have  a  low  probability  of  reaching  the 
neighborhood  of  u,  before  leaving  the  level.  In  this  case 
the  distribution  of  Lorentzians  comprising  the  Doppler 
broadened  line  can  be  considered  almost  static.  This 
behavior  is  characteristic  of  transitions  between  short 
lived  atomic  levels.  If,  however,  the  level  lifetime  is  long 
compared  to  the  velocity  diffusion  time,  then  the  molecule 
will  suffer  several  velocity-changing  collisions  within  its 
lifetime  and  the  probability  of  reaching  the  neighborhood 
of  w,  and  interacting  with  the  field  is  enhanced.  This 
behavior  is  representative  of  transitions  between  long- 
lived  rotation-vibration  transitions.  Velocity-changing 
collisions  therefore  induce  saturation  across  the  entire 
velocity  distribution  and  qualitatively  explain  the 
presence  of  the  Doppler  saturation  component.  The 
appearance  of  Doppler  saturation  in  the  CO,  levels  and 
its  apparent  absence  in  the  Ne  2p4  level1,  reflects  the  fact 
that  the  collisionally  induced  cross-relaxation  rate  is  much 
faster  than  the  slow  vibrational  decay  rate  for  the  low 
lying  levels  of  CO,  while  the  collisionally  induced  cross- 
relaxation  rate  is  slower  than  the  decay  rates  for  the 
pertinent  levels  of  Ne. 

Collisions  that  change  only  the  velocity  are  but  a 
partial  specification  of  the  total  cross-relaxation  process. 
In  molecular  systems  with  rotational  fine  structure,  a 
molecule  with  velocity  f,  in  rotational  level  j  can  migrate 
to  a  different  velocity  f,  within  the  same  level  j  via  two 
collisions,  the  first  of  which  carries  the  molecule  into  any 
of  the  other  rotational  levels  j',  and  the  second  one 
carrying  the  molecule  back  into  level  j.  Spectral  cross- 

1  The  upper  3«,  level,  which  was  not  probed  in  [7],  may  contain 
a  nonnegtigible  Doppler  component  due  to  its  longer  lifetime. 
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Fig.  3.  Saturation  behavior  of  spectral  lines  with  different  broad¬ 
ening  properties.  Dotted  curve  is  the  unsaturated  spectral  profile. 
Solid  curve  is  the  saturated  profile  for  a  monochromatic  wave  at 
line  center. 
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relaxation  effects  in  rotation-vibration  levels  must 
therefore  include  both  the  elastic  velocity-changing  and 
inelastic  rotational  level-changing  collision  processes. 

A  probe  of  the  adjacent  P( 20)  transition  with  L, 
oscillating  on  the  P(18)  transition  [Fig.  2(e)]  shows  that 
collisional  coupling  causes  the  entire  rotational  distribu¬ 
tion  to  saturate  wher.  one  pair  of  levels  is  being  saturated. 
The  Doppler  saturation  within  the  pair  of  levels  being 
stimulated  and  the  saturation  of  the  rotational  distribution 
are  thus  quite  similar  phenomena.  Both  are  cross-relaxa¬ 
tion  or  equivalently  cross-saturation  effects  brought  about 
by  collisions.  In  the  former  case  the  cross  saturation 
occurs  between  different  portions  of  the  same  spectral 
line  while  in  the  latter  case  cross  saturation  occurs  between 
different  portions  of  the  rotational  distribution. 

The  pressure  dependence  of  the  mixed  saturation 
behavior  in  C02  was  measured  using  traveling  wave 
saturation  in  an  oscillator  amplifier  experiment  shown  in 
Fig.  4.  Here,  lasers  L,  and  L,  supply  the  saturating  field 
and  the  probing  field,  respectively,  for  the  amplifier  LA. 
A  flowing  premixed  ratio  of  50  percent  CO*  and  50  per¬ 
cent  N,  gases  dc  excited  along  the  1.3-m-discharge  length 
constituted  the  amplifying  medium.  The  three-mirror 
configuration  for  both  lasers  was  chosen  to  avoid  the 
strongly  divergent  beam  from  the  curved  mirror  member 
of  a  two-element  grating  resonator. 

A  sample  of  the  DG  data  for  the  P(18)  (00°1)-(10°0) 
transition  is  shown  in  Fig.  5.  Figs.  6  and  7  plot  the  pressure 
dependence  of  the  hole  width  (half-width  at  half-maximum) 
in  the  DG  aud  of  the  ratio  8  between  the  line  center 
amplitude  of  the  hole  component  to  the  amplitude  of  the 
Doppler  component  in  the  DG  under  the  following  condi¬ 
tions:  7.  =  1  W/cm*,  Ir  =  20  mW/cm*,  and  amplifier 
current  of  30  mA.  The  nonzero  intercept  of  the  hole  width 
plot  at  zero  pressure  is  consistent  with  the  observed 
frequency  stability  of  the  two  oscillators*  and  it  is  clear 

*  Observation  of  the  beat  frequency  signal  between  the  two 
lasers  generated  in  a  ue:Au  detector  showed  a  short-term  jitter  of 
approximately  1  MHz. 


/f 


from  Fig.  7  that  individual  rotation-vibration  transitions 
assume  greater  inhomogeneous  character  with  decreasing 
pressure. 

We  will  now  present  a  rate  equation  formulation  of 
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collision  effects  in  molecular  rotation-vibration  transi¬ 
tions.  Our  goal  is  to  obtain  an  expression  for  the  optically 
induced  change  in  the  frequency  profile  of  the  inversion. 
The  theory  will  predict  the  observed  pressure  dependence 
of  the  saturation  behavior  and  will  elucidate  the  effect 
of  collisions  towards  Lamb  dip  formation  in  molecular 
lasers. 


Rate  Equations 


An  exact  treatment  of  the  collision  problem  for  COi 
requires  the  inclusion  of  both  rotational  and  vibrational 
energy  exchange  collisions.  To  simplify  the  analysis  we 
will  only  treat  the  more  rapid  and  therefore  more  impor¬ 
tant  rotational  coupling  process  explicitly.  This  is  done  by 
assuming  a  single  rotation-vibration  manifold  of  levels 
for  the  upper  and  lower  states.  A  straightforward  extrapo¬ 
lation  of  the  results  of  the  simple  model  will  account  for 
the  effects  of  the  slower  coupling  to  other  vibrational  levels. 

The  dynamic  processes  that  will  affect  the  saturation 
response  of  a  single  transition  within  the  rotation-vibration 
spectrum  are  shown  in  Fig.  8.  Denote  the  vibration  and 
rotation  indices  by  the  pair  (a,  j.)  where  a  =  2  for  the 
upper  vibrational  level  and  a  =  1  for  the  lower  vibrational 
level.  For  a  single  linearly  polarized  field  we  can  write 
the  following  equations  for  the  population  distribution  of 
the  level  pair  being  stimulated. 


f  [n.,»]  =  w.,.(v)  -  t n.tM  ±  «(w,  v) 

-  ^nWl(t>)J  _  [/__*'  <$.,(»,  v')nml.(v) 

-  f  dv'  S..(v',  u)n.,.(u')J 

-  £  T  [  dv '  S.,(j.v,  j'.v')n.i.(v) 

L*'-- 

-  f  dv'  S.rijW,  •  (1) 


For  clarity  the  equations  have  been  written  so  that 
collision  processes  are  expressed  in  terms  of  their  respective 
rates.  In  (1),  v  is  the  velocity  component  along  the  field 
propagation  direction,  unprimed  rotational  indices  repre¬ 
sent  the  two  levels  connected  by  the  field,  and  primed 
rotational  indices  represent  all  other  levels.  The  pump 
rates  wti.(v)  are  assumed  to  equal  waij(v)  where  /(»)  — 
(Vw  «)”*  exp  —  (v/u)2  and  \mu  =  kT.  Equations  for 
rotational  levels  not  connected  by  the  field  are  obtained 
from  (1)  by  setting  R(u,  v)  to  zero. 

The  stimulated  emission  rate  R(w,  v)  can  be  obtained 
from  the  rate  equation  approximation  of  the  density 
matrix  equations  modified  to  include  molecular  motion. 
For  a  traveling  wave  field 


Et(r,  0  -  Re  !*#  exP  i(±K-r  ~  «0) 

R(u,  v)  is  given  by  the  expression 

R(o>,  v)  -  ^  (IJe'LAw  -  «o  =F  Kv).  (2a) 


Fig.  7.  Ratio  of  the  line  center  amplitude  of  the  hole  component  to 
the  Doppler  component  in  the  DO. 


8.  Processes  affecting  saturation.  Only  the  two  rotational 
levels  connected  by  the  field  are  shown.  1 — lower  vibrational 
quantum  number,  S — upper  vibrational  quantum  number, 
) — rotational  quantum  number,  w— pump  rate,  y — vibrational 
decay  rate,  St — elastic  velocity— changing  scattering  rate, 
Sr—  rotational  level-changing  scattering  rate,  R — radiative 
transition  rate. 


For  a  standing  wave  field 

E(r,  t )  -  £+(r,  0  +  E-(r,  0, 

R(u,  v)  becomes 

*(«,  v)  -  £  (|)V[L+(«  -  «„  -  Kv) 

+  L.(«  -  wo  +  Kv)].  (2b) 

In  (2)  pJ  =  (2),  +  l)"1  |(1;,  \p-e\  2;a)|5,  w0  is  the  line 
center  frequency.  L.  (x)  =  [x2  +  (y')]’1,  and  /  is  the  total 
phase  interruption  rate  whose  component  contributions 
we  will  not  separate.  Two  physical  processes  have  not 
been  included  in  (1)  and  (2).  First,  spatial  diffusion  of 
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excited  states  through  the  finite  extent  of  a  laser  beam  at 
low  pressures  competes  with  other  vibrational  relaxation 
processes  [10].  The  rate  equations  with  the  diffusion  terms 
included,  however,  are  not  analytically  soluble.  For  our 
purposes  it  is  adequate  to  realize  that  the  effect  of  diffusion 
is  to  decrease  the  available  time  for  interaction  with  the 
field  and  can  therefore  be  treated  in  an  approximate 
manner  by  augmenting  the  vibrational  decay  rate  y. 
Second,  it  is  known  that  collisions  can  cause  shifts  and 
asymmetries  in  the  molecule-field  interaction  while  R («,  v) 
is  devoid  of  such  properties.  Since  our  primary  interest 
does  not  encompass  such  effects,  which  are  expected  to  be 
small  anyway,  they  will  be  neglected  in  the  following 
analysis. 

We  assume  the  following  elastic  and  rotational  collision 
models 


S,(v,  v')  =  S.f(v)  (3) 

Sr(jv,  jV)  =  x,Srj(v).  (4) 

Xi  is  the  rotational  Boltzmann  factor  given  by 

g,  exp  ~[E(j)/kT]/ 23  Qi  exp  ~[E(j)/kT], 

i 

The  elastic  model  requires  the  molecule  to  scatter  into  a 
Boltzmann  distribution  of  final  velocities  regardless  of  its 
initial  velocity  while  the  rotational  scattering  model  in 
addition  ignores  A;  selection  rules.  The  latter  assumption 
is  consistent  with  the  results  of  rotational  relaxation 
measurements  in  the  pertinent  levels  of  CO*  [11].  For  the 
above  models  the  elastic  and  rotational  scattering  terms 
in  the  rate  equations,  denoted  by  n,(»)|,  and  n,(u)|,, 
respectively,  become 

=  -£.[«,(«)  -  f(v)  f  "  dv'  (5a) 

=  -&[(1  -  xt)n,(v)  -  x,f(v) 

■  £  [  dv'  n,.(f/)l  •  (5b) 

To  obtain  an  explicit  equation  for  the  inversion 
An(a)  =  nlt,(v)  —  nUl(v),  we  make  the  following  simpli¬ 
fying  assumptions:  SIt  =  S,„  Si,  =  St„  and  yt  =  yt. 
These  approximations  are  conveniences  that  do  not  affect 
the  manner  in  which  collisions  influence  saturation  but 
do  simplify  the  manipulation  of  the  equations.  An  addi¬ 
tional  simplification  can  be  obtained  by  noticing  that 
under  typical  CO.  laser  plasma  conditions  (T  =  400°K) 
the  ratio  g3/gi  is  near  unity  for  transitions  near  the  peak 
of  the  distribution  (g^/gi  =  0.95  for  a  P(18)  transition). 
We  therefore  set  the  ratio  to  unity. 

Using  the  above  approximations  and  the  following 
definitions 


Nj(v)  =  A  n(v) 

Wj(t>)  =  w3i,(v)  -  wu,(v) 

in  which  J  represents  the  pair  of  rotational  quantum 


numbers  (ji,  h)  On  h  ±  •),  (3)  and  (5)  combine  to 
give  the  follovdng  equation  for  the  inversion  density  N  j(v) : 

Nj(v)  =  Wj(v)  -  yNj(v)  -  v)Nj(v ) 

-  -  f(v)  /“  dv'  A^')] 

-  Sr[(l  -  Xj)Nj(v)  -  xjKv)-Z  J  dv’  (6) 

To  obtain  the  steady-state  solution  Nj(v),  we  must 
evaluate  the  sum  in  (6)  over  all  levels  not  directly  coupled 
to  the  field.  This  is  done  in  Appendix  I.  The  result  is 

a(*-W/0)  +  M/(0  f  dv'  NA» ') 

NAv)  =  - (7) 

1  +  wJR(u> v) 

The  quantities  y„,  y,,  and  A  are  defined  by  the  following 
relations 


y,r  =  s,  +  -1 — ST 

(8) 

1  +  XjS , 

y>  =  y  +  S,  -f  (1  —  xj)Sr 

(9) 

A  -  1  +  1  -1  “  *' 

(10) 

XJ  1+  y 
^  XjSr 

y,  -  7.r  =  Ay. 

(11) 

y„  is  the  collisional  cross-relaxation  rate. 

Before  we  proceed  to  a  detailed  consideration  of  (7),  a 
very  important  inf*_„nce  of  rotational  coupling  on  the 
saturation  process  can  be  seen  by  disregarding  molecular 
motion.  Setting  v  =  0,  dv'  Nj(v')  =  Nj  and  using 
the  line  center  value  of  R(u>)  in  (7),  we  obtain 

N°, 

Nj - (12) 

1  +  ~AE\ 

The  quantity  El  =  yy'ih/p)*  is  the  saturation  parameter 
for  the  J  transition  calculated  assuming  two  isolated 
levels.  The  A -fold  increase  in  the  saturation  parameter 
over  the  value  for  an  isolated  transition  results  from  the 
fact  that  collisional  coupling  among  the  rotational  levels 
funnels  the  available  inversion  in  adjacent  rotational 
levels  into  the  pair  being  stimulated  [12].  In  the  limit 
y/xjS,  «  1  the  factor  A  approaches  (l/xj),  which  is  15.5 
for  a  P(18)  transition  at  T  =  400°K. 

Returning  to  (7)  we  can  obtain  an  exact  solution  for 
N j(v)  by  integrating  the  equation  to  extract  the  total 
inversion  density  dv'  N j(v').  This  solution,  however, 
involves  the  function  Jr.  dv  /(u)[i  +  (2 /y,)R(u,  t>)]-1, 
which  cannot  be  analytically  evaluated.  To  obtain 
analytically  tractable  expressions,  we  confine  our  attention 
to  the  regime  of  strong  Doppler  broadening  (y’/Ku)  «  1 
and  weak  saturation  (2/y, )R(u,  v)  «  1.  Equation  (7)  is 
evaluated  in  Appendix  II  in  the  above  limits.  We  find  that 
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the  differential  inversion  (Dl)  (unsaturated  inversion- 
saturated  inversion)  naturally  separates  into  two  terms 
with  different  velocity  dependences 

&N j(v,  oi )  =  6N L(v,  oi )  +  &Nd(v,  oi).  (13) 


of  the  field  and  interact  with  it  after  undergoing  n  velocity¬ 
changing  collisions,* 


Pn 


dv'  S.( v,  v') 


(18) 


The  predominently  Lorentzian  term 


Pt  and  r  become 


SNl(v,  oi)  ca 


E2N°f(v)- 


{ 


L± 

L.  +  L- 


}■ 


(14) 


which  represents  saturation  of  the  velocity  distribution 
over  a  velocity  width  2y'/K  and  is  small  compared  to 
the  full  distribution  width  2u,  is  the  well  known  Lorentzian 
hole  expected  for  a  Doppler  broadened  line.  The  additional 
Doppler  component 

«•*•*>“  is?  «p  [-feO’lfe) 

(15) 

which  indicates  saturation  across  the  entire  velocity 
distribution,  is  the  result  of  collisions  that  change  the 
molecular  velocity.  The  upper  and  lower  cases  in  braces 
of  (14)  and  (15)  correspond  to  traveling  wave  and  standing 
wave  fields,  respectively. 

The  degree  of  mixed  saturation  behavior  is  best  char¬ 
acterized  by  the  ratio  r  defined  by 


_  *Nj&iSt£. 

SND(0,oio) 


(16) 


From  (14)  and  (15)  we  have 

r=Ku  Ay 
= 


(17) 


The  ratio  r  is  thus  independent  of  field  in  the  weak 
saturation  limit.  In  the  limit  y/y„  =  <*>  the  DI  assumes 
the  pure  Lorentzian  form.  In  the  opposite  limit  y/y„  =  0 
the  DI  assumes  a  pure  Doppler  form.  Collisions  can 
therefore  restore  homogeneity  to  a  line  with  strong 
Doppler  broadening. 

To  appreciate  the  physical  significance  of  (17)  it  is 
helpful  to  derive  the  quantity  r  by  an  alternative,  intuitive 
method  for  pure  elastic  scattering  in  the  following  way. 
The  ratio  r  is  approximately  equal  to  the  ratio  Pi/P% 
where  Pt  is  the  field  interaction  probability  for  a  molecule 
whose  Doppler  shifted  emission  frequency  coincides  with 
that  of  the  field  (|(«  —  w0)  T  Kv\  <  y')  and  Pt  is  the 
interaction  probability  for  a  molecule  whose  Doppler 
shifted  emission  lies  far  from  the  field  frequency 
(|(w  —  «„)  T  Kv |  >  y').  Pi  evaluated  at «  =  «0  and  v  =  C 
is  approximately  R( o>0,  0)/[y,  +  R(o>0,  0)J.  Pt,  however, 
consists  of  an  infinite  sum  of  terms,  P,  =  p„,  in 

which  p,  is  the  probability  that  a  molecule  of  velocity  v 
will  migrate  to  the  neighborhood  |(«  —  w0)  =F  Kv\  <  y' 


~  (&)  i  Np 

-  VrKu\yJ  h\y,rl 


(19) 


Except  for  a  numerical  factor  of  ir,  the  above  result  is 
identical  to  the  elastic  limit  of  (17).  We  can  now  interpret 
r  as  follows.  y'/Ku  is  the  probability  that  a  molecule  will 
scatter  into  the  interaction  region  of  width  2 y'/K  in 
velocity  space  out  of  the  total  available  width  2u  of  the 
velocity  distribution  while  S./y  represents  the  effective 
number  of  ways  that  the  above  scattering  process  can 
occur  due  to  all  possible  collision  routes. 

The  presence  of  a  set  of  closely  spaced  rotational  levels 
adds  additional  features  to  this  picture.  The  cross-relaxa¬ 
tion  rate  must  now  be  generalized  to  include  rotational 
collisions  since  population  can  now  be  exchanged  among 
different  portions  of  the  velocity  distribution  via  rapid 
exchanges  with  other  rotational  levels.  The  form  (8)  of 
the  rotational  contribution  to  7,r  is  a  consequence  of  the 
fact  that  a  velocity  shift  in  a  given  rotational  level  induced 
by  rotational  scattering  requires  a  minimum  of  two 
collisions  ( vj  — »  v'j'  -*  v"j).  This  assertion  is  easily  demon¬ 
strated  by  considering  the  elastic  example  again.  If  we 
impose  the  condition  that  a  minimum  of  two  collisions  be 
required  to  transport  a  molecule  into  the  vicinity  of  the 
field,  then  the  summation  in  (19)  must  be  replaced  by 


The  cross-relaxation  rate  is  therefore  reduced  from  S,  by 
a  factor  in  the  denominator  which  is  analogous  to  the 
rotational  coupling  factor  (1  +  y/xjSr)~l  appearing  in 
the  rotational  term  of  (8).  Because  this  factor  is  essen¬ 
tially  unity  in  the  pressure  range  of  a  few  tenths  of  a  torr, 
and  since  rotational  thermalization  is  rapid  in  COi,  the 
rotational  scattering  process  constitutes  an  important 
contribution  to  the  spectral  cross-relaxation  or  spectral 
hole-filling  process. 

More  important  than  its  contribution  to  the  cross- 


’  The  probability  that  a  molecule  will  not  scatter  into  the  vicinity 
of  the  field  in  a  collision  is  given  by  (\[yMdv'S,(v,  »')  where  the 
proper  limits  of  integration  exclude  the  field  interaction  region.  In 
the  strong  Doppler  umit,  however,  inclusion  of  the  narrow  inter¬ 
action  region  is  of  little  consequence. 
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relaxation  process,  which  restores  homogeneity  to  the 
line,  however,  is  the  opposite  tendency  for  rotational 
scattering  to  augment  the  inhomogeneity  of  the  line  by 
the  factor  A.  This  effect  results  from  an  A -fold  increase  in 
the  vibrational  decay  rate  y  and  is  explained  by  the  fact 
that  the  effective  decay  rate  in  the  presence  of  rotational 
coupling  is  neither  the  vibrational  decay  rate  y,  nor  the 
total  rotational  scattering  rate  (1  —  Xj)Sr,  but  is  the 
net  rate  out  y,  minus  the  net  rate  in  ytt,  which  is  just  Ay 
from  (11).  Since  the  rotational  contribution  to  7„r  is 
generally  smaller  than  the  elastic  contribution,  the  net 
effect  of  strong  rotational  coupling  is  to  augment  the 
inhomogeneity  of  the  rotation-vibration  transition.  The 
extent  to  which  rotational  collision  effects  actually 
influence  saturation  is  determined  by  the  rotational 
coupling  factor  (1  +  y/xjSr)~l  (8)  and  (10).  The  tight 
coupled  and  decoupled  limits  are  given  by  y/xjS ,  =  0 
and  y/xjSt  =  ® ,  respectively. 

It  is  now  obvious  what  the  effects  of  collisions!  coupling 
among  the  various  vibrational  levels  will  be.  First  the 
factor  A  will  be  multiplied  by  the  effective  number  of 
vibrational  levels  coupled  to  the  upper  and  lower  laser 
levels  and  second  a  term  whose  magnitude  depends  on  the 
ratio  between  the  rates  for  vibration-vibration  transfer 
and  vibrational  decay  must  be  added  to  yor. 

Lamb  Dip 

The  laser  intensity  can  be  calculated  by  requiring  that 
the  saturated  gain  a(w,  w.)  be  equal  to  a  cavity  distributed 
ksas paianifetsf  J.  The  gain  n(i*,  fora  LK4.pU.r  brwJUxrjd 
line  is  «(«,  «,)  =  dv  Nj(v,  w.)  <r(w  —  Kv).  The  field 
independent  quantity  a(u  —  Kv)  -  (hu/I)R(u,  v)  is  the 
Stimulated  emission  cross  section  per  molecule. 

'(“  -  *">  -  '••{l,  +  l) 

with  (t0  =  y'(fiu/ct0)(p/h)2.  Using  (14)  and  (15)  for  a 
standing  wave  field  we  find 

«(«. «.)  ^  ao(w)  —  j 

•{l  +  \  [1  +  (y')'L(o>  -  «„)]}•  (22) 

If  t  is  the  mirror  transmittance,  the  laser  intensity 
/( w  —  w0)  =  (i)c(<0/2)Z?J<  becomes 

«»■-  «*>  -  rri^b  <23) 

where 

V  =  — ~  (24) 


L'(w  -  wo)  =  (y')*£(w  -  wo).  (25) 


The  strength  of  the  Lamb  dip  is  controlled  by  the 
parameter  y.  We  see  that  the  dip  depth  reflects  the  degree 
of  mixed  broadening  in  the  line.  This  is  perfectly  clear 
since  the  Lamb  dip  is  a  manifestation  of  the  Lorentzian 
holes  burned  into  the  inversion.  When  cross  relaxation  is 
weak  (r  »  1),  the  DI  assumes  strong  Lorentzian  character 
and  the  depth  parameter  y  approaches  1.  This  is  the  limit 
most  favorable  for  Lamb  dip  formation  and  is  closely 
approached  in  atomic  gas  lasers  where  the  short  level 
lifetimes  (7"’  ^  10'*  s)  render  the  effects  of  velocity¬ 
changing  collisions  less  important.  As  cross-relaxation 
effects  become  more  prominent  (r  «  1)  due  to  longer  level 
lifetimes,  the  DI  assumes  a  stronger  homogeneous 
character  and  the  Lamb  dip  weakens.  In  contrast  to  atomic 
levels,  the  relatively  long  vibrational  lifetimes  ( y~l  ca  10“*  s) 
of  molecular  levels  allow  collisional  processes  to  dominate 
the  saturation  behaviour.  Formation  of  the  Lamb  dip  in 
the  COa  laser  at  pressures  between  0.1  torr  and  1.0  torr 
where  it  has  been  observed  is  made  possible  through  the 
influence  of  strong  rotational  coupling,  which  effectively 
increases  the  inhomogeneous  character  of  the  transition. 

Collisional  effects  in  saturation  of  atomic  transitions 
have  also  been  treated  theoretically  by  density  matrix 
approaches  [13J  and  [14].  Our  formulation  and  that 
of  [13]  and  [14]  share  the  common  feature  of  identical 
elastic  collision  models.  Although  the  saturated  inversion 
is  not  calculated  in  the  self-consistent  field  approach  [13], 
the  dependence  of  the  Lamb  dip  factor  on  the  atomic 
collision  parameter*  given  by  [13,  eq.  (125),  [128)],  in  the 
limit,  of  tipiai  u  p,T  and.  lower  state  lifetimes,  'jgTtvn, 
exactly  with  the  dependence  predicted  by  (24)  in  the 
limit  of  no  rotational  coupling,  Sr  -  0.  A  calculation  of 
the  SeCuud  jrder  pul<irisa.tiuu  using  the  lunualisiu  of  1 15] 
also  shows  exact  agreement  with  the  elastic  limit  of  (14) 
and  (15).  The  method  adopted  in  [14]  allows  an  exact 
solution  of  the  density  matrix  equations  for  a  traveling 
wave  perturbation,  ihe  phenomenon  of  Doppler  satura¬ 
tion  again  appears. 

Comparison  of  Experimental  Result  With  Theory 

V*  idle  tilt  quantity  r  ueovirilmig  li.e  o»fciiraliv»Il 

was  calculated  by  assuming  interaction  of  the  molecules 
with  a  single  field,  a  probe  of  the  saturated  gai»'  involves 
the  simultaneous  interaction  of  two  resonant  fields.  The 
relationship  between  r  and  the  measured  gain  quantity 
8  =  6aL(u0,  w0 )/6aD(u0,  &>o)  will  therefore  be  determined 
by  the  details  of  the  nonlinear  coupling  between  the  two 
fields. 

Dienes  [15]  has  shown  that  in  the  absence  of  collisions 
both  the  magnitude  and  the  frequency  behavior  of  the 
coupling  is  dependent  upon  the  level  degeneracies.  Using 
the  theory  of  [15]  for  a  P(18)  ( J  =  18  — >  J  =  19)  transition 
and  orthogonal  fields  we  have  found  that  the  frequency 

4  The  —I  exponent  in  the  expression  for  8  appearing  in  [13, 
eq.  128b]  should  be  absent. 


1 4.0 

dependence  of  the  coupling  is  not  changed  appreciably 
from  the  classical  Lorentzian  result  but  the  amplitude  of 
the  coupling  is  reduced  from  the  classical  value.  In  terms 
of  the  inversion  quantity  &N L(v,  oi.),  the  coupling  can  be 
expressed  approximately  as 

SaL(oir,oi.)  cvc,  f  dv  &Nl(v,  ».)»(&),  —  Kv) 


where  c,  ^  0.07.  Using  (14)  we  have 


Notice  that  the  hole  in  the  DG  is  twice  the  width  of  the 
hole  in  the  DI  because  the  probe  beam,  which  interacts 
with  molecules  in  a  frequency  width  2y'  about  o>„  interacts 
with  the  hole  of  width  2y'  over  a  region  of  width  4y'.  This 
result  and  the  slope  of  Fig.  6  gives  a  value  for  the  phase 
interruption  rate  of  y'  =  12  MHz/torr,  which  agrees  with 
values  obtained  from  measurements  of  the  Lamb  dip 
width  [6J. 

Collisions,  which  result  in  a  change  of  the  DI  from  the 
form  6Nl(v,  oi)  to  5Nl(v,  oi)  +  JNc  (v,  o>),  introduce 
complications.  The  form  of  the  coupling  in  this  case  is 
difficult  to  predict  because  a  two-field  theory  including 
degeneracy  and  collision  effects  is  not  available.  We  can 
speculate,  however,  that  the  relation  between  r  and  s 
would  fall  within  the  rather  large  interval  bounded  by  the 
limits  in  which  the  magnit  Je  of  the  coupling  C2  for  the 
Doppler  component  is  either  the  same  as  for  the  Lorentzian 
component  or  is  equal  to  1  (no  coupling) . 

SaD(oir,  oi.)  cactj  dv  SND(f),  oi,)a(o\  —  Kv). 

From  (15) 

"•>  -  c’{tk){ 70  aW,  © 

We  therefore  might  expect  2s  <  r  <  2<w. 

To  compare  with  experimental  results,  we  can  deduce 
the  quantity  y  from  the  data  of  Fig.  7  assuming  r  =  2s. 
This  result  is  compared  in  Fig.  9  with  values  of  y  obtained 
by  measuring  the  depth  of  the  Lamb  dip  in  a  separate 
experiment  using  a  flow  laser  with  a  50/50  mix  of  COa 
and  Na  gases.4  There  is  good  qualitative  agreement 
between  these  two  independent  determinations  of  y,  The 
other  limit  r  =  2 as  produce  values  of  y  that  are  so  large 
as  to  fall  beyond  the  scale  in  the  figure.  Reliable  data 

*  The  authors  are  indebted  to  H.  T.  Powell  who  made  these 
meMuremente. 
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Fig.  9.  Comparison  of  the  parameter  y  as  derived  from  the  DG 
data  and  as  measured  from  the  Lamb  dip. 

above  0.5  torr  could  not  be  obtained  because  the  Lorentzian 
component  became  too  weak  to  measure  accurately. 

It  is  difficult  to  obtain  a  precise  fit  of  the  data  to  all  the 
relaxation  parameters  appearing  in  r  for  two  reasons. 
First,  the  diffusion  contribution  to  y  is  complicated  and 
depends  on  the  exact  beam  geometry,  which  in  the  case 
of  two  overlapping  beams  is  not  well  defined.  Second, 
dissociation  of  the  binary  COa  +  Na  mixture  used  in  the 
amplifier  results  in  an  average  of  the  collision  parameters 
for  the  various  species  present.  We  can,  however,  estimate 
a  value  for  yet.  The  mean  L,  beam  radius  of  2  mm  inside 
La  gives  a  value  for  7  of  7  X  103  s-1  at  0.2-torr  pressure. 
From  rotational  relaxation  measurements  [11],  ST  ~  107 
s"1  torr-1.  The  rotational  coupling  factor  1  +  y/xjS,  is 
therefore  unity  and  A  becomes  1  /x}  ^  15.5  for  J  =  18 
at  T  =  400°K.  Using  Ku  =  23  MHz,  7'  =  12  MHz/torr 
and  r  =  0.7  at  0.2  torr,  ye,  is  calculated  from  (17)  to  be 
7c,  —  9  X  104.  Since  7.,  ~  S.  +  Sr  >  Sr  2  X  10*  s“‘ 
at  0.2  torr,  the  deduced  value  of  7,,  is  at  least  a  factor 
of  2  smaller  than  expected.  This  discrepancy  is  not 
surprising  because  we  have  not  included  collisional 
coupling  between  the  various  vibrational  levels  of  COa. 
It  is  known  that  collisional  coupling  among  the  (00°ps) 
levels  [16]  and  between  (10°0)  and  the  lower  (v^O)  levels 
[17],  [18]  are  rapid.  These  levels  should  still  be  partially 
collieionally  coupled  to  the  (00°1)  and  (10°0)  levels  at 
our  operating  pressures.  The  factor  A  must  therefore  be 
multiplied  by  the  effective  number  of  vibrational  levels 
coupled  to  the  ones  being  stimulated  by  the  laser  field. 
Increasing  A  brings  the  value  of  7or  deduced  from  our 
measurements  into  closer  agreement  with  the  expected 
value. 

The  major  limitation  in  the  present  study  is  the  in¬ 
ability  to  extend  the  range  of  measurements  in  pressure 
due  to  laser  instabilities.  With  more  stable  oscillators  it 
may  be  possible  to  separate  the  contributions  of  the 
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elastic  and  rotational  contributions  to  the  cross-relaxation 
parameter  and  study  collisional  decoupling  effects 
through  changes  in  A.  The  difficulty  of  these  measure¬ 
ments,  however,  lie  in  the  fact  that  these  various  effects 
are  coupled  and  that  r  can  no  longer  be  measured  accu¬ 
rately  when  it  becomes  appreciably  greater  than  unity. 


Appendix  I 


We  first  write  the  steady-state  equations  for  levels 
J'  7^  J.  They  are 


0  =  Wj  '(v)  yNj.(v)  -  S.  Nj-(v)  -  f(v)  J  dv'  JV,<(t/)] 


-  Sr  (1  -  xJNj.fy)  -  xj.ffy)  X  f  dv’  Nk(v’)  ■ 

k*J'  J-® 


By  first  integrating  over  v,  rearranging  the  sum  over  k  and 
then  summing  over  J'  ^  J  we  obtain 


X)  Wj.  +  (1  -  Xj)Sr  f  dvNjip) 

X  [  dvNJ.(v)  = 

J‘*J  J-OD 


y  +  XjSr 


The  right  side  of  the  above  expression  can  be  expressed 
solely  in  terms  of  the  parameters  of  the  level  pair  J  if  we 
assume  that  a  Boltzmann  distribution  among  the  rota¬ 
tional  levels  does  exist  in  the  absence  of  fields.  In  this 
case  =  JVy-(p)|r  =  0  and  (3)  of  the  text  reduces 

to  v>j-1(v)  =  yNj ,f(v)  =  y(xj./xj)Nj1(v).  The  sum 
X-'w  WJ'  thus  becomes  yN°(l  —  Xj)/ij.  Substituting 
/"»>  dv  Nj-(v)  into  (6)  with  Nj(v)  -  0  and  re¬ 
arranging  terms  gives  (7). 


Appendix  II 

Integration  of  (7)  gives  the  following  expression  for  the 
total  inversion  density 

dv’  Nj(v')  =  — ^ - 

1  -  ^  F(o>) 

in  which  F(co)  =  J"„  dv  /(u)[l  +  (2/y,)R(u,  a)]”1.  Using 
the  strong  Doppler  broadening  and  weak  saturation 
assumptions  results  in  the  approximation 

F(&>)  ~  1  —  £(w). 

For  a  traveling  wave  field 

*(«)  ~(\4  Ku)~l  exp  [-  (w  -  wn/Ku)2}{n/y,){P/h)2E\ 

£(&>)  for  a  standing  wave  is  exactly  twice  the  traveling 
wave  value.  Substituting  the  above  value  of  /"„  dv'  Nj(v’) 
in  (7)  gives 

™ 


5Nj(v,  w)  a*  v)  +  (1  -  B)j 

The  quantity  B  equals 

A(y/y ,)(1  +  y„r[l  -  {(w))/[y,  -  y0,  +  Yor  {(«)]). 

Use  of  the  weak  saturation  condition  in  expanding  (1  —  B) 
gives  (15).  Since  weak  saturation  implies  small  changes 
in  inversion,  5ND(v,  to )/N°(v)  =  1  —  B  <K  1  imply  B  ~  1. 
This  result  immediately  gives  (14). 


Conclusion 

We  have  measured  the  frequency  profile  of  the  saturated 
gain  in  rotation  vibration  lines  of  the  COa  molecule  at 
pressures  where  individual  lasering  transitions  display  the 
Lamb  dip.  Each  transition  saturates  in  a  mixed  inhomo¬ 
geneous-homogeneous  broadened  manner  due  to  spectral 
cross  relaxation  induced  by  collisions  that  change  the 
molecule’s  velocity. 

A  rate  equation  treatment  shows  that  in  COa  this  mixed 
saturation  behavior,  characterized  by  the  mixed  broaden¬ 
ing  parameter  r,  which  is  the  ratio  of  the  amplitude  of  the 
Lorentzian  component  to  the  amplitude  of  the  Doppler 
component  in  the  differential  inversion,  is  a  balance 
between  the  three  collision  processes  of  elastic  scattering, 
rotational  thermalization,  and  vibrational  decay.  The 
depth  of  the  Lamb  dip  in  the  laser  intensity,  which  reflects 
the  influence  of  these  various  collision  processes,  is  given 
by  the  expression  y  =  [1  +  4/r]"1. 
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